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Abstract 

We review recent experimental and theoretical work on ultrasmall metallic grains, 
i.e. grains sufficiently small that the conduction electron energy spectrum becomes 
discrete. The discrete excitation spectrum of an individual grain can be measured by 
the technique of single-electron tunneling spectroscopy: the spectrum is extracted 
from the current- voltage characteristics of a single-electron transistor containing the 
grain as central island. We review experiments studying the influence on the discrete 
spectrum of superconductivity, nonequilibrium excitations, spin-orbit scattering and 
ferromagnetism. We also review the theoretical descriptions of these phenomena in 
ultrasmall grains, which require modifications or extensions of the standard bulk 
theories to include the effects of level discreteness. 
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1 Introduction 



One of the most fundamental features of quantum mechanics is the fact that 
the energy spectrum of a system of particles confined to a small region is 
discrete or quantized, with the typical spacing between energy levels increas- 
ing with decreasing system size. In atomic and nuclear physics, spectroscopic 
techniques for measuring and analyzing such discrete spectra have for decades 
been a major source of detailed information on the forces between the parti- 
cles and the correlations which they experience. In condensed matter physics, 
however, it has been much more difficult to spectroscopically study the dis- 
crete spectrum of an individual sample, since system sizes were typically so 
large that discrete eigenenergies could not be resolved on the energy scale set 
by the temperature. 

This changed in the course of the last 15 years due to advances in micro- 
fabrication techniques, which made it possible to study individual systems of 
mesoscopic or nanoscopic dimensions, whose characteristic length scales range 
from a few /im down to a few nm. In the early 1990's, semi-conductor de- 
vices were used to fabricate the first "quantum dots" , i.e. droplets of charge 
confined to a 2-dimensional region so small (radius of order 50 nm) that dis- 
crete levels in the conduction electron spectrum [1,2] could be resolved at 
dilution rcfridgcrator temperatures in the 10-100 mK range. The technique by 
which this was done is called single- electron-tunneling spectroscopy: the dot is 
connected to two leads via electrostatically-defined tunnel barriers to form a 
so-called single-electron transistor (SET), the current- volt age characteristics 
of which are measured and analyzed. Under certain conditions the conduc- 
tance shows well-defined resonances which can be associated with tunneling 
through discrete eigenstates of the dot. Since quantum dots exhibited various 
features familiar from atomic physics, such as energy shells that feature magic 
numbers and are filled according to Hund's rules, etc., they are often aptly 
referred to as "artificial atoms" (see e.g. the collection of reviews in [3], in 
particular that by Kouwenhoven et al. [4]). 

In the mid 1990's, a similar advance was achieved with metals [5-10], when 
Ralph, Black and Tinkham (RBT) succeeded to perform single-electron-tun- 
neling spectroscopy on individual ultrasmall metallic grains (of radii r < 5 nm 
and mean level spacings d > O.lmeV): by attaching such a grain via oxide tun- 
nel barriers to two leads they constructed a single-electron transistor having 
the grain as central island, and showed that a well-resolved, discrete excitation 
spectrum could indeed be extracted from the conductance. This opened up a 
new frontier in the study of electron correlations in metals, since the ability to 
resolve discrete energy levels allows the nature of electron correlations to be 
studied in unprecedented detail. During the last several years, single-electron- 
tunneling spectroscopy of ultrasmall metallic grains has been used to probe 
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superconducting pairing correlations in Al grains [6,10], nonequilibrium exci- 
tations [10-13] and spin-orbit interactions [5,7,14-16] in normal grains, and 
ferromagnetic correlations in Co grains [17,18]. A brief overview of some of 
the early experiments is given in Ref. [9] by Ralph et al, a more thorough one 
in the thesis of Black [8] . 

Although tunnel-spectroscopic studies of metallic grains are similar in spirit 
to those of semiconductor quantum dots, there are a number of important 
differences: 

(1) Metals have much higher densities of states than semiconductors (because 
the latter have smaller electron densities and effective masses), hence 
metals require much smaller sample sizes (< 10 nm) before discrete levels 
become resolvable. 

(2) Consequently, metallic grains have much larger charging energies than 
quantum dots, which is an advantage when fluctuations in electron num- 
ber are to be minimized. On the other hand, this also implies that for 
quantum dots the electron number can be varied over a much larger 
range than for metallic grains, which is useful for analyses of statistical 
properties. 

(3) For metallic grains, the ability to study a range of materials throughout 
the periodic table, including samples doped with impurities (which tend 
to trap or deplete electrons in quantum dots) and alloys, allows some 
control over the strength and type of electron interactions to be studied. 
In particular, one may study superconductivity and itinerant ferromag- 
netism. 

(4) For metallic grains, the tunnel barriers to the leads are insulating oxide 
layers and hence insensitive to applied voltages, whereas for quantum 
dots they are electrostatically defined and hence tunable. Tunability is 
usually an advantage, but not always: for example, nonequilibrium effects 
are easier to study quantitatively for a metallic grain than a quantum 
dot, because for the latter a large source-drain voltage lowers the tunnel 
barrier in poorly-controlled ways. 

(5) Spin effects are easily probed in metallic grains by applying a magnetic 
field and studying the Zeeman-splitting of time-reversed pairs of states, 
which for nm-scale metallic grains is a much stronger effect than that of 
the applied field on orbital properties. In contrast, in semiconductors the 
latter are dominant over spin effects. 

(6) For the same reason, spin-orbit effects are more easily studied in metallic 
grains than in quantum dots. 

A particularly interesting feature of experiments on ultrasmall metallic grains 
is that they probe the ways in which finite-size effects modify a system's 
characteristic correlations relative to their bulk properties, both due to meso- 
scopic fluctuations and via level discreteness. Such modifications arise when 
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the new energy scale characterizing the spectrum's discreteness, namely the 
single-particle mean level spacing d = 1/M{£f) ~ 1/Vol [where M{e) is the 
density of states per spin species], becomes comparable to the energy scale 
characterizing the correlations in bulk systems (such as the energy gap in su- 
perconductors or the Kondo temperature in magnetic alloys). Such quantum 
finite size effects in metals had attracted considerable attention in the past, but 
could hitherto only be studied in ensemble-averaged quantities (for reviews, 
see [19,20]). Spectroscopic studies of discrete spectra of individual grains yield 
significant new information. Let us briefiy mention the most important ex- 
amples, organized according to the section in which they will be discussed in 
detail later. 

RBT were able to determine the number parity (even or odd) of a given grain 
[Sec. 3], by studying the evolution of the discrete spectrum in an applied 
magnetic field: For an odd grain, the ground state energy was observed to 
Zeeman-split in a magnetic field iJ, as expected for a spin- 1/2 Kramers dou- 
blet; moreover, for pure Al grains, the Lande g factor extracted from the size 
of the splitting was close to the expected value of ^i?"'''^ = 2. In contrast, for an 
even grain the ground state is a nondegenerate spin singlet, and accordingly 
no ground state splitting was observed. 

Parity effects were also observed in RBT's experiments on largish (r > 5 nm) 
Al grains [6,10,21]: an even grain had a distinct spectroscopic gap (^ d) but 
an odd grain did not, which is clear evidence for the presence of supercon- 
ducting pairing correlations in these grains [Sec. 4]. The spectroscopic gap 
for even grains was driven to zero by an applied magnetic field, hence the 
paramagnetic breakdown of pairing correlations could be studied in detail. 
The corresponding theory was worked out by Braun, von Delft, Ralph and 
Tinkham [22,23]. 

In RBT's smallest grains (r < 3 nm), however, no such distinct spectroscopic 
gap could be discerned. This observation revived an old but fundamental ques- 
tion [Sec. 5]: What is the lower size limit for the existence of superconductivity 
in small grains? , and thereby stimulated a large number of theoretical inves- 
tigations [21-44]. Anderson [45] had addressed this question already in 1959, 
arguing that "superconductivity would no longer be possible" when the mean 
level spacing d becomes larger than the bulk gap, to be denoted by A, because 
then A looses its special significance as gap in an otherwise continuous spec- 
trum. RBT's new experiments stimulated a number of theoretical attempts 
to quantitatively describe the crossover from the bulk limit <^ A, where 
superconductivity is well-developed, to the fiuctuation-dominated regime of 
d ^ A, where pairing correlations survive only in the form of weak fluctua- 
tions. Describing this crossover constituted a conceptual challenge, since the 
standard grand-canonical mean-field BCS treatment of pairing correlations 
[21-28] breaks down for d > A. This challenge elicited a series of increasingly 
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sophisticated canonical treatments of pairing correlations [33-41], based on a 
simple reduced BCS-Hamiltonian for discrete energy levels, which showed that 
the crossover is completely smooth, but, interestingly, depends on the parity 
of the number of electrons on the grain, as pointed out by von Delft et al. [21]. 
Very recently, the main conclusions of these works were confirmed [39] using 
an exact solution of the reduced BCS model, discovered by Richardson in the 
context of nuclear physics in the 1960s [46-54]. (The existence of this solution 
came as a surprise - in the form of a polite letter from its inventor - to those 
involved with ultrasmall grains, since hitherto it had apparently completely 
escaped the attention of the condensed-matter community.) 

An interesting finite-size effect was also revealed in nonequilibrium grains 
[Sec. 6], for which BRT observed the excitation spectrum to consist of dusters 
of resonances, with the spacing between clusters comparable to free-electron- 
estimates of the single-particle mean level spacing d, but the spacings between 
subresonances of the same cluster much smaller than d. A theory for this effect 
was developed by Agam et al. [11,13], who showed that clusters of resonances 
can be caused by nonequilibrium excitations on the grain, provided that meso- 
scopic fluctuations of the matrix elements of the electron-electron interaction 
are sufficiently strong. Such fluctuations are neglected in the so-called "or- 
thodox model" that is commonly used to describe single-electron transistors, 
but are expected, according to the general theory of disordered interacting 
electron systems (e.g. [13]), to become significant for sufficiently small grains, 
and can be described using random matrix theory. The experimental observa- 
tion of such clusters thus constitutes a beautiful and direct illustration of the 
importance of mesoscopic fluctuations in ultrasmall grains. 

A further example of such mesoscopic fluctuations was published by Salinas 
et al. [14], who experimentally studied the effect of the spin-orbit interaction 
[Sec. 7] in Al grains doped with Au, and by Davidovic and Tinkham [15,16], 
who studied Au grains. These authors observed effective g'^^ factors signifl- 
cantly smaller than g^^"^*^ = 2, as expected, since the spin-orbit interaction 
mixes states with opposite spin. Interestingly, the measured g^^ factors were 
also found to vary from one Kramers doublet to the next, which is a clear 
signature for mesoscopic fluctuations. Very recently, a theory for the statistics 
of the fluctuations of g'^^ was worked out by Matveev, Glazman and Larkin 
[55] and by Brouwer, Waintal and Halperin [56]. 

During the last year, Gueron et al. [18] published the flrst set of detailed exper- 
imental results on ferromagnetic Co grains [Sec. 8]. Among the novel features 
found for these are a strong asymmetry between the tunneling probabilities 
for spin-up or spin-down electrons, a larger-than-expected density of low-lying 
excitations, and hysteretic behavior of the excitation spectrum as function of 
an applied magnetic fleld. The properties of the hysteresis loop reflect the 
changes in the direction of the grain's magnetic moment as the applied field is 
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ramped. The latter features hold great promise for being a potential tool with 
which to study the dynamics of magnetization reversal in individual nm-scale 
ferromagnets. No detailed theory for these experiments exists at present. 

An ultrasmall grain containing a single magnetic impurity can in principle be 
used to study the effect of level discreteness on the Kondo effect [Sec. 9]. A 
theory for such a "Kondo box" was worked out by Thimm, Kroha and von 
Delft [57]. They found that when the mean level spacing in the grain becomes 
larger than the Kondo temperature, the Kondo resonance is strongly affected 
in a way that depends on the parity of the number of electrons on the grain, 
and that should be detectable by tunnel-spectroscopic measurements. 

The present review summarizes the above developments. The guiding principle 
in the choice of topics was to focus on the experimental studies of discrete spec- 
tra in ultrasmall metallic grains that have been carried out since 1995 using 
single-electron-tunneling spectroscopy, and on those theoretical developments 
that were directly inspired by them. The notation to be used throughout is 
introduced in Sec. 2, which describes how tunneling spectroscopy works in 
practice and reviews the well-known orthodox theory for single-electron tran- 
sistors (also reviewed in. e.g., Refs. [58-60]), emphasizing those features that 
have special relevance for ultrasmall grains. The contents of Sees. 3 to 9, which 
for the most part can be read independently of each other, were already out- 
lined above. The main conclusions of each section are summarized concisely 
in Sec. 10, which also gives an outlook towards directions for further work. 

For each subject, an attempt was made to give a detailed account of the 
arguments and calculations that have a direct bearing on understanding ex- 
perimental data, and to qualitatively explain the theoretical ideas required 
for their interpretation. Theoretical developments beyond those of direct rele- 
vance to experiment are usually either summarized or only briefly mentioned, 
but seldom reproduced in detail. The fact that more than a third of the re- 
view is devoted to superconductivity in ultrasmall grains is a reflection of the 
number of papers that have appeared on this subject since 1995. 
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2 Single-electron-tunneling spectroscopy 

The experiments in which Ralph, Black and Tinkham (RBT) succeeded to re- 
solve discrete energy levels of an ultrasmall metaUic grain used the technique of 
single-electron-tunneling spectroscopy. In Sec. 2.1 we explain the idea behind 
this technique and illustrate its capabilities by showing some representative 
data sets (without delving into the interesting physics contained in this data, 
which will be done in later sections). In Sec. 2.2 we estimate theoretically un- 
der which conditions discrete states should be resolvable, in Sec. 2.3 develop a 
theoretical description for ultrasmall single-electron transistors (SETs), and in 
Sec. 2.4 discuss experimental details concerning fabrication and measurement 
techniques. 



2.1 Ultrasmall single- electron transistor 

In the first generation of experiments of 1995 [5-7], a grain made from Al (a su- 
perconducting material) was connected to two metal leads via high-resistance 
tunnel junctions, with capacitances Cl and C^, say. In the next generation of 
1997 [10], the grain was also coupled capacatively to a gate, with capacitance 
Cg. The resulting device, schematically depicted in Fig. 1(a), has the struc- 
ture of a SET, with the grain as central island. The circuit diagram for an 
SET is shown in Fig. 1(b). Applying a bias voltage V between the two leads 
causes a tunnel current / to flow between the leads through the grain, via in- 
coherent sequential tunneling through the tunnel junctions. The CTirrcnt can 
be influenced by changing the gate voltage Vg (hence the name "transistor" ) . 
which tunes the electrostatic potential on the grain and thereby also its aver- 
age number of electrons N. (For devices without a gate these two quantities 
cannot be tuned and instead have some sample-dependent, fixed value. For 
such devices, set Cg = in all formulas below.) 

The physics of SETs had been clarified in the early 1990s [3] through extensive 
studies of lithographically defined SETs of mesoscopic size, i.e. with micron- 
scale central islands. The fundamentally new aspect of RBT's work was that 
their SETs, made by a novel fabrication technique (described in Sec. 2.4), were 
nanoscopic in size: they had ultrasmall grains with radii between 15nm and 
2nm as central islands, which were thus several orders of magnitude smaller in 
volume than in previous experiments. This had two important consequences: 

(1) The grain's charging energy Eq = e^/2C was much larger than for meso- 
scopic SETs, ranging roughly between 5 and 50 meV (where C = Cl + 
C'R + Cg). Ec is the scale that determines the energy cost for changing N 
by one. Since for ultrasmall grains it far exceeds all other typical energy 
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Fig. 1. (a) Schematic cross section of the ultrasmall SETs studied by RBT in [10], 
and (b) the corresponding circuit diagram. 

scales of the SET, such as those set by the bias voltage {V < 1 mV), 
the temperature (T < 4.2 K) and the bulk superconducting gap for Al 
(Abuik = 0.18 meV), fluctuations in electron number are strongly sup- 
pressed. 

(2) Discrete eigenstates of the conduction electron energy spectrum became 
resolvable - their mean level spacing d ranged from 0.02 to 0.3 meV, which 
is in order-of-magnitude agreement with the free-electron estimate of c? = 
27r^^^/(m/cFVol) for the single-particle level spacing. Such d-values are 
much larger than k-^T for the lowest temperatures attained (around T ~ 
30mK), but on the order of Abuik- However, the number of conduction 
electrons for grains of this size is still rather large (between 10^ and 10^). 

Since the two scales Eq and d differ by at least an order of magnitude, they 
manifest themselves in two distinct and easily separable ways in the low- 
temperature I-V curves of the devices, shown in Fig. 2 for a series of different 
values: 

(1) When V is varied on a large scale of tens of mV [Fig. 2], the I-V curves 
have a typical "Coulomb- staircase form" characteristic of SETs: zero cur- 
rent at low \V\ (the "Coulomb blockade" regime), sloping steps equally 
spaced in V, and step thresholds sensitive to Vg. This proves that the 
tunnel current flows only through one grain. The maximal width of the 
flat step of zero current around |eV^| = is governed, in order of magni- 
tude, by Ec and typically varies between 5 and 50 mV. As Vg is varied, 
the I-V curves periodically repeat, with a period e/Cg. 

(2) When V is varied on the much smaller scale of a few mV around the 
threshold of the Coulomb blockade regime and the temperature is suffi- 
ciently low (T d), the I-V curves have a step-like substructure, shown 
in Fig. 3(b) [see also Fig. 5(b)]. As first pointed out by Averin and Ko- 
rotkov [61], such small steps in the I-V curve are expected to occur when- 
ever the voltage drop across one of the tunnel junctions equals the thresh- 
old energy at which the rate for tunneling across that junction into or out 
of one of the grain's discrete energy eigenstates becomes nonzero, since 
this opens up another channel for carrying current across that junction. 
Correspondingly, the differential conductance (d//dV^) curves contain a 



10 



-60 -40 -20 20 40 60 

V (mv) 

Fig. 2. Current-voltage curves for an ultrasmall SET [10] at 50 mK, artificially offset 
on tfie vertical axis for a set of equally-spaced values of Vg between —1.2 and 1.8 V. 
The I-V curves display Coulomb-staircase structure on a bias voltage scale of tens of 
mV. By fitting these to standard SET theory [62], the SET's basic parameters can be 
determined: i?L = 3.5 MO, = 0.2 Un, = 3.5 aF, Cr = 9.4 aF, Cg = 0.09 aF, 
Eq = 46 meV. The grain radius and mean level spacing are estimated as r ~ 4.5 nm 
and d ~ 0.45 meV, using assumptions stated in Sec. 2.4.2. 

series of fine peaks, see Fig. 3(b). Under certain conditions described in 
Sec. 2.3.4(a), the distances between these peaks directly reflect the en- 
ergy differences between cigcncnergies of same-A^ eigenstates of the grain. 
Such conductance curves thus directly yield the grain's fixed-N excitation 
spectrum, or more precisely, the set of energy differences 

for those eigenstates | a) at of the A^-electron grain that are accessible final 
states for a tunneling process that removes or adds an electron to the grain 
if its initial ground state has (iV + 1) or (A^ — 1) electrons, respectively. 

The magnetic-field dependence of the fixed- excitation spectrum can be 
obtained by simply tracing the motion of the conductance peak positions as 
a magnetic field is turned on (at fixed V^g). Examples are given in Fig. 8 in 
Sec. 3.1, Fig. 10 in Sec. 4.1, Fig. 31 in Sec. 7.1 and Fig. 34 in Sec. 8.1. If a 
gate is present, two further very interesting options exist: Firstly, by tuning Vg 
by an amount large enough (~ Ec/e) to change N by one unit, the influence 
on the spectrum, of the parity of the number of electrons on the grain can be 
studied. Secondly, by tuning Vg such that the Coulomb blockade regime is 
large or small, so that the V-threshold at which current begins to fiow is large 
or small, nonequilibrium effects can be maximized or minimized, respectively, 
depending on whether one chooses to study them or not. 
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Fig. 3. (a) Cartoon of an d// dy- measurement: The number of available transport 
channels through discrete states (here three) is determined by the bias voltage V. 
(b) The current and differential conductance as functions of bias voltage V for 
one of RBT's ultrasmall grains; beyond the Coulomb-blockade threshold at around 
5.5 mV, the current displays fine steps and the conductance fine peaks, on a voltage 
scale of a few mV, reflecting the grain's discrete eigenspectrum. 

2.2 Conditions under which discrete states are resolvable 



Having shown in Fig. 3 an example of the data by which RBT demonstrated 
that discrete energy levels of an ultrasmall metallic grain can be resolved, let 
us take a step back and recall the theoretical arguments, mainly due to Averin 
and Korotkov [61], for why this should be possible. This requires comparing 
various characteristic energy scales of a 3D grain and their dependence on its 
size, such as the mean level spacing d, the charging energy Eq, the Thouless 
energy £'Thouiess, and the amount of level broadening due to inelastic relaxation 
processes (with rate F™^') and tunnehng (with rate Ftun)- 

The mean level spacing d = 1/A/'(£f) can be estimated using the free-electron 
expression for the density of states at the Fermi surface of a 3D grain: 



27v^ff _ l.SOeV.nm^ 
~ mfcpVol ~ A;fVo1 ' 

where, for example, kp — 17.5 nm~^ for Al and kp — 12.1 nm"~^ for Au. For the 
sake of order-of-magnitude estimates below, we shall use kp — 10 nm~^ and 
crudely assume the grain to have linear dimension r, volume and contact 
area (via tunnel junctions) with each lead. 

Next we consider the charging energy Ec = e^/2C. The capacitances Cl 
and C^i to the left and right leads are each of order Cnr^, where Cn is the 
capacitance per unit area of the tunnel junctions, typically of order Cn — 
0.05aF/nm^ [63]; if a gate is present, Cg is typically an order of magnitude 
smaller. We therefore have 
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O.SeV.nm^ 



, • (3) 



Although the mean level spacing increases faster (r~^) with decreasing size 
than the charging energy they become comparable only for grains of 

almost atomic size (r ~ 2 A). For nm-scale grains (say r = 15 - 2 nm), d is 
still significantly smaller (0.04 - 20 meV) than Ec (4 - 200 meV). Moreover, 
the total number of conduction electrons oc ep/d is still very large (~ 
300 000 - 500, for ep ~ 10 eV). As long as r is substantially larger than 
the Thomas- Fermi screening length Atf — I/^f? nm-scale grains can still be 
regarded as metallic. 

Since d -C Ec for nm-scale SETs, the large-scale features of their current- 
voltage characteristics [e.g. Fig. 2] can be correctly described [61] using contin- 
uous spectra for both island and leads and adopting the standard "orthodox" 
theory [58] for Coulomb-blockade phenomena in metallic SET's. For additional 
fine-structure due to a discrete spectrum to be resolvable, three more condi- 
tions must be met: 

(i) To avoid thermal smearing, T <^ d is required. 

(ii) r^, the tunneling rate out of a given discrete state \a) on the grain into 
lead r (=L,R), must be small enough that the tunneling- induced level widths 
hT"^ don't cause neighboring levels to overlap, i.e. hT"^ -C d] but this condition 
is equivalent to one which must hold anyway for Coulomb blockade effects to 
occur [59,60], namely Gr <^ 2e'^/h, where G,,. = e'^Af{ep)Tl^^ is, roughly, the 
total conductance across junction r, and F^^^^ the average of the FJ^'s near ep. 

(iii) For an excited state on the grain with energy e above the ground state 
to be resolvable, its inelastic relaxation rate F™^^ must be small enough that 
the corresponding line-width is less than the level spacing, /iF^°'^' < d. The 
question when this is satisfied depends (neglecting phonons) on the grain's 
so-called "dimensionless conductance" gdim, defined for finite systems by [64] 

ddim = -E'Thouless/ d, (4) 

^-DdiffA^ = (0.25 fcp/tr/'^^) meV.nm^ (diffusive) , 

^Thouless — \ (5) 

hvp/{a2r) = (38 kp / a r) meV .nm^ (ballistic) . 

Here the Thouless energy -Exhouiess is the inverse time for an electron near ep to 
travel once across the system, D^is = 'f^F^tr/3 is the diffusion constant, Itr the 
transport mean free path, and a a geometrical constant. Eqs. (5) are order-of- 
magnitude estimates[^ for diffusive or ballistic 3D grains, for which gdim oc r 



The expressions for g^im and -Exhouiess given by various authors differ by factors 
of order unity. (For i^dim we followed [64], and for -^Thouless [55] for the diffusive and 
[65,15] for the ballistic case.) Note that the ballistic definition for -^Thouless is similar 
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or r^, respectively. It was shown in Refs. [66,64] that for grains so smallQthat 
9 dim ^ 1, the spectrum of low-lying excitations consists of 5-function like peaks 
with ^ir^^^i < d that can definitely be resolved [64], whereas for larger grains 
with (yf ^ 1, the condition /tF™*^' < d can be met as long as e < -EThouiess- For 
grains with d <^ -Exhouiess -C Eq, the fine-structure due to level discreteness 
will thus be smeared out for e ~ Eq, but should be clearly resolvable for e < 
several units of d. (For a more detailed discussion of F™''^ including the effect 
of phonons, see Sec. 6.3.) 

For future reference, we remark here that the parameter ^fdim also controls the 
importance or not of corrections to the orthodox model, which can be shown 
[67-71,13] to be small as l/f/dim [see Sec. 6.1.3]. 

Finally, note that the amount of nonequilibrium effects on the grain is con- 
trolled by the ratio of F[y^ to F™'^'. If F^^^ ^ F™*^', nonequilibrium effects are 
negligible, since an electron tunneling into an excited state of the grain has 
plenty of time to relax before the next tunneling event. Conversely, nonequi- 
hbrium effects should become very important for F[^^ ^ F™'^'. Indeed, they 
have been observed directly [10-13] and are discussed in Sec. 6. 

2.3 Theoretical description of an ultrasmall SET 

In this section we set up a general formalism for describing transport through 
ultrasmall metallic grains, in order to show explicitly how level discreteness 
causes steps in the tunneling current and to clarify precisely what information 
can be extracted from the latter. Since ultrasmall grains may, due to their 
discrete states, be regarded as "metallic quantum dots", we can carry over 

to the diffusive one, except that the grain size determines the effective mean free 
path (/tr 3r/2a). The proper choice of a depends on the precise geometry and 
the amount of disorder present [13]; in practice, it is somewhat of a "fudge factor". 
For a hemispherical grain, Davidovic and Tinkham [15] used a = 3 (the number of 
dimensions); they argue that if the electronic motion in the grain is ballistic and 
surface scattering diffusive, an electron's mean-free path is approximately 2r, but to 
diffuse through the entire particle volume, it should scatter from the surface about 
three times [15]. For a pancake-shaped grain with thickness z, radius r and diffusive 
dynamics, Agam et al. [11] argued that a oc (r/z)'^. 

^ For a ballistic, hemispherical Al grain, for example, Eqs. (2) to (5) yield g^ixa = 
(16/o)(r/nm)^, which, for a = 3, gives ^dim < 1 if r < 0.4 nm. An experimental 
estimate [11] of ^dim for two ballistic Al grains with Vol ~ 40 and 100 nm^ yielded 
values of ^dim ~ 5 in both cases, as discussed in Sec. 6.1.3. If these grains are 
assumed to be hemispherical, this would imply very large a factors of ~ 24 and 
42; this might imply that an assumption of pancake-shaped grains would be more 
appropriate, for which a does contain a large factor {r/z)"^, cf. footnote 3 and, in 
Sec. 6.1.3, footnote 18. 
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much of the formaUsm developed in the hterature for quantum dots; indeed, 
in spirit and notation we shall very closely follow a review by Schoeller entitled 
"Transport Theory of Interacting Quantum Dots" [72] . Although this material 
may be familiar to many readers, it is needed here to establish the notation 
to be used throughout this review. 

2.3.1 Charging energy and Coulomb blockade 

We start by deriving the charging energy for the grain, using the so-called "or- 
thodox" or "Coulomb-blockade" model [58], which assumes the electrostatic 
potential to be homogeneous on the grain. This assumption was argued to be 
reasonable as long as the Thomas- Fermi screening length Atf ^ I/Zcf is very 
much shorter than the grain's linear dimensions; a more rigorous condition 
is gdim ^ 1 [67-71,13]. Nevertheless, deviations from the predictions of the 
orthodox model can be expected in principle, and indeed were observed in 
some of RBT's smallest grains with gidim — 5 (see Sec. 6). 

Consider the SET shown in Fig. 1(b), and let £'pot(-/Vex) denote the electrostatic 
work required to add A'cx excess electrons with a total charge of Qcx = ^N^x 
(with e < 0) to a grain with initial random off-set charge Qo, while the time- 
independent voltages Vl, and Vg of the left and right leads and the gate 
electrode, respectively, are held fixed. Within the orthodox model, Epot{Nex) — 
lQr'^'^° ^Q^iQ)^ where V{Q) is the electrostatic potential of the grain for 
given grain charge Q. It is determined by Cr[Vr — V{Q)] = Qr, where Qr is 
the screening charge on capacitor r (= L,R,g) and Q-^ + + Qg = —Q- 
(In the absence of a gate one has Qg = and hence should set Cg = 0.) 
Using the definitions C = Ci^ + Cr -|- Cg and = I]r=L,R,g C'rK, we obtain 
V{Q) = {Qb + Q) /C and thus 

^pot(Ac.) = gex(gD + Qo)/C + QIJ2C = e^DiVex + ^^cA^ex ■ (6) 

Since the first term on the r.h.s. is linear in the number of excess electrons, 
= (fc + Qo)/C can be viewed as the electrostatic potential on the grain 
(the subscript D, for "dot", is used to conform to Schoeller's notation [72]). 
For gated devices, Qo may be absorbed into by shifting Vg by an amount 
—Qo/Cg, but for gateless ones, it must be treated as a fit parameter. The 
second term on the r.h.s. of Eq. (6) represents the Coulomb interaction energy 
of the iVex excess electrons due to their mutual repulsion, and its scale is set 
by the charging energy Ec = e^/2C. 

Epot is often written as Ec{Ncx — ""d)^ ~ -E-c^d' where = —eVY)/2Ec, since 
this makes it clear that for a given value of n^, the system will adjust A^^cx to 
be the integer closest to n^, in order to minimize Epot- The particle number 
on the grain can thus be controlled in discrete units by varying via the 
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gate voltage Vg. Since the electrostatic energy difference 



5£;±t(A^ex) = £^pot(iVex ± 1) - £;pot(iVex) (7) 

between grains with Ncx ± 1 or A^^cx excess electrons vanishes when is tuned 
to lie half-way between A^ex ± 1 and iVex, half-integer values of n^, are called 
"degeneracy points" . Transport is possible through a grain whose ground state 
has A^ex electrons only if, roughly speaking, mm[SEp^^{Nex)] ^ maxf/cBT", 
i.e. only if the "Coulomb barrier" presented by Eq. (7) can be overcome by the 
temperature or bias voltage. If both of these are small (fceT, \eV\ <^ Eq), a 
complete suppression of transport through the grain, the so-called "Coulomb 
blockade", occurs far away from the degeneracy points; in particular, the 
grain's low-temperature, linear-response conductance shows "Coulomb oscil- 
lations" as function of V^, i.e. a scries of peaks, with a uniform n^-spacing 
of 1, i.e. V^-spacing of e/Cg. In between degeneracy points, low-temperature 
transport is possible only with a large bias voltage, which will, in general, lead 
to nonequilibrium effects. 



2.3.2 General Hamiltonian 

To describe transport through the grain, we shall adopt a Hamiltonian H 
Hi, + Hk + -f^D + -f^tun of the following rather general form: 



Hr = Y^i^kr + eV;.)4arCfcar {t = L, R) , (8) 

kcr 

^D = E(e^DA^ex,a + ^a)|«)(«| , (9) 

a 

^tun= Y'^kla 4arClcr-D + (h-C.) , (10) 

r=L,R kla 

= E E ^L,aa'cLl«)(«1 + (h-C), (11) 

r=L,R k(T,aa' 

TL,aa'-YTklA^\ciaBW) ■ (12) 
I 

Here c^^^ creates an electron in lead r, in a single-particle state with spin a, 
kinetic energy Skr (measured relative to the Fermi energy e^p^ of that lead) and 
electrostatic energy eV^. \a) is a many-body eigenstate of the isolated grain 
in the absence of tunneling, with definite total electron number N^- The sum 
on a in Eq. (9) is not restricted to a fixcd-A'^a Hilbert space, but is over all 
eigenstates in the isolated grain's Fock space. The eigenvalue eVj;)Ncx,a + 
of \a) has been split into two parts: the first contains, via Vj), all dependences 
on the external voltages arising from Epot of Eq. (6), while Sa is the voltage- 
independent remainder and includes the interaction contribution (EcN^^^ ^) to 
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Epot. In the absence of many-body correlations, e.g. for a grain made from a 
normal metal, one could take, as Averin and Korotkov did [61], 

^normal ^ Epot{Ne.) + E ^/^dcLdQ^^D, (13) 

la 

l«) = \{n^l}) , = EcN^^ + E si^^ ntl ■ (14) 

Jo- 
Here cj^j^ is a creation operator for one of the grain's single-particle states with 
energy ejo-o, and each many-body cigcnstatc \a), with eigenenergy eVi^Ncx^a + 
Sa, is specified by a set {n^^^} of single-particle occupation numbers. How- 
ever, since such a single-particle description fails, e.g., for superconducting or 
ferromagnetic grains, we shall treat the \a) as general many-body eigenstates 
below. 

The tunneling term Htun has been written in two forms: Eq. (10) describes 
the tunneling of electrons between single-particle states in the leads and grain, 
with matrix elements T^^^ for tunneling events involving lead r; and in Eq. (11), 
Htnn is represented in the \a) eigenbasis, with corresponding matrix elements 
'^kaaa' giveu by Eq. (12). When doing perturbation theory in Htun, it turns 
out [72] that the T^;^ matrix elements always occur in the combination 

Ott 

^ E ^M'. - ekr) ^ Sw ri , (15) 

k 

which represents the total tunneling rate from lead r across barrier r into 
the single-particle state \la) on the grain. On the right-hand side we have 
neglected all off-diagonal terms, since for I ^ V the sum on k would involve 
matrix elements with randomly varying phases that average to zero, and also 
the a; dependence, which typically is of order uj je-p (since sets the energy 
scale for changes in the leads' density of states). 

2.3.3 Tunneling current and master equation 

The operator for the current carried by electrons tunneling between lead r 
and the grain, — —edNr/dt (defined to be positive if Nj. increases), is 

/. = ie[Nr, ^] = ie E na,aa' cLl«)(«1 + (h-c) . (16) 

ka,aa' 

Its steady-state expectation value = TrL j^ Q,(/3/r) in general requires knowl- 
edge of the system's full density matrix p. However, since the Hamiltonian is 
quadratic in the lead degrees of freedom, the latter can be integrated out, so 
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that the remaining degrees of freedom are described by the reduced density 
matrix P = TrL^j^(/3) for the grain. Its diagonal elements = {a\P\a) give 
the probability to find the grain in state \a). They satisfy, and may be found 
by numerically solving, a normalization condition and master equation of the 
form 

^Pa/ = 1 , = ^ i^aa'Pa' " ^a'aPa) for Cach a, (17) 

a' a'^a 

where T^aa' (for a 7^ a') is the total transition rate from initial state \a') to 
final state jo;). Furthermore, the current Ir can be written as 



^^=eEEp(sL:?-sw)^"'' (18) 

aa' p=l 

where S|^^f i^^aa^) is that part of the total rate S^q/ that involves the coherent 
transfer of a total of p electrons onto (from) the grain from (onto) lead r. Of 
course, charge conservation ensures that /l = /j^. 

Eqs. (17) and (18) have intuitively plausible forms, but the calculation of 
the rates T^aa' and S^^f is in general a highly non-trivial task. The standard 
strategy is to perform an expansion in powers of Htun, the most general and 
systematic formulation of which is the so-called real-time diagrammatic ap- 
proach developed by Schoeller, Schon and Konig [72-78]. Fortunately, for our 
present purpose of analysing tunnel-spectroscopic measurements on ultrasmall 
grains, we may restrict ourselves to the simplest possible situation: the tun- 
nel barriers are so large and the tunneling current so small that it suffices 
to consider only sequential tunneling ^ of electrons, described by lowest-order 
perturbation theory in i^tun- For this case one has [77,78] 



^aa' = E E Kl' , (19) 

r=L,Rp=± 

where E^^, = E^J^^I can be calculated using the golden rule. 



Kl' = X E fi^kr) \n:,^,JS{S^ - + eVu - Skr - eVr) (20) 

= /(^:a-^:a'-eK)r;+ , (21) 

^^aa' = Y ^ - fiSkr)] \TL,aa'?K^a ' ^a' ' eV^ + S^r + eK) (22) 

ka 



^ An exception will be encountered in Sec. 6.2, where cotunneling needs to be 
considered. 
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Fig. 4. Schematic illustration of the threshold conditions [see Eq. (26)] for bot- 
tleneck rates to become non-zero, (a) eVr > Sa'^^ ~ + ^atr+iGjv ^^"^ 
(b) -eVr > £a - £g'^^ - eVB for '^1~Gn+i^ for r = L and the case that the left 
barrier is larger than the right {T^ ^ T^). The long-dashed line indicates the equi- 
librium, V = chemical potential of the L and R leads relative to which eVr and 
~ =t eVo are drawn, using upward or downward arrows for positive and neg- 
ative energies, respectively. Note that excitation energies are always positive: (a) 
S^~^^ — £q + eVj) for particle-like excitations produced when an electron tunnels 
onto the grain; and (b) 8^ — £q'^^ — eVj) for hole-like excitations when it tunnels 
off. Filled circles depict the final electron configuration \a) after such a tunneling 
process, i.e. the electron number in the initial ground state before tunneling is drawn 
to be even (AT) in (a) and odd (iV + 1) in (b). 

= f{S^-Sa' + eVr)Tl-, , (23) 

r;+' = ErL , r;;, = Y.ri \{aMa')\' . (24) 

1(7 la 

Here f{E) = l/{e^/''^'^ +1) is the Fermi function and K = K - is the 
voltage drop (electrostatic potential difference) between lead r and the grain. 
For the standard case of a "symmetric" circuit with Vj^ = —V^ = V/2, which 
we shall henceforth adopt, it is given by 



L/R 



±y (Cr/l + Cg/2) - FgCg - Qo] /C , (25) 



which reflects the capacitative division of the bias voltage V{= V^l ~ = 
14 — ^r) across the two barriers. Eqs. (21) and (23) show that when V is 
swept, at a fixed gate voltage Vg and at temperatures much lower than the 
typical spacing d between the eigenenergies £a, a rate S^^, will be switched 
from "off" (exponentially small) to "on" (of order F^^,) each time eK passes 
through a threshold at which one of the inequalities 

± eVr >Sa- Sa' (26) 
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becomes true. Intuitively speaking, this occurs each time the energy gained 
by an electron when leaving (entering) lead r, namely ±eVr, becomes greater 
than the energy needed to enter (leave) the grain while inducing the transition 
I a') — ^ I a), namely Sa ~^a' ±eVD- Figs. 4 (a) and (b) illustrate this condition 
for the rates and T,^~ , respectively, and the case that the initial state is 
the ground state, jo;') = \G). 

The standard expressions for a mesoscopic (i.e. large) SET with a normal 
island, described by ^normal {'^•i): can be recovered as follows from 

the above formulas: firstly, one uses a factorized form for the island's density 
matrix, with — PmPun (and Y^a^f — !)■ Secondly, one makes the 
replacement 



where SEp^^ [Eq. (7)] is the change in electrostatic energy for the transition 
A^'cx Nex^l. For the second line we took the continuum limit for the energies 
£;o-D) and assumed that = F*" is independent of / and a. 

2.3.4 Eoctracting the fixed-N excitation spectrum 

The very large charging energies (£"0 ^ 5meV, ^ T, eV) of ultrasmall grains 
ensure that for a fixed value of gate voltage Vg, at most two kinds of charge 
states are involved in electron transport through the grain, i.e. Pa — unless 
Na equals or + 1 (say). To get a non-zero current, one needs non-zero 
rates for tunneling both onto and off the grain, through different barriers, say 
E'""'' , and S'" ~ , . To be specific, these will be of type 

EL+, E^- for ey > and E^- , E^+ for < [with e < 0], (29) 

since the direction of electron tunneling is always from the higher toward the 
lower chemical potential. One can distinguish two situations, depending on 
whether (a) only one or (b) both barriers are bottlenecks that limit transport. 

(a) One bottleneck. — If all the tunneling-on rates E''+ are generally, say, much 
smaller than any of the tunneling-off rates E^ ~ , then barrier r is the bottleneck 
for transport. The probability to find iV or -|- 1 electrons on the grain 
will then be close to unity or zero, respectively. Moreover, for T ^ d (and 
neglecting nonequilibrium effects, to be discussed in Sec. 6), the iV-electron 
ground state [Gat) will be the overwhelmingly most probable, i.e. P©^ ~ 1 and 



f{±ei,B + SE^ot T eK)/(T£;.D)F 



la 



= 2rW(£F) 



-(5E^,,TeVr) 



(27) 



(28) 
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Pq ~ for all other a. In this case the current through barrier r [Eq. (18)] 

is simply ~ eX^ajv+i (^o^+iGjv)- ■'■^ ^ swept at fixed V^, the current 
will thus show a steplike increase and the differential conductance a peak each 
time a rate '^at^^^Ct^ switched on, i.e. each time eVr increases past one of the 
threshold values E^^^ — Sq of (26), see Fig. 4(a). Denote this threshold value 
by and the corresponding threshold bias voltage [related to it via (25)] 
by Va- The voltage differences Vaa = Va — Va between any two conductance 
peaks can directly be translated into one of the excitation energies of the 
fixed-{N + 1) excitation spectrum, namely (for bottleneck barrier r =L/R) 

SN+l _ ^ ^ ^ya^ _ ya^^ ^ ±eV^^{C^/^ + C,/2)/C . (30) 

The case in which all rates S*"" are much smaller than any rate [Fig. 4(b)] 
is completely analogous to the one just discussed. Then the current is ^ 
—eJ2aN (^a^G]v+i)' ^'^'^ fixed- excitation energies are given by (for 
r =L/R) 

- 8^ = -e (Vl% - ^iJr) = TeV^a{Cn/L + C,/2)/C . (31) 

Eqs. (30) and (31) form the basis of RBT's tunnel-spectroscopic measurements 
of discrete spectra: a conductance curve showing a series of peaks as function of 
voltage can be replotted as function of energy by rescaling the horizontal axis 
using the voltage-to-energy conversion factors^ e(Cp^/L + Cg/2)/C in Eqs. (30) 
or (31). The rescaled curves, which we shall henceforth generically call "fixed- 
N excitation spectra", where is the number of electrons in the final state 
after the bottleneck tunneling process, then allow the corresponding excitation 
energies to be read off directly from the peak spacings. 

(b) Two bottlenecks. — The situation is more involved if the tunneling-on and 
tunneling-off rates E^"*" and are comparable in magnitude so that both 
barriers are bottlenecks: Then some of the Paj^ and PaN+i ^ill be compa- 
rable too, implying the same for the total probabilities to find or + 1 
electrons on the grain. The conductance will then show peaks corresponding 
to the switching-on of both and S'' ~ rates, at both eVr and eVr' values 
corresponding to both S^^^ — £a' £^ — 8^,^^ energy differences. It would 
thus be rather difficult to extract a purely fixed- A^ or fixed- (A^ + 1) excitation 
spectrum from the resulting mixed set of peak spacings (though in principle 
not impossible: if a magnetic field is used to switch the leads from being su- 
perconducting to normal, peaks due to or Y7 ~ shift by different amounts, 
see Eq. (34) in Sec. 2.4.3; or, if a gate is present, peaks due to S^"*" or ~ 
move in opposite V^-directions if Vg is changed, see Eq. (35) in Sec. 2.4.4). 

^ If no gate is present (Cg = 0) or if Cg ^ Cl, (as in the device of Fig. 2), this 
factor reduces to eCj^/L/(CL + Cr). 
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2-4 Experimental details 



This section is devoted to experimental details such as device fabrication and 
characterization, and what effects superconducting leads or a gate voltage have 
on the tunneling spectra. In particular, we explain how one can tell whether 
all conductance peaks correspond to tunneling across the same grain, whether 
a given peak has barrier L or R as bottleneck, and how the corresponding 
capacitance ratios can be determined very accurately [this knowledge is needed 
in Eqs. (30) or (31)]. 



2.4-1 Device fabrication 

The devices of RBT shown schematically in Fig. 1(a) were fabricated as follows 
[10]: First electron-beam lithography and reactive-ion etching were used to 
make a bowl-shaped hole in a suspended silicon nitride membrane, with an 
orifice between 5 and 10 nm in diameter [79]. The gate electrode was formed 
by evaporating 12 nm of Al onto the flat [bottom in Fig. 1(a)] side of the 
membrane. Plasma anodization and deposition of insulating SiO were then 
used to provide electrical isolation for the gate. Next an Al electrode which 
fills the bowl-shaped side [top in Fig. 1(a)] of the nitride membrane was formed 
by evaporation of 100 nm of Al, followed by oxidation in 50 mTorr O2 for 45 s 
to form a tunnel barrier near the lower opening of the bowl-shaped hole. A 
layer of nm-scale grains was created by depositing 2.5 nm of Al onto the 
lower side of the device; due to surface tension the metal beaded up into 
separate grains less than 10 nm in diameter [80,81]. In approximately 25% 
of the samples [determined as those exhibiting the typical SET "Coulomb- 
staircase" structure exemplified by Fig. 2], a single grain formed under the 
nm-scale tunnel junction to contact the top Al electrode. Finally, after a second 
oxidation step to form a tunnel junction on the exposed surface of the grain, 
a lower electrode is formed by evaporating 100 nm of Al to cover the grain. 

The resulting device is an SET whose central island is a nm-scalc grain fully 
coated by an insulating oxide layer, sitting on the lower lead electrode, encir- 
cled by a gate electrode, with the bowl part of the upper lead electrode at a 
fixed distance above it hke "an STM tip cast in concrete" . 

Recently, Davidovic and Tinkham succeeded in contacting leads to an ultra- 
small grain with a radius as small as r ^ 1 nm, using a somewhat different 
fabrication technique, described in detail in [82,15]. 
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Fig. 5. Current-voltage characteristics of an Au grain with radius r ~ 4.5 nm, from 
[15]. (a) I-V curve of a typical device at 4.2K. (b). I-V curve of the device at 30 mK. 
(c) First current step (a) of (b), at three different refridgerator temperatures. The 
curves between the points are fits to the Fermi distribution. After correcting for 
the capacitive division of voltage, the full-width-half-maximum (FWHM) of the 
peak, shown in inset, is linear with the refridgerator temperature, with a slope of 
FWHM/ZcbT = 3.8, which is close to the expected slope of 3.5. (d) Lineshape of the 
level when the leads are superconducting (H = 0) compared to the lineshape when 
the leads are normal (H = 1 T). 



2.4-2 Device characterization 



Only those devices are studied in detail whose I-V characteristics display the 
large-scale Coulomb-staircase [Fig. 2] expected for a SET, since its presence 
indicates that the tunneling current flows through a single grain only (as op- 
posed to several). The basic parameters of the SET (i?L, R^, Cl, C^, Cg and 
Ec) can be determined by fitting this Coulomb-staircase to the standard ex- 
pressions for mesoscopic SETs (with Cg — for devices without a gate), as 
described in detail by Hanna and Tinkham [62] . To estimate the grain radius, 
RBT crudely assumed [5-7,10] the grain to be a hemisphere with radius r and 
Vol = 27rr^/3 (actually their grains were more pancake- shaped) and deter- 
mined r by equating the (larger of the) grain-to- lead capacitances to nr^Cn] 
here Cq is the capacitance per unit area of the tunnel junctions. In Ref. [5], 
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Fig. 6. Current- Voltage characteristics for an Al grain connected to superconducting 
and normal leads [5], without a gate electrode (i.e. Cg = 0). The device parameters 
were Cl = 4.9 ± 0.5 aF, Cr = 8 ± 1 aF, Ri^/R^ = 8 ± 2, i?L + = 9 ± 1 MQ, 
Qo/e = 0.20 ± 0.03, and Eq ~ 6 meV. (a) / vs. V curves; the A^-lead curve shows 
steps, the 5-lead curve (displaced 10 pA in I) shows spikes, reflecting the BCS 
quasiparticle density of states in the leads [cf. Eq. (33)]. (b) and (c) d// dV vs. V 
for positive and negative V , with the 5-lead data shifted in V , as labeled, so as 
to align the maxima of dl/ dV with the A/^-lead data. For ease of comparison, the 
amplitude of the 5-lead data is reduced by a factor of 2 and offset on the d// dV 
axis in (b) and (c). 

RBT estimated that Cn — 0.075aF/nm^ for the AI2O3 layers insulating their 
grains from the leads and gates. More detailed subsequent studies [63] gave a 
slightly revised value of Cn — 0.050aF/nm^. Finally, the grain's single-particle 
mean level spacing d can be estimated using Eq. (2). 



2.4-3 Superconducting leads 

In many of RBT's samples, the leads are made from Al. They thus become 
superconducting below = 1.21K, but can be driven normal by turning on 
a magnetic field (the critical field is very small, about if ~ 20 mT, since 
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RBT's leads are in effect disordered tfiin films). Figs. 6 [and also 5(d)] show 
how the current and conductance curves change when the leads switch from 
normal to superconducting upon removal of this magnetic field: each current 
step is shifted to a higher \V\, relative to the normal- lead data, and takes the 
form of a spike, with, correspondingly, a region of negative dl / dV. This can 
readily be understood by redoing, for superconducting leads, the golden-rule 
calculations of Sec. 2.3.3 (see App. A): at if = and temperatures T ^ 
(the BCS gap in lead r), Eqs. (21) and (23) for the rates S'^ are modified to 

Ka' = fi^a - Sa' T eK) K,^{Sa - 8^' T eVr) T^^, , (32) 
E\/^E'' - A2 for |£;|>A,, 



7v;,qp(^) = < 



for < A^ 



(33) 



The novel feature in Eq. (32) is the appearance of N'r,qp{E), the BCS quasi- 
particle density of states in lead r, which reflects the fact that a quasiparticle 
is created each time an electron enters or leaves a lead. A detailed analysis 
[5] showed that the shape of the spikes in the superconducting-lead data in 
Fig. 6(a) [and also 5(d)] quantitatively agrees rather well with the prediction 
of Eq. (32), up to a slight amount of broadening near the threshold (cf. Fig. 3 
of [5]). A possible reason for this broadening was suggested by Levy-Yeyati et 
al. [83]: as a discrete state is swept past the threshold, its life-time broadening 
becomes anomalously large, since this broadening is proportional to the leads' 
density of states and hence reflects the BCS peak in the latter. 

Despite this slight amount of anomalous broadening, the conductance peaks in 
the dl / dV curve still are much sharper for superconducting than for normal 
leads, as is evident in Fig. 6(b). This is a very useful feature that can be 
exploited to enhance spectroscopic resolution. 

The shift in V"-thresholds for superconducting versus normal leads can be 
used to determine which barrier (r = L or R) acts as bottleneck for a given 
conductance peak: For superconducting leads, the threshold value for eVr at 
which a rate S^^ is switched on, namely ±eV^ > [E^ ~^a') + A^. [by Eq. (32)], 
is shifted by 5{eVr) = iA^ relative to the threshold of Eq. (26) for normal 
leads, reflecting the added energy cost for removing or adding an electron from 
or to a superconductor. By (25), this corresponds to a shift in bias voltage of 



Whereas the sign of SV always equals that of V [as follows by combining 
Eqs. (34) and (29), with e < 0], its magnitude evidently depends on r, allowing 
the bottleneck barrier r to be identified. Incidentally, once r is known, Eq. (29) 
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can be used to determine whether bottleneck transitions for a given sign of 
bias voltage involve tunneling onto or off the grain. 

Figs. 6(b) and (c) illustrate these ideas: they show two kinds of shifts, whose 
magnitudes 0.28 and 0.52 mV agree well with the predictions of Eq. (34) for L 
and R, using the sample parameters listed in the figure caption. In Figs. 6(b), 
all peaks are shifted by the same amount of 0.28 mV, thus L is the bottle- 
neck barrier, with ^a~G,^_^_^ bottleneck rates [by Eq. (29), since eV < 
for Figs. 6(b)]. Thus, this is an example of the "one-bottleneck" situation 
described in Sec. 2.3.4(a), from which fixed- A?" excitation spectra can be ex- 
tracted. In Figs. 6(c), both types of shifts occur, implying that this is the "two- 
bottleneck" situation described in Sec. 2.3.4(b): the peaks between —9 and 
— 12 mV are shifted by —0.52 mV and hence have bottleneck rates Tj^~r , 
whereas the peaks between —12 and —14 mV with shifts of —0.28 mV have 
bottleneck rates ^at+iGN- Fiii^lty) the fact that only two values of SV shifts 
were observed confirms that all current steps are due to tunneling through the 
same grain. 



2.4-4 Effects of gate voltage on tunneling spectrum 

In devices with a gate electrode, Vg can be used to controUably change the av- 
erage number of electrons on the grain, and to maximize or minimize nonequi- 
librium effects, as desired. This is illustrated in Fig. 7(a), which shows a series 
of differential conductance curves for a range of different gate voltages, for the 
same device as in Fig. 2 (whose caption gives the device parameters) . The bot- 
tleneck processes in its various regimes are depicted schematically in Fig. 7(b), 
which one can construct by the following arguments: First note in Fig. 7(a) 
that as Vg is increased, the extent of the Coulomb blockade region at low \V\, 
in which there are no conductance peaks, decreases, goes to 0, and then in- 
creases. At this zero crossing the SET is at a Coulomb blockade degeneracy 
point, at which, say, the and (A^-|- l)-clectron ground states are degenerate. 
When passing this point while increasing V^, which lowers el/o (since e < 0) 
and hence [by Eq. (6)] favors larger A^ex, the grain's average electron number 
changes from A^ to A^-|-l. Moreover, nonequilibrium effects (discussed in detail 
in Sec. 6) are weak close the the zero crossing but become stronger as Vg is 
tuned away from it, since their strength is governed, roughly speaking, by the 
size of the threshold voltage at which the first peak occurs. If one wants 
to make an "ideal" measurement of the "isolated" grain's properties in which 
the grain is disturbed as little as possible by the measuring process, one has 
to minimize nonequilibrium effects by tuning Vg close to a zero crossing. 

Next, note from Eq. (29) that the left half of Fig. 7(a) (where eV > 0) is 
governed by rates of the type and Ef;"*" , , and the right half by 

, and Ef"*" , , as illustrated in the left and right halves of Fig. 7(b), 



26 




-3-2-10123 

(a) V (mV) 




Fig. 7. (a) dl/dV vs. bias voltage V for the same sample [10] as in Fig. 2, plotted 
with different vertical offsets for a set of V^-values ranging from 75 mV (bottom) 
to 205 mV (top). Ah data are taken at T = 50 mK and H = 0.05 T (to drive 
the Al leads normal), (b) Schematic depiction (using the conventions of Fig. 4) of 
the processes acting as bottlenecks in the four quadrants of (a), for being odd, 
with rates: (i) '^a^ Gn+i (^PP^'^ quadrant), (ii) '^^n+i Gm (uPP^'^ lower left 
quadrants), (iii) ^ (upper and lower right quadrants), (iv) ^ (lower 

right quadrant). Solid (dashed) arrows depict bottleneck tunneling transitions into 
the lowest- (highest) energy final states accessible for the chosen value of V, and 
filled circles represent the electron configuration of the lowest-energy final state. 
Increasing Vg decreases eVo (recall: e < 0) and shifts the drawn grain levels down; 
the degeneracy point between the N- and {N + l)-electron ground states, at which 
the Coulomb blockade regime has shrunk to zero, is reached when the topmost 
non-empty level of the final ground state coincides with the equilibrium, V = Q 
chemical potential of the L and R leads, indicated by the horizontal long-dashed 
line. 
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respectively. To find out which of these rates (r = L or R) are bottlenecks to 
transport for a given V^, consider the change SV in the position of a conduc- 
tance peak upon a small change 6Vg in gate voltage; these changes are related 
by the condition that the corresponding threshold voltage Vr in (26) remain 
unchanged [by Eq. (25)]: 

SV _ Cg sv _ -c. 



A peak that has barrier L or R as bottleneck will thus move with slope 
5V/5V^ > or < 0, respectively. In Fig. 7(a), all peaks indeed do shift uni- 
formly with V^, which is another proof that all are due to tunneling through 
the same grain. In the upper left and lower right quadrants of Fig. 7(a), slopes 
of both signs occur (with magnitudes 0.009 and —0.025, consistent with the 
capacitances listed in the caption of Fig. 2); each of these quadrants thus cor- 
responds to a "two-bottleneck" situation [see Sec. 2.3.4(b)], in which both of 
the rates shown in the corresponding halves of Fig. 7(b) act as bottlenecks. 
In contrast, in the upper right and lower left quadrants all peaks move with 
slope 5V/5V^ > 0, implying that only the left barrier acts as bottleneck (con- 
sistent with the fact that Ri^/R^ ~ 18 1); each of these quadrants thus 
corresponds to a "one-bottleneck" situation [see Sec. 2.3.4(a)], from which ex- 
citation spectra can be extracted. The upper right quadrant has ^, as 
bottleneck rate and hence yields the fixed- excitation spectrum [cf. Eq. (31)], 
while the lower left quadrant has S^^^^ ^, as bottleneck rate and hence yields 
the fixed-(A^ -|- 1) excitation spectrum [cf. Eq. (30)]. 

Excitation spectra such as the latter two are the central measurement results of 
SET tunneling spectroscopy. One can proceed to extract from them all sorts of 
interesting information about the grain and the correlations which its electrons 
experience. For example, a very remarkable property of the fixed-(A^ -|- 1) 
excitation spectrum in the lower left quadrant of Fig. 7(a) is the sizeable 
distance between the first and second peaks (counted from \V\ — 0). As will 
be discussed in more detail in Sec. 4.1, this "spectroscopic gap" is evidence for 
superconducting pairing correlations on the grain: roughly speaking, it reflects 
the gap of 2Ag].ain (the pairbreaking energy cost) in the excitation spectrum of 
a superconducting grain with an even number of electrons. This also implies 
that A^ is odd. The behavior of these excitation spectra in a magnetic field, 
discussed in Sees. 4.1, 4.5 and 4.6, confirms these conclusions. 



2.4-5 Capacitance ratios 

The accuracy of tunnel-spcctroscopic measurements depends on the accuracy 
with which one knows the capacitance ratios {Cn/-^ + Cg/2)/C [or C'j^/l/(C'l + 
Cr) if no gate is present]. These occur, e.g., in the voltages Vi^/b. [Eq. (25)] that 
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enter in the threshold condition (26), or as voltage-to-energy conversion factors 
in Eqs. (30) and (31). These capacitance ratios can of course be determined 
from the large-scale Coulomb-blockade I-V curves of the grain. However, two 
additional methods are available, which are more accurate (typically to within 
1% [18]), since they exploit the properties of the discrete tunneling spectrum 
itself: 

(i) If the leads are made from a superconducting material, then the ratios 
(C'r/l + Cg/2)/C can be read off directly from the SV shifts [Eq. (34)] in peak 
positions that occur when superconductivity in the leads is switched off by 
switching on a magnetic field. 

(ii) One can compare the positions of two tunneling peaks occurring for 
opposite signs of the bias voltage but involving the same transition, e.g. 
|Gjv) — > \aN+i) (the case IGat+i) — \oin) is analogous). If eV > 0, this transi- 
tion will [by Eq. (29)] occur by tunneling across barrier L with rate S^^, and 
if eV < across barrier R with rate Sf^^. Denote the bias voltages at which 
the two conductance peaks occur by and V"'^, respectively, and use 
Eq. (25) to define from these two corresponding V"^-voltages, denoted by V{^'^ 
and By Eq. (26), these in fact are equal, V{^'> = K''^{^ - <^G )■ 
Using Eq. (25) to rewrite this equality in terms of and V"''^, one imme- 
diately finds that |\/"'>/\/"'<| = (Cl + Cg/2)/(CR + Cg/2). Adding unity to 
this or to its inverse gives the desired capacitance ratios. 



2.4-6 Ground state energy differences are currently not measurable 

According to Eq. (26), the threshold values of e\/" and eVa at which the 
first conductance peak occurs as \V\ is increased from should allow one, in 
principle, to also determine the ground state energy difference -Ec^vii ~ ^Gn 
between grains with adjacent electron numbers. In practice, however, it has so 
far not been possible to do this accurately: by Eq. (25), eV,°' depends through 
e{VgCg + Qo)/C on the gate voltage and the random offset charge Qo, and in 
the devices studied so far, cQq/C could not be determined with sufficiently 
high accuracy, i.e. with an uncertainty smaller than the scale of the grain's 
mean level spacing d, typically ~ 0.1 meV, for the following reason: 

Qo can of course be determined reasonably accurately by studying the large- 
scale Coulomb oscillations of the I-V curve that occur as functions of V"g at 
fixed y, a procedure that is well-estabhshed for mesoscopic SETs, for which 
indeed it has been possible to measure the ground state energy difference 
between a superconducting island with an even or odd number of electrons 
[62,84-88]. However, a complication arises for the nanoscopic grains of present 
interest, due to the smallness of their gate capacitances (typically ~ 0.1 aF): 
to sweep through one Coulomb oscillation, the gate voltage Vg must be 
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swept through a range so large (namely e/Cg IV) that during the sweep, 
RBT routinely observed small "rigid" shifts of the entire tunneling spectrum 
at random values of Vg. (Similar spontaneous shifts are discussed, e.g., in [88].) 
Likewise, shifts are also observed when eV is swept over a range of order Eq. 
These shifts presumably are due to single-electron changes in the charges con- 
tained in other metal grainsQ in the neighborhood of the grain of interest. 
These changes produce sudden shifts in the electrostatic potential energy of 
the grain, i.e. shifts in cQq/C by a few %, which are comparable in magnitude 
to the mean level spacing d and hence spoil any attempts to determine cQq/C 
with an uncertainty smaller than d. — The and Qq dependencies are sub- 
tracted out, however, when one considers the distances Vaa between peaks, as 
in Eqs. (30) and (31). Moreover, to measure only the first few discrete states, 
the range over which V needs to be varied typically is small enough that (with 
some luck) no shifts occur. 



' This interpretation is bolstered by the fact that the shifts usually occur at ap- 
proximately periodic values of V^. 
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3 Normal grains in an applied magnetic field 



In this section we describe how the abihty to resolve individual eigenenergies 
allowed RBT to determine the parity (even or odd) of the number of elec- 
trons, say N, in the V — ground state of a normal-state metal grain [5] 
[Sec. 3.1]. We also explain why orbital magnetism is negligible in ultrasmall 
grains [Sec. 3.2]. 



3.1 Breaking of Kramers-degeneracy by applied magnetic field 

In the absence of an applied magnetic field {H = 0), a normal-state grain 
will have time-reversal symmetry. For an even-N grain, the many-electron 
wave function for the ground state will be a spin singlet, in order that the 
orbital energy be minimized. In contrast, the ground state of an odd-A?" grain 
for H = necessarily is two-fold degenerate, by Kramers' theorem, forming 
a Kramers doublet. When H is turned on, this doublet is Zeeman-split by 
±^IJ,BgH. Therefore, for an even-N grain at small H, the lowest-lying tun- 
neling excitations correspond to transitions from the even-grain ground state 
singlet to the odd-grain ground state doublet, i.e. to two states split by H, 
so that the lowest-K conductance peak should exhibit Zeeman splitting in an 
applied field. On the other hand, for an odd-N grain with T <^ fi^g H/kB, the 
odd-grain ground state will be the lower-energy state of the Kramers doublet; 
the lowest-lying tunneling excitation will thus consist only of a single transi- 
tion from this odd-grain ground state to the even-grain ground state singlet, 
so that the lowest-T^ conductance peak should not split into two as a function 
of H. 

RBT observed both kinds of behavior, first in several different ungated Al 
grains [5,7], in each of which N has a fixed, random value, and subsequently 
also in a given gated Al grain, in which could be tuned via Vg [10,9]. 
Bulk Al superconducts, but the present grains were so small (r < 3 nm) that 
superconducting pairing correlations on the grain were negligible, since the 
mean level spacings (~ 0.5 meV) were substantially larger than the bulk gap 
A 0.18 meV for Al). Figs. 8(a) and (b) show the lowest-K conductance 
peaks for an even-A^ and an odd-iV grain, respectively. The lack of splitting of 
both of the prominent peaks in Fig. 8(b) indicates that the first two even-grain 
excited states for this grain are both spin singlets. The small peak in Fig. 8(b) 
(visible below the second large peak) , which moves to lower V with increasing 
H, is attributable to the noncquilibrium occupation of the higher-energy level 
of the odd- grain initial-state Kramers doublet. 

By measuring the difference in V between the Zeeman-split peaks and con- 
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Fig. 8. dl/dV vs. V at 50 mK and H = 0.03, 1, 2, 3 T from bottom to top, for 
two different samples [5]. The first transition above the Coulomb-blockade threshold 
exhibits Zeeman splitting in (a) but not in (b), hence the electron number N of the 
V = ground state is even for grain (a) and odd for grain (b). 

verting them to energy using Eqs. (30) or (31), RBT determined the (jf-factors 
for the states in Fig. 8(a) and (b) to be ^ = 1.87 ± 0.04 and ^ = 1.96 ± 0.05, 
respectively. In the former case the deviation from the value expected for pure 
Al, namely g^^^^ = 2, is significant, and can be attributed to spin-orbit scat- 
tering from the surface or impurities [20,61], a subject that will be discussed 
in detail in Sec. 7. 



3.2 Why orbital diamagnetism is negligible in ultrasmall grains 



Of course, a magnetic field in principle also couples to the orbital motion of 
the conduction electrons. Orbital effects in spherical and cylindrical supercon- 
ductors whose dimensions are smaller than the penetration depth were first 
considered by Larkin [89]. However, in grains as small as those of RBT, or- 
bital diagmagnetic effects are negligible [14], as can be seen from the following 
argument [90]: Orbital diamagnetism becomes important when the splitting 
which it produces in the eigenenergies of two time-reversed states (analogous 
to a Zeeman splitting) becomes comparable to the mean level spacing. This 
can be shown to occur [see Sec. 7.4.1] when the flux enclosed by an electron 
traversing a "closed" trajectory corresponding to a discrete quantum state 
exceeds one flux quantum $o {— hc/2e), i.e. when the magnetic fleld exceeds 
the scale -f/orb = ^o/^typj where A^yp is the typical (directed) area covered 
by the electron's trajectory during the "period" of its motion, h/d. Since 
the number of bounces off the grain's boundaries during this time is roughly 
5'dim = -EThouiess/c?, where ^/ii^Thouiess IS the time to cross the grain once, the 
directed area is r^^/g^, where the square root accounts for the fact 

that the direction of motion after each bounce is random [55]. It follows that 

H b ~ ^° = 2067 T.nm^ 
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Thus, /forb grows with decreasing grain size. Using d from Eq. (2) and the 
balhstic estimate of Eqs. (4) and (5) for Qdim, with a = 3, we find that hemi- 
spherical Al grains with radii of (say) r ~ 3 or 5 nm have i/oib ~ 19 or 7 T, 
respectively. If larger values are used for a, as would be appropriate for more 
pancake-shaped grains [cf. footnotes 3 and 4], iforb would be even larger. 

We may thus conclude that orbital diagmagnetic effects only begin to play a 
role for largish grains (> 5 nm), and then only for the highest fields (of 7 T) 
studied by RBT. 
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4 Superconductivity: experiment and phenomenological theory 

Among RBT's most striking experimental results are ttiose for Al grains: tliey 
found a significant spectroscopic gap indicative of superconducting pairing 
correlations in the fixed- excitation spectra of largish Al grains (r > 5 nm) 
if A^ is even, but not if it is odd (see Figs. 9 and 10 below); however, in their 
smallest grains (r < 3 nm) no such parity-dependent gap could be discerned. 

These results are of rather general interest and significance, since they invite 
reconsideration of an old but fundamental question: what is the lower size 
limit for the existence of superconductivity in small grains? Anderson [45] 
addressed this question already in 1959: he argued that if the sample is so 
small that its electronic eigenspectrum becomes discrete, "superconductivity 
would no longer be possible" when its mean level spacing d becomes larger than 
the bulk gap, which we shall denote by A. Heuristically, this is obvious (see 
Fig. 12 below): A/d is the number of free-electron states that pair-correlate 
(those with energies within A of ep), i.e. the "number of Cooper pairs" in the 
system; when this becomes < 1, it clearly no longer makes sense to call the 
system "superconducting" . 

Although Anderson's answer is correct in general, it generates further ques- 
tions: What, precisely, does "superconductivity" mean in ultrasmall grains, for 
which many of the standard criteria such as zero resistivity, Meissner effect 
and Josephson effect, are not relevant^? How is one to modify the grand- 
canonical BCS theory to obtain a fixed- A^ theory appropriate for ultrasmall 
grains, whose charging energy suppresses number fluctuations? What hap- 
pens in the regime c? > A in which superconductivity has broken down? Is 
the breakdown parity dependent? How is it influenced by a magnetic field? 
Sees. 4 and 5 are concerned with providing detailed answers to these and 
related questions. 

Sec. 4 is devoted to the experiments themselves. We analyze and qualitatively 
explain them in the framework of a phenomenological theory by Braun et al. 
[22,23,91], which offers a simple intuitive picture for visualizing the pairing 
correlations and the changes these incur when the grain size is decreased. Fur- 
ther theoretical developments, inspired by RBT's experiments but not directly 
concerned with their interpretation, will be discussed in Sec. 5. 

Sec. 4 is organized as follows: subsection 

^ For an isolated nm-scale grain, (i) its resistivity is not defined, since electron 
motion is ballistic and the mean free path is boundary-limited; (ii) the grain radius 
is smaller than the penetration depth, so that no Meissner effect occurs; and (iii) 
the electron number is fixed, so that the order parameter cannot have a well-defined 
phase. 
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Fig. 9. Excitation spectra for the same sample [10] as Figs. 2 and 7(a), measured 
at T = 50 mK and H = 0.05 T (to drive the Al leads normal), for four different 
V^-values, corresponding to different values for the grain's average electron number 
(from top to bottom: N + 1, N, N, N — 1). The curves are artificially offset on the 
vertical axis and each is labeled by the associated bottleneck tunneling rate 11^='= 
depicted in Fig. 7(b), the bottleneck barrier being r = L in this case. Plotted is 
dl/ dV vs. energy, where the latter is given by |eF|(CR + Cg/2)/C = 0.73|ey| 
[see Eqs. (30) and (31)]; the voltage-to-energy conversion factor reflects the voltage 
drop across barrier L. The sizeable spectroscopic gap between the first two peaks 
in the middle two curves, and its absence in the top and bottom curves, reflects 
the pairbrcaking energy cost in the excitation spectrum of a superconducting grain 
with an even number of electrons, and implies that is odd. 



(4.1) presents the experimental results of RBT; 

(4.2) proposes a phenomenological model for an isolated ultrasmall grain; 

(4.3) discusses how pairing correlations can be visualized in a fixed- system and 
explains w^hcn and in what sense it can be called "superconducting" ; 

(4.4) presents a generalized variational BCS approach for calculating the eigen- 
energies of various variational eigenstates of general spin | s) ; 

(4.5) discusses how an increasing magnetic field induces a transition from a pair- 
correlated to a normal paramagnetic state; 

(4.6) presents theoretical tunneling spectra of the RBT type, which are in quali- 
tative agreement with RBT's measurements; 

(4.7) explains how RBT's experiments give direct evidence for the dominance of 
purely time-reversed states in the pairing interaction; 

(4.8) discusses various parity effects that are expected to occur in ultrasmall 
grains. 
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Fig. 10. Magnetic field dependence [10] of excitation spectra such as those of Figs. 7 
and 9 and taken for the same grain, at (a) f« 110 mV and (b) Vg ~ 180 mV [lower 
left and upper right quadrants of Figs. 7(a), respectively]. Each line represents a 
distinct conductance peak in the dl/ dV curve and traces how its energy changes 
with H. Upward-moving peaks are broader and less distinct than downward-moving 
ones (for reasons poorly understood) and can be followed only for a limited range 
of H before they are lost in the background. The distances between lines directly 
give the grain's (a) fixed-(A'^ + 1) and (b) fixed-A'^ excitation spectrum, where 
is odd; the corresponding rate-limiting processes are schematically illustrated in 
Fig. 7(b,ii) for (a) and in Fig. 7(b,iii) for (b). The vertical dashed lines indicate the 
first four level-crossing fields Hg^g' (defined in Eq. (71) and assigned by comparison 
with Fig. 17, see Sec. 4.6), namely ifo,i =4T, i?i/2,3/2 = 4.25T, ifi,2 = 5.25T and 
^3/2,5/2 = 6.5T with uncertainty ±0.13T (half the iJ-resolution of 0.25T). 



4-1 A gap in the excitation spectrum 



RBT's results for largish Al grains (r>5 nm) are exemplified by Figs. 9 and 

10: if the number of electrons on the grain in the final state after the bottleneck 
tunneling process is even (middle two curves of Fig. 9), the excitation spectra 
display a spectroscopic gap between the first two conductance peaks that is 
significantly larger than the mean spacing between subsequent peaks, whereas 
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no such gap occurs if the final-state electron number is odd [top and bottom 
curves of Fig. 9). In other words, even-A^ excitation spectra (number parity 
p = 0) are gapped, whereas odd-A^ excitation spectra {p = 1) are not. This 
is even more clearly apparent when comparing Figs. 10(a) and 10(b), which 
give the magnetic-field (H) dependence of an even-N and odd-N excitation 
spectrum, respectively. However, in their smallest grains (r < 3 nm) no such 
clear parity-dependent spectroscopic gap could be discerned. 

BRT interpreted these observations as evidence for superconducting pairing 
correlations in their larger grains, using notions from the BCS theory of su- 
perconductivity: in an even grain, all excited states involve at least one broken 
Cooper pair, i.e. two BCS quasiparticles, and hence lie at least 2A above the 
fully-paired BCS ground state; in contrast, in an odd grain all states have at 
least one unpaired electron, i.e. at least one quasiparticle, and hence no sig- 
nificant gap exists between ground- and excited states. Fig. 7(b) is a cartoon 
illustration of this interpretation, if one imagines two electrons drawn on the 
same energy level to represent a "Cooper pair" (making this cartoon precise 
will be one of the main goals of this chapter): in Figs. 7(b,ii) and 7(b,iv) the 
final electron number is even, and all final excited states (reached via dashed 
arrows) have one less "Cooper pair" than the final ground state (reached via 
the solid arrow); in contrast, in Figs. 7(b,i) and 7(b,iii) the final electron num- 
ber is odd, and the final ground and excited states have the same number of 
"Cooper pairs" . 

The approximately hnear if-dependence of the excitation spectra in Fig. 10 
was attributed by RBT to the Zeeman energy shifts of discrete levels in a 
magnetic field, as discussed in Sec. 3.1. Indeed, the method described there 
for determining the grain's number parity applies here too: the fact that the 
lowest state in Figs. 10(a) or (b) does not or does display Zeeman splitting, 
respectively, confirms that N is odd. The reduction of the spectroscopic gap 
in Figs. 10(a) therefore is purely due to Zeeman energy shifts and has nothing 
to do with the reduction of the BCS gap parameter due to pair-breaking that 
occurs in bulk samples in a magnetic field [92]. A detailed discussion of the 
spectra's magnetic field dependence will be given in Sees. 4.5 and 4.6. 

For completeness, it should be remarked that a spectral gap in ultrasmall 
superconducting grains was observed as long ago as 1968 by Giaever and 
Zeller [80,81], who studied tunneling through granular thin films containing 
electrically insulated Sn grains. They found gaps for grain sizes right down to 
the critical size estimated by Anderson (radii of 2.5 nm in this case), but were 
unable to prove that smaller particles are always "normal" . RBT's experiments 
are similar in spirit to this pioneering work, but their ability to focus on 
individual grains makes a much more detailed study possible. 



37 



4-2 A model for ultrasmall grains with pairing correlations 

In this section we construct a model for an isolated ultrasmall grain with 
pairing correlations, using phenomenological arguments valid for the regime 
(i < A. The model, which we shall call "discrete BCS model", allows a rather 
detailed qualitative understanding of the measurements of RBT [22,23] and 
hence is "phenomenologically successful" . For d ^ A it is unrealistically sim- 
ple, however, and should rather be viewed as a toy model for studying how 
pairing correlations change as the grain is made smaller and smaller. 

4-2.1 A simple reduced BCS interaction plus a Zeeman term 

Following the philosophy of the "orthodox theory" for Coulomb blockade phe- 
nomena mentioned in Sec. 2.3.1, we assume that the only effect of the Coulomb 
interaction is to contribute an amount -Ec-^ex [cf- Eq. (6)] to the eigenenergy of 
each eigenstate of the grain. Since the charging energy is huge (5 to 50 meV) in 
ultrasmall grains, this term strongly supresses electron number fluctuations, so 
that, to an excellent approximation, all energy eigenstates will also be number 
eigenstates.[^ Since EqN"^^ is constant within every fixed- Hilbert space, we 
shall henceforth ignore it, with the understanding that the model we are about 
to construct should always be solve in a fixed- Hilbert space (and that the 
use of grand-canonical approaches below, after having dropped EqN^^, simply 
serves as a first approximation to the desired canonical solution). 

The only symmetry expected to hold in realistic, irregularly-shaped ultrasmall 
grains at zero magnetic field is time-reversal symmetry. We therefore adopt a 
single-particle basis of pairs of time- reversed states |j±), enumerated by a dis- 
crete index j. Their discrete energies ej are assumed to already incorporate the 
effects of impurity scattering and the average of electron-electron interactions, 
etc. As simplest conceivable model that incorporates pairing interactions and 
a Zeeman coupling to a magnetic field, we adopt a Hamiltonian H = ifo + -f^red 
of the following reduced BCS form: 

Ho= J2 i^j - ^ - (^h)c]^Cj^ , i^red = -Xd^cl+cl_Cj_Cj+ . (37) 

Here —ah = a^nsgH is the Zeeman energy of a spin a electron in a magnetic 
field H, and we shall take h > below. Models of this kind had previously 
been studied by Richardson [46-54], Strongin et al. [93], Miihlschlegel et al. 
[94,95] and Kawataba [96,97]. The first application to RBT's grains for h = 
was by von Delft et al. [21] and for /i 7^ by Braun et al. [21-23]. 

^ An exception occurs at a degeneracy point, where -E'pot(-Acx) = Epoti^ex + 1) 
[cf. Eq. (6)]; the resulting complications will not be considered here. 
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Due to level repulsion the e/s will, to first approximation, be uniformly spaced. 
Unless otherwise specified, we shall for simplicity always (except in Sec. 5.4) 
take a completely uniform spectrum with level spacing d. Fluctuations in 
the level spacings have been studied with methods of random matrix theory 
[24,39], with qualitatively similar results (see Sec. 5.4). For a system with a 
total of electrons, where the electron number parity p = A^mod2 is equal 
to for even and 1 for odd A^, we use the label j = for the lowest-lying 
non-doubly-occupied level (with occupation number p) in the T = Fermi 
sea, which we shall denote by |Fjv). We choose the Fermi energy at = 
write[^ 

£. = jd + (1 - p)d/2 , (38) 

thereby taking the doubly-occupied and empty levels of [Fat) to lie symmetri- 
cally above and below (see e.g. Fig. 13 below). The parameter /i in Eq. (37) 
is, in g.c. theories, the chemical potential, whose value^° determines the av- 
erage particle number. For canonical theories, which make no reference to a 
chemical potential, fi is not needed and can be dropped (i.e. set equal to 0). 

The pairing interaction is of the reduced BCS form, in that it scatters a pair 
of electrons from one pair of time-reversed states into another. It is taken to 
include only states whose energy separation from the Fermi energy lies within 
the cutoff given by the Debye frequency: \ej\ < uj-d. Experimental support 
for the sufficiency of neglecting couplings between non-time-reversed pairs of 
states]^ i.e. of using only a reduced BCS-Hamiltonian, is given in Sec. 4.7. The 
pair-coupling constant in Eq. (37) is written as Xd, where A is a dimensionless 
parameter independent of the grain's volume, to make it explicit that both 
Hq and Hj-ed make extensive (oc Vol) contributions to the ground state energy 
(since the number of terms in each sum in Eq. (37) scales with A^, and 
d oc Ep/N). The "bulk gap" of the model, obtained by solving the standard 
BCS gap equation [Eq. (57)] at T = in the bulk limit, thus is 

A = cjD/sinh(l/A) . (39) 

To be precise, by "bulk limit" we shall always mean rf/A — >■ and N ^ oo 
while the product Nd is kept fixed, and use dJ2j ^ I dej. 

This convention differs slightly from that used in [21-23], namely ej = jd+so. The 
latter is a little less convenient, resulting, e.g., in a p-dependent chemical potential 
for the variational BCS ground states discussed below, /Up*-^^ = + {p — l)d/2, 
whereas (38) results simply in ^p*^^ = 0. 

A theoretical motivation for the simple form assumed for the pairing interaction, 
based on random matrix theory, may be found in Ref. [13] and is briefly mentioned 
at the beginning of Sec. 6.1.3. 
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An applied magnetic field will completely penetrate an ultrasmall grain, since 
its radius (typically r < 5nm) is much smaller than the penetration length 
of 50 nm for bulk Al. The Zeeman term in Eq. (37) models the fact that 
RBT's measured tunnel spectra of Fig. 10 evolve approximately linearly as a 
function of magnetic field, with g factors betweenf^ 1.95 and 2 (determined 
from the differences between measured slopes of up- and downward-moving 
lines). Deviations from g = 2 probably result from spin-orbit scattering, known 
to be small but nonzero in thin Al films [98,99], but neglected below (where 
g = 2 is used). Furthermore, orbital diamagnetism is also negligible, just as 
for thin films in a parallel magnetic field [98,99] but in marked contrast to 
bulk samples where it causes the Meissner effect: as explained in Sec. 3.2, 
a magnetic field only begins to significantly affect the orbital motion of the 
electrons once it exceeds the scale i^orb — "^oA^^/^dim [90,55]; but according to 
the estimates mentioned after Eq. (36), RBT's ballistic grains with r < 5 nm 
have Horh ^ 7 T, hence orbital effects should set in only near the highest fields 
of 7 T studied by RET. Indeed, some larger grains do show slight deviations 
from //-linearity [6] for large fields, which probably reflect the onset of such 
orbital effects [100]; however, these are much smaller than Zeeman effects in the 
grains of present interest, and will be neglected here. Thus, the model assumes 
that Pauli paramagnetism due to the Zeeman energy completely dominates 
orbital diamagnetism, similarly to the case of thin films in parallel magnetic 
fields [98,99]. 

Intuitively speaking, it is clear that the discrete BCS model introduced above 
contains all ingredients necessary to make contact with the spectra of Fig. 10: it 
is formulated in terms of discrete levels, it contains a pairing interaction which 
is known, from bulk BCS theory [101,92], to cause a gap in the excitation 
spectrum, and it contains a Zeeman term that will cause eigenenergies to 
linearly depend on an applied magnetic field. Indeed, we shall see in Sec. 4.6 
that it can be used to obtain a rather detailed qualitative understanding of 
the spectra of Fig. 10. 



4-2.2 Choice of numerical values for model parameters 

When doing numerical calculations for this model, some choices must be made 
for the numerical values of its parameters (though slight changes in their 
values will not change the results qualitatively). We shall follow the choices 
made by Braun et al. [23], since these led to reasonable agreement between 
experimental and theoretical excitation spectra. For the Debye frequency they 



Claims of smaller g factors made in Ref. [6] are wrong, the result of confusing 
different orbital states as Zeeman-split spin states. This was made clear in Ref. [10], 
where it was observed that upward-trending Zeeman states can have significantly 
smaller amplitude than downward-trending states, making them difficult to observe. 
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used the textbook value [102] for Al of ud = 34meV. Making an appropriate 
choice for the "bulk gap" A is less straightforward, since its experimental 
value for systems of reduced dimensionality often differs from that of a truly 
bulk system, presumably due to (poorly-understood) changes in the phonon 
spectrum and the effective electron-phonon coupling. For example, for thin 
Al films [93,103] it is known that Athinfiim — 0.38 meV, which is about twice 
as large as the gap of a truly bulk system, Abuik = 0.18 meV. (This increase 
in A is not universal, though; e.g., for Nb A is smaller in thin films than in 
the bulk.) Since ultrasmall grains are in many ways analogous to thin films 
in a parallel magnetic field [see Sec. 4.5], Braun et al. adopted the thin-film 
value for grains too, i.e. used A ~ 0.38meV. These choices imply that the 
dimensionless pair-coupling constant A = [sinh^^(c<jD/A)]~^ [cf.Eq. (39)] has 
the value A = 0.194. (In Sec. 4.6 we shall see, a posteriori, that the choices A = 
0.34 and A = 0.189 would have been slightly more appropriate.) Finally, for 
those numerical calculations that are explicitly cut-off dependent, Braun et al. 
smeared the cutoff of the BCS interaction over two single-electron levels; this 
smooths out small discontinuities that would otherwise occur in rf-dependent 
quantities each time the energy \sj\ of some large-|j| level moves beyond the 
cutoff cud when d is increased. 

Note that the above way of choosing A lumps into a single phenomenological 
constant all the poorly-understood effects of reduced dimensionality [93] on 
the phonons that mediate the attractive electron-electron interaction. Study- 
ing these effects in detail would be interesting in its own right, but would 
require systematic investigations with grains of well-controlled shapes and 
sizes. For the case of RBT's irregularly-shaped grains, using a phenomeno- 
logical coupling constant seems the best one can do. Note, though, that the 
precise value of A is not very important as long as all energies are measured 
in units of A (as we shall do for all numerical calculations), since most of 
the A-dependence is thereby normahzed away. Therefore, the shght difference 
between the A-values proposed above and those used in [35,91,39] (namely 
0.224) hardly matters. 



4.2.3 Some general properties of the eigenstates - the blocking effect 

The eigenstates of the discrete BCS model of Eq. (37) have some simple but 
general properties that are worth stating at the outset. 

Firstly, every eigenstate of H will also be an eigenstate of the number operator 
N = ct^c^.^, since [H, N] = 0. 

Secondly, since the interaction only involves levels within the cutoff energy 
cud of £f, the dynamics of those lying outside this range is trivial. We shall 
thus ignore them henceforth and focus only on the remaining set of interacting 
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levels, denoting this set by /. 



Thirdly, singly-occupied levels do not participate in the pair scattering de- 
scribed by H: "unpaired" electrons in such levels are not scattered to other 
levels, hence the labels of singly-occupied levels are good quantum numbers. 
Moreover, every unpaired electron Pauli-blocks the scattering of other pairs 
into its own singly-occupied level, i.e. it restricts the phase space available 
to pair scattering and thereby weakens the amount of pairing correlations, as 
we shall see in detail later. This was called the "blocking effect" by Soloviev 
[104], who discussed it extensively in the early 1960's in the context of nuclear 
physics. The eigenstates \a) and corresponding eigenenergies £a, of H thus 
have the following general forms: 

|q;) = |*„,S) = n4l*n>' (40) 

U n 
l*n)= E V'(il,---,Jn)n^iJVac), (41) 



This describes N = 2n + b electrons, b of which are unpaired and sit in a set 
B of singly-occupied, blocked levels, making a contribution Esih) to the total 
eigenenergy. The remaining n pairs of electrons, created by the pair operators 
&j = c]_,_c]_, are distributed among the remaining set U — I\B of unblocked 
levels, with wavefunction ■0(ji, . . . ,_;„) = ^j^i\B denotes a sum over all 
unblocked levels in I). The corresponding state is an eigenstate of the 
pair number operator and a Hamiltonian Hu involving only pair operators: 

E&>,|*n) = Hu\^n) = Sn\^n) , (43) 

j 

Hu ^jZi^isj - l^)5ii - Xd]blbj. (44) 

Each eigenstate |\E'„, B) may be visualized as a coherent superposition of eigen- 
states of Ho that all lie in the same fixed- Hilbert space, and in all of which 
each pair of unblocked (j G U), time-reversed levels |j±) is either doubly 
occupied or empty. This is illustrated in Figs. 11(a) and (b), which schemat- 
ically depict the exact ground states for even and odd N, respectively. The 
odd ground state has a single blocked level, at the Fermi energy, containing 
an unpaired electron. The latter somewhat weakens pairing correlations rela- 
tive to the even ground state and hence leads to parity effects, which will be 
extensively discussed in later sections. 
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Fig. 11. A cartoon depiction of the exact ground state for a reduced BCS Hamilton- 
ian, for N being even (a) or odd (b): they are coherent superpositions of eigenstates 
of Hq (whose respective amphtudes are not depicted) that aU have the same A''; 
the leftmost of these is (a) the even or (b) the odd Fermi ground state |F^), whose 
Fermi energy is indicated by the wavy hne. 

A useful measure for the amount of energy gained by \a) via its correlations 
is its "condensation energy" relative to the uncorrelated state \ck)o, 

= ^a- o{a\H\a)o, where |a)o = U clp)o , (45) 

ieB 



and \U)o is the "Fermi ground state" in U, for which the n pairs occupy the 
n lowest-lying levels in U. 

Note that Hu is /i-independent, since the total Zeeman energy of any pair 
of electrons is zero. Hence the full /i-dependence of the eigenenergies resides 
in the rather trivial contribution Ssih) of the blocked levels, which is a very 
important and useful simplification. 

Diagonalizing would be trivial if the 6's were true bosons. However, they 
are not, and in the subspace spanned by the set U of all non-singly-occupied 
levels instead satisfy the "hard-core boson" relations. 



by = 0, [b^, b],] = 5,y(l - 2b]b^), [b]b^, b],] = 5jj,b] , (46) 

which reflect the Pauli principle for the fcrmions from which the 6's are con- 
structed. In particular, b^ = imphes that only those terms in (41) are 
non-zero for which the indices ji, . . .jn are all distinct. 

The task of finding the eigenstates is thus highly non-trivial. Nevertheless, 
an exact solution does exist. Unbeknownst to most of the condensed- matter 
physics community, it was found and studied extensively by Richardson in the 
1960's and will be presented in Sec. 5.1. Throughout the present section 4, 
however, we shall use more well-known approaches based on the variational 
wavefunctions introduced by BCS [101], and that had been used to study the 
discrete BCS model before Richardson's solution was revived towards the end 
of 1998. 
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4-3 Canonical characterization of pairing correlations 

Since the discrete BCS model has the standard reduced BCS form, the most 
natural first step toward understanding its T = properties is to use BCS-like 

variational wavefunctions (or cquivalently Bogoljubov's mean-field approach), 
and indeed this will be done in Sees. 4.4 to 4.8. However, the limitations of 
such an approach should be realized at the outset: the spectra measured by 
RBT are excitation spectra for a grain with a fixed electron number N, and 
hence should be calculated for a grain with definite electron number N (i.e. 
completely isolated from the rest of the world, e.g. by infinitely thick oxide 
barriers). In contrast, the variational wavefunctions of BCS [Eq. (47) below] 
do not have the fixed- form [Eq. (41)] which any true cigcnstate should have, 
but instead are formulated in a grand- canonical (g.c.) framework (as is the 
Bogoljubov mean- field approach to which they are intimately related). 

When considering a truly isolated superconductor such as a perfectly insulated 
grain (another example would be a superconductor levitating in a magnetic 
field due to the Meissner effect), one therefore needs to address the following 
question, which will be the main theme of the present section: how is one to 
incorporate the fixed-N condition into BCS theory, and how important is it to 
do so? This issue is well understood and was discussed at length in the early 
days of BCS theory (Rickayzen's book gives a beautiful discussion [105]), in 
particular in its application to pairing correlations in nuclei [106, p. 439] (see 
also the general remarks in [107]). Nevertheless, for pedagogical reasons the 
arguments are well worth recapitulating in the present context. 

We shall first remind the reader that the use of a g.c. framework is only a 
matter of convenience, since the essence of the pairing correlations that lie 
at the heart of BCS theory is by no means inherently g.c. and can easily be 
formulated in canonical language [21]. We then show how standard BCS theory 
fits into this scheme, point out that the differences between results obtained 
using g.c. and canonical wavefunctions are negligible for ^ A, and conclude 
that for the purposes of gaining a phenomenological understanding of the 
experimental data, standard grand-canonical BCS theory should be sufficient. 
Nevertheless, the fundamental question of how to improve on this theory, in 
order to achieve a truly canonical description and to properly treat fluctuation 
effects, which become important for > A [108,109,28], is interesting and 
important in its own right and will be addressed at length in Sec. 5. 

For simplicity, throughout the present subsection 4.3 we shall consider only the 
even ground state in the thermodynamic limit (in which even-odd differences 
are negligible), so that U — I and blocking effects need not be worried about. 
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4-3.1 The grand- canonical BCS wavefunction 

Conventional BCS theory describes the pairing correlations induced by an 
attractive pairing interaction such as H^-ed of Eq. (37) within a g.c. ensemble, 
formulated on a Fock space of states in which the total particle number N is 
not fixed. This is illustrated by BCS's famous variational ground state Ansatz 

|BCS) = lliuj + e'^^ Vjh"^^ |Vac) , with u] + v'j = 1, (47) 

3 

where the variational parameters Uj and Vj are real and is a phase (which, 
it turns out, must be j'-independent, for reasons discussed below). |BCS) is 
not an eigenstatc of and its particle number is fixed only on the average 
by the condition (iV)Bcs = -A) which determines the g.c. chemical potential 
Likewise, the commonly used g.c. definition 

Age = A d E(^i)BCS = A d ^ UjVj e'^^ (48) 
j j 

for the superconducting pairing parameter only makes sense in a g.c. ensemble, 
since {hj) would trivially give zero when evaluated in a canonical ensemble, 
formulated on a strictly fixed- A^ Hilbert space of states. (We shall use the term 
"pairing parameter" instead of "order parameter" , since the latter carries the 
connotation of a phase transition, which would require the thermodynamic 
limit N —>■ oo, which is not applicable for ultrasmall grains). 

4-3.2 Canonically meaningful definition for the pairing parameter 

A theory of strictly fixed-A^ superconductivity must evidently entail modifi- 
cations of conventional BCS theory. However, these are only of technical, not 
of conceptual nature, since the essence of the pairing correlations discovered 
by BCS can easily be formulated in a canonically meaningful way, including a 
definition for the pairing parameter. We shall now attempt to explain, in in- 
tuitive, non-technical terms, how this may be done (our discussion is indebted 
to that of Rickayzen [105]). Readers with a preference for rigor may consult 
Sec. 5 for a corroboration, using Richardson's exact solution, of the arguments 
presented below. 

Let \G) be the exact even ground state of the system, depicted in Fig. 11(a). 
As explained in Sec. 4.2.3, it is a coherent superposition of eigenstates of Hq 
that all have the same A^ and in all of which each pair of time-reversed levels 
|j±) is either doubly occupied or empty. Due to this coherent superposition, 
|G) entails strong pairing correlations, whose essential properties may be un- 
derstood by investigating how they modify the correlators 
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^ {h\h^) , v] ^ C,, = {b%) , u'^ ^ {b^b]) , (49) 
relative to the form these take on for the Fermi ground state IFtv): 



(a,)F = 6.,{vf)F , (v^)^ = 9{-e,) , (uj)^ = 9{e,) . (50) 

Cij{= C*j) is the matrix element for the interaction to be able to scatter a 
pair of electrons from level j to i, and w| and are the probabilities to find 
level j doubly occupied or empty, respectively. The pairing correlations in 
|G) must be such that ^^red lowers the ground state energy below that of the 
uncorrelated Fermi sea iF^v) by an amount that is extensive {(x N (x d~^) 
in the thermodynamic limit. Clearly, this requires that (-ffred)G — (-^^red)F is 
negative and extensive, i.e. that 



XdY, [Cij - (C,j)f] ^ a ^ ^ 2Re{Cij) oc (and positive). (51) 

ij i j<i 

In the second expression we neglected the diagonal terms, since their number 
is so small (only oc A^) that XdJ^ji^] — is at best of order unity in the 

thermodynamic limit. For Eq. (51) to hold, |G) must have two properties: 

(i) the number of Ci/s that differ significantly from zero (i.e. are of order 
unity) should scale like iV^, i.e. one power of per index [105, p. 167]; 

(ii) most or all of the Cij for i < j should have the same phase, since a sum 
over random phases would average out to zero. 

Since a suitable pairing parameter should vanish in the thermodynamic limit 
unless both these conditions hold, the definition 



AL = (A dY - (4c,-+) (cl_c,-_)) (52) 

(or its square root) suggests itself, where the subscript emphasizes that (52) 
is meaningful in a canonical ensemble too, and we subtracted]^ the "normal- 

This subtraction was suggested to us by Moshe Schechter, who pointed out that 
then Eq. (52) has a natural generalization to position space: it is the spatial average, 
■^can = (A d)'^ J dri dr2 .7^(ri, 12), of the two-point function 

^{ri,r2) = (V'l(n)^l(n)V'-(''2)V'+('^2)) - (V'l(n)'0+(^2))('0l(n)V'-('^2)) 

(with ipa{r) = VoP^/^ 6''^'^ '^fco-)' '^hich evidently measures the amplitude for the 
propagation of pairs as opposed to uncorrelated electrons. Other definitions for a 
canonically meaningful pairing parameter have been suggested [21,35,23], such as 
Ac?^^- UjVj or Xd'^jKbjbj) — (ct_j_c^.^)(ct_c^._)]-^/^ , but these focus only on require- 
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state contribution to Cij." If (i) and (ii) hold, Acan will take on a finite value; 
its relation to a gap in the spectrum will become clear below. In the bulk limit, 
Acan can be shown [see Sec. 5.3] to reduce to the "bulk pairing parameter" A 
of Eq. (39). 



4-3.3 Redistribution of occupation probability across s-p 

Now, property (i) can be realized if all C^j in a finite (d-independent) range 
of £j's and e/s around the Fermi energy differ significantly from zero; the 
width of this range will evidently determine the magnitude of Acan (provided 
(ii) also holds), which conversely can be viewed as a measure of this width. 
But a nonzero Cij evidently requires both blbj\G) ^ 0, implying {yj)G 7^ 
and ('Uj)G 7^ 0, and also {G\blbj ^ 0, implying (^1)0 7^ and {uj)G 7^ 0. The 
product {ujVj)G must thus be different from zero [in contrast to {ujVj)^ = 0] 
for all Ej within a finite range around ep (cf- Fig. 12). This can be achieved by 
smearing out the sharp steps of the ^-functions of (vjOf and {uj)p, so that {vj)G 
[or ('Uj)g] is nonzero also for a finite range of Ej above [or below] Ep. In other 
words, for |G) some occupation probability must be redistributed (relative to 
IFat)) from below to above ep, as illustrated in Fig. 11. This redistribution, 
which was called pair-mixing in [21,23], frees up phase space for pair scattering 
and so achieves a gain in interaction energy (provided (ii) also holds) that more 
than compensates for the kinetic energy cost incurred thereby. 

Furthermore, note that properties (i) and (ii) imply, even without detailed cal- 
culations, that the spectrum will be gapped. Consider, for example, a "block- 
ing excitation" that disrupts pairing correlations by having \j+) definitely 
occupied and |j — ) definitely empty; since pair-scattering involving level j is 
blocked, the energy cost is 

{ej -I,)- [{ej - i^)2{b%)) - A d E (^1^.- + (53) 
= {e, -/,){! -2vf)+XdY: i^j + Cj,) , (54) 

in which the restriction on the sum refiects the blocking of scattering events 
involving level j. Since the first term of (54) is positive definite (particle- 
hole symmetry ensures that {\ — Vj) < if £j — < 0) and the second of 
order Acan, the excitation energy will be finite even for d — > 0, implying 
the existence of a gap of order Acan- Similarly, "phase-breaking excitations" 
that violate the fixed-phase condition (ii) are gapped too: for example, if 



ment (i) and fail to incorporate requirement (ii). A quantity very similar to Eq. (52) 
was recently proposed in Eq. (55) of Ref. [30], namely (Ad)^ Z^jJC'jj — {Cij)\=o]. 
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Fig. 12. A cartoon depiction of why "superconductivity breaks down" when the 
sample becomes sufficiently small. Vertical Unes are drawn at each single-particle 
energy Sj, spaced with a mean level spacing d corresponding to (a) a "large" grain 
(d < A); (b) a "small" grain {d ~ 0.25 A); (c) an "ultrasmall" grain ~ A). 
The lines' height represents the function UjVj = ^A^/(ej + A^) of standard bulk 
BCS theory, to illustrate the energy range (of width A around £f) within which 
pairing correlations are strongest. Loosely speaking, the number of single-particle 
levels in this regime, namely A/d, corresponds to "the number of Cooper pairs" 
of the system. Evidently, this number becomes less than one when d > A as in 
(c), so that it then no longer makes sense to call the system "superconducting" [cf. 
Sec. 4.4.2]. 



(C'ij)excited = -(Cii)ground for a given j and all j), the energy cost is 
-Xd Ei(^j)[(Cij+C*j)excited-(Cjj+Cji)ground], which is at least of order 2Acan. 

We see, therefore, that the essence of pairing correlations can readily be formu- 
lated in a canonical framework: (i) a redistribution of occupation probability 
across occurs, such that each level j in a finite range around Sp has a finite 
probability of both being doubly occupied or empty, and (ii) any two compo- 
nents of the ground state wavefunction that differ only by the exchange of a 
pair of electrons between two levels i and j have the same phase. 

Since pairing correlations with these properties are the microscopic property 
at the heart of all manifestations of "superconductivity" , it seems reasonable 
to call a sample "superconducting" as long as it exhibits pairing correlations 
with measurable consequences. And by this criterion the gap observed in the 
even grains of RBT certainly qualifies. 
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4-3.4 Gauge symmetry breaking 

Note that property (ii) will be preserved under the gauge transformation 
Cja- Q^^'^ Cja, i.e. Cij — > e'^^^'^'i~'^'j^ Cij, only if all (p'j are equal, say 0^- = (f)'. 
Property (ii), and likewise the pairing parameter Acan- tlicrcfore (a) are not 
gauge invariant "locally" in j-spacc, but (b) are gauge invariant globally. These 
are obvious consequences of the facts that (a) a correlated fixed-N state con- 
sists of a phase- coherent superposition of many different components, and 
hence cannot be invariant under arbitrary changes of the phases of individual 
components; and that (b) all of these components contain the same number 
of electrons and hence under a global gauge transformation all pick up the 
same phase factor e'^*^ . Obviously, global gauge symmetry can therefore never 
be broken in a canonical ensemble. In contrast, the breaking of global gauge 
symmetry by the g.c. pairing parameter Age of Eq. (48), which transforms as 
Age — > e^"^' Age, is an inevitable consequence or artefact of its g.c. definition 
[105, p. 142]. 



4-3.5 Making contact with standard BCS theory 

One of the breakthrough achievements of BCS was, of course, to propose a 
simple variational ground state which has precisely the properties (i) and (ii) 
described above: when evaluating the correlators of Eq. (49) using |BCS) of 
Eq. (47), one finds 

(%)bcs = (^j)lcs = > (Cy )bcs = UiViUjVj e-'<^i-'t>i) , (55) 

and also (A^3^^)bcs = |Agep. The definite-phase requirement (ii) can thus be 
implemented by choosing all the phases to be the same, say (pj = (j) for all 
J, thereby breaking local gauge invariance (usually one simply takes = 0); 
and requirement (i) is fulfilled automatically when minimizing the expectation 
value {H)bcs w.r.t. uj and Vj, since this does yield smeared-out step functions, 
namely [101,92] 

v] = I [1 - (sj - i,)/Ej] , Ej = ^(£,.-/.)2 + |AgeP . (56) 

Here we neglected terms that vanish for d — > 0, and Age is determined by the 
famous gap equation (for T — 0), 



\ej\<u>u ■> 

The BCS wavefunction instructively illustrates some of the general properties 
discussed above. Firstly, the product u^jV^, shown in Fig. 12, has a bell-shaped 
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form with a well-developed peak around ep of width ~ |Agc|, illustrating that 
pairing correlations are strongest within a region of width |Agc| around the 
Fermi surface. Secondly, the energy of a blocking excitation [Eq. (54)] reduces 
to {sj — — 2v|) + 2ujVj\Agc\ = Ej, which is just the well-known energy of 
the Bogoljubov quasiparticle state 7]+|BCS), where 

jj^ = UjCj„ - avje"^c]_„ . (58) 
Thirdly, an example of a phase-breaking excitation is 

7]+7]-|BCS) = i-vje-"^ + Ujb]) [] {ui + Vje%l)\Ye.c) , (59) 
which has (Cjj)excited — —UiViUjVj and energy 2Ej. 

It should be appreciated, however, that BCS chose a grand- canonical con- 
struction purely for calculational convenience (as is made clear on p. 1180 of 
their original paper [101]): the trick of using a factorized form of commuting 
products in (47), at the cost of A^-indefiniteness, makes it brilliantly easy to 
determine the variational parameters Uj and Vj. In fact, BCS proposed them- 
selves to use the projection of |BCS) to fixed N as the actual variational 
ground state, namely [105] 

|PBCS) = Jd(f) e-'^^l[{uj+e^'^ vjb]) | Vac) (60) 

-J^(n»,)(i:^^])"'V-) (61) 

(PBCS for Projected BCS), which is of the general form of Eq. (41). In the 
bulk limit {d/A <^ 1), however, it is completely adequate to use |BCS): firstly, 
the relative error which its factorized form causes, by taking the occupation 
amphtude of level j to be independent of that of level i, scales like 1/N [105, 
pp. 150,163]; and secondly, the fluctuations in its particle number, (AA'^^)bcs = 
(-^^)bcs — N"^ = Y^ji'^UjVjY, are equal to nA/d in the bulk limit, in which 
the relative fluctuations (AA^^)bcs/^^ oc dA/e-p^ therefore vanish. Thus, bulk 
results obtained from |PBCS) or |BCS) are essentially identical. In fact, Braun 
[35,91] checked by explicit calculation that the functions (^1)0, ('1^|)pbcs and 
v| are practically indistinguishable even for d/A as large as 0.5 [see Sec. 5.3]. 
Signiflcant differences do develop between them once d/A increases past 0.5, 
however, as will be discussed in Sec. 5. 

To end this section, note that Fig. 12 offers a very simple intuitive picture 
for why pairing correlations weaken with increasing level spacing until, in 



50 



Anderson's words [45], "superconductivity is no longer is possible" when d > 
A: an increase in level spacing implies a decrease in the number of levels 
within A of e-p for which u^Vj differs significantly from zero, i.e. a decrease in 
the number of pairs with significant pairing correlations. This number, namely 
A/d, can roughly speaking be viewed as the "number of Cooper pairs" of the 
system, and when it becomes less than one, as in Fig. 12(c), it no longer makes 
sense to call the system "superconducting" . However, this should not be taken 
to imply that pairing correlations cease altogether in this regime; remnants of 
them do persist, in the form of fluctuations, up to arbitrarily large d/A, as 
will be discussed in detail in Sec. 5.3. 



4-4 Generalized variational BCS approach 

In the next several sections we review the generalized variational BCS ap- 
proach used by Braun et al. [35,23,91] to describe the paramagnetic break- 
down of superconductivity in nm-scale grains in a magnetic field. This theory 
produces theoretical excitation spectra that are in good qualitative agreement 
with the measurements of BRT shown in Fig. 10 and thereby yields the most 
direct confirmation available of the relevance to experiment of the discrete 
BCS model. Moreover, it sheds considerable light on how "superconductivity 
breaks down" (more precisely, how pairing correlations weaken) with increas- 
ing d and h: As mentioned in the previous paragraph, in grains with d ~ A 
(bulk gap), near the lower size limit [45] of observable superconductivity, the 
number of free-electron states with strong pairing correlations (those within 
A of e-p) is of order one. Thus, even in grains in which a spectral gap can 
still be observed, pairing correlations are expected to become so weak that 
they might be destroyed by the presence of a single unpaired electron [21]. 
This can be probed directly by turning on a magnetic field, since its Zeeman 
energy favors paramagnetic states with nonzero total spin. 

The theory reviewed below exploits analogies to thin films in a parallel mag- 
netic field [98,99], but explicitly takes account of the discreteness of the grain's 
spectrum. Since in RBT's experiments the temperature T = 50mK is much 
smaller than all other energy scales (d, A) , we shall neglect finite-temperature 
effects and set T = 0. In Sec. 4.4 the eigenenergies Sa of the grain's lowest- 
lying eigenstates \a) are calculated approximately using a generalized g.c. 
variational BCS approach that goes beyond standard mean-field theory by 
using a different pairing parameter A^ for each \a). The are then used to 
discuss various observable quantities, such as /i-dependent excitation spectra 
(Sec. 4.6), direct experimental evidence for the dominance of purely time- 
reversed states in the pairing interaction (Sec. 4.7), and various parity effects 
(Sec. 4.8). 
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Fig. 13. Cartoon of four typical variational states, labeled using the notation of 
Eq. (63) for (a-c) and Eq. (64) for (d). They represent (a) the even ground state 
|0); (b) the odd ground state |^); (c) the spin-| ground state ||); (d) a spin-| 
excited state ||, 2). The single-particle levels are drawn for h = 0, with the chemical 
potential half-way between levels and 1 for even systems (a), but exactly on level 
for odd ones (b,c,d). The ellipses joining states on the same level are meant to 
represent a "Cooper pair", and signify its being empty or doubly occupied with 
amplitude {uj + Vjbj); solid (dashed) ellipses are used for levels that would be 
completely filled (empty) in the absence of pairing correlations. 

The reasons for deciding to calculate the excitation spectra, despite their fixed- 
N nature, within a grand- canonical framework are as follows: Firstly, its sim- 
plicity. Secondly and perhaps most importantly, the exact eigenenergies have 
the general form = En + ^sih) [Eq. (42)], in which all h-dependence resides 
in the exactly known contribution Ssih) from the blocked levels. The choice of 
approximation scheme therefore only affects 8^ which determines the h — Q 
properties of the spectrum, such as the size of the zero-field spectral gap, 
etc., but not the qualitative features of the ^-dependence. In particular, this 
means that all of the analysis below could easily be "made exact" by simply 
replacing the g.c. approximations for 8^ by the exact values from Richardson's 
solution. However, this is expected to cause only slight quantitative differences, 
since, thirdly, canonical calculations (mentioned after Eq. (61) and discussed 
in Sec. 5) yield very similar results to g.c. ones as long as d/A < 0.5, which, 
by inspection of Fig. 10, does seem to be the case for the grain in question 
(the analysis of Sec. 4.6 yields d/A ~ 0.67). 



The generalized variational Ansatz 

The Zeeman term in the Hamiltonian of Eq. (37) favors states with a nonzero 
total z-component of the total spin, s — \ J2ja ^^]a^ja (henceforth simply 
called "spin"). Increasing h will thus eventually lead to a series of ground 
state changes to states with successively larger spins. In general, therefore, we 
are interested in pair-correlated states with nonzero spin, and in particular in 
their eigenenergies. Following Braun et al. [22,23,91], we now show how this 
can be calculated variationally, using the following general BCS Ansatz for a 
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state \s, B) with N = 2n + 2s electrons and a definite total spin s > (first 
introduced by Soloviev for application in nuclei [104]): 



k> 5) = n 4 h^^""^ + ^t''^) iVac). (62) 
ieB j 

If the spin is nonzero, it is built up by placing 2s unpaired spin-up electrons 
in a set S of 6 = 2s single-particle levels [cf. Eq. (40)] while the remaining 
single-particle levels have BCS-like amplitudes to be either empty (uf'^^) or 
doubly occupied by a pair (vj^'^^), with (u^^'^^y + {Vj^'^^y = 1. The product 

thus constitutes a g.c. approximation to the state of Eq. (41). The 
superscript U over products (and over sums below) indicates exclusion of the 
singly occupied levels in B, for which u^^'^\ v^^'^^ are not defined. 

More specifically, in a given spin-s sector of Hilbcrt space the following two 
types of specializations of Eq. (62) were studied in detail (p = 2s mod 2): 



s-l+p/2 U 

k)= n 4nK + ^j&])|Vac). (63) 

i=-s+p/2 j 

\s, k) = cf,_i+p/2+fc)+C(s_i+p/2)+|s) . (64) 

\s) is the spin-s state with the lowest energy, i.e. the "variational spin-s ground 
state", obtained by placing the 2s unpaired electrons as close as possible to ep 
[Fig. 13(b,c)], in order to minimize the kinetic energy cost of having more spin 
ups than downs. \s,k) is a particular type of excited spin-s state, obtained 
from |s) by moving one electron from its topmost occupied level (s — 1 +p/2) 
upwards by k units of d into a higher level {s — l+p/2 + k). These constructions 
are illustrated in Fig. 13, of which (a) and (b) represent the variational ground 
states of a grain with an even or odd number of electrons, respectively. 

The orthogonality of the wavefunctions, {s, B\s' , B') = Sss'Sbb', implies that 
the variational parameters Vj'^'^^ and Wj*'^^ must be found anew for each (s, B) 
(hence the superscript), by minimizing the variational "eigenenergies" 



£f^2^{h,d) = {s,B\H\s,B) (65) 
= -2sh + E(£, - /.) + E [2{ej - ii){v^t^^f + Ad(^f ^))^] (66) 

U 



u) >v) 



which we use as approximations to the exact eigenenergies £|'g^*(/i, d). The 
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term is not extensive and hence neglected in the bulk case where only effects 
proportional to the system volume are of interest. Here it is retained, since 
in ultrasmall systems it is non- negligible (but not dominant either) [22,23]. 

Q£BCS 

Solving the variational conditions ^-^Sy = in standard BCS fashion yields 

(.f^V = i(i - + Kb]'^') , 6- ^ - - xdiv^-'^r , (67) 



where the "pairing parameter" ^ is defined by the relation 



u -1 u 



^ Ad E v^^^^^lf'^ , or 1 = E ' (68) 



which in the limit d/A — > reduces to the standard bulk T = gap equation. 
Note that it is /i- independent, because it involves only unblocked levels j e U, 
which are populated by pairs with zero total Zeeman energy. Note also that in 
Eq. (67) the \d{v'f'^^y' shift in usually neglected because it simply renor- 
malizes the bare energies, is retained, since for large d it somewhat increases 
the effective level spacing near e-p (and its neglect turns out to produce a 
significant upward shift in the £f^{h, dys, which one is trying to minimize). 

The chemical potential /i is fixed by requiring that 



2n + 2s = {s, B\N\s, B)^2s + 2 Y^ivj^'^^f- (69) 

j 

In contrast to conventional BCS theory, the pairing parameter A^ ^ can in 
general not be interpreted as an energy gap and is not an observable. It should 
be viewed simply as a mathematical auxiliary quantity which was introduced 
to conveniently solve the variational conditions. However, by parameterizing 
Vj^'^'' and uf'^\ A^ b does serve as a measure of the pairing correlations 
present in \s,B): for vanishing A^^ the latter reduces to an uncorrelated 
paramagnetic state |s, B)q with spin s and energy £^°^, namely 

|s,S)o=n4+n&]|0), with £% = ,{s,B\H\s,B)^, (70) 
and the condensation energy Eg°^^ = Sf ^^ — ^ of B) reduces to zero. 



54 




Level spacing M Level spacing ^^/A 



Fig. 14. Properties of (a,b) spin-s ground states \s) = \s,0) [Eq. (63)] and (c,d) 
spin-^ excited states \^,k) (for k = 0, ...,4) [Eq. (64)], as functions of d/A 
(i.e. decreasing grain size), calculated for A = 0.194. (a) The pairing parameters 
As{d)/A, which vanish at critical level spacings 

d^cs of 2.36, 0.77, 0.44, 0.31,. . . 
for s = 0, 1/2, 1,3/2, . . ., respectively, (c) The pairing parameters Ai/2,fe) together 
with their k ^ oo limit, which equals Aq (dashed line). (b,d) show relative energy 
densities (since normalized by d/A^ oc Vol~^) at /i = for both correlated (solid) 
and uncorrelated (dashed) states, the latter obtained by setting A^^^ = in the 
former, (b) shows [Sf^^ - S^^^)d/A^ (solid) and {S^ - Eli^djA^ (dashed), the 
energy differences of \s) and |s)o relative to the uncorrelated spin-p/2 Fermi sea 
|p/2)o. (d) shows (ffCS_ _ £^^^)d/A^ (sohd) and (f 0/2,fe " ^i/i)^/^^ (dashed), the 
energy differences of \\-,k) and |^,A;)o relative to the uncorrelated spin-^ ground 
state ||,0)o = |^)o- Solid and dashed lines meet at the critical level spacing d^'^^ at 
which As^k becomes and the condensation energy E'^^^^ = Sf^^ ^ ^sk vanishes. 

4.4-2 General numerical solution - illustration of the blocking effect 

The simultaneous solution of Eqs. (67), (68) and (69) is a straightforward nu- 
merical exercise which Braun and von Delft performed [22,23], for the sake 
of "numerical consistency", without further approximations. (Analytical so- 
lutions can be found only in the limits d <^ A and d ^ ^s,b, see App. A 
of [23].) The numerical results are summarized in Fig. 14, which shows the 
pairing parameters A^^^ [Figs. 14(a,c)] and energies Sf jp [Figs. 14(b,d), solid 
hues] of some selected variational states |s, -B), as well as the energies Sg^^ 
of the corresponding uncorrelated states \s,B)o [Figs. 14(b,d), dashed lines]; 
both go 

Q are plotted relative to the energy E^^j^ of the uncorrelated 
spin-p Fermi sea |p/2). The results have a number of salient features: 

(i) In the bulk limit d/A — > 0, all of the pairing parameters A^^^ reduce to 
A, as expected, and the energy differences Sf^ — E^i^ between the corre- 
lated states |s, -B) and the uncorrelated Fermi sea |p/2) reduce to — |A^/o? = 
— ^A^(£f)A^, which is the standard bulk result for the condensation energy. 

(ii) Each Ag^B in Figs. 14(a,c) decreases with increasing d. This reflects the 
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fact that with increasing d, the number of pair-correlated states within A of Sp 
decreases [cf. Fig. 12 and the last paragraph of Sec. 4.3.5], so that the amount 
of pairing correlations, for which ^ is a measure, decreases too. 

(iii) Each Ag^B vanishes abruptly at a critical level spacing (whose precise 
numerical value depends sensitively on model assumptions such as the value 
of A and the use of uniformly-spaced levels [24]). For d > d^%^ no pairing 
correlations exist at this level of approximation, so that that the condensation 
energy E^"^"^ (difference between solid and dashed lines) vanishes and the solid 
and dashed lines in Figs. 14(b,d) meet. 

(iv) In Figs. 14(a), the pairing parameters A^ for the spin-s ground states 
decrease rapidly with increasing s at fixed d (and d^^^ < d^P^ if s > s'). 
[This is a generalization of a parity effect discussed by von Delft et al. [21], 
who studied only ground state pairing correlations and found that these are 
weaker in odd (s = 1/2) grains than in even {s = 0) grains, A1/2 < Aq, cf. 
Sec. 4.8.] This tendency is a direct consequence of the blocking effect described 
in Sec. 4.2.3 and is independent of model details: larger s means more unpaired 
electrons, more terms missing from the sum , less correlated pairs and hence 
smaller Ag^B- 

(v) As d, increases the blocking effect described in (iv) becomes stronger, i.e. 
the difference between the various A^ for different s becomes more pronounced, 
since then the relative weight of each term missing in the sum increases. 
The blocking effect is most dramatic in the regime d/A e [0.77, 2.36] in which 
Ao 7^ but A^^o = 0. This is a regime of "minimal superconductivity" [22,23], 
in the sense that all pairing correlations that still exist in the even variational 
ground state |0) (since Aq 7^ 0) are completely destroyed by the addition of a 
single electron or the flipping of a single spin (since A^^o = 0). 

(vi) Considering the spin-| excited states of Figs. 14(b,d), one finds 
that the larger k, the longer the pairing correlations survive with increasing d: 
the critical spacings dfj^j^ increase with k, approaching the value d^'^^ of the 
spin-0 case as /c — > 00; correspondingly, the larger k, the larger the d-value at 
which the condensation energies £^2% [differences between sohd and dashed 
lines in Fig. 14(d)] vanish. The intuitive reason why the amount of pairing 
correlations in an excited \s,k) increases with k is of course quite simple: 
the further the unpaired electron sits from the Fermi surface where pairing 
correlations are strongest, the less it disrupts the latter (since UkVk becomes 
very small for large k, see Fig. 12). In fact, the state ||, A; — 00) will have just 
about the same amount of pairing correlations as the even ground state ]0) 
(Ai/2,jt-^oo ^ Ao). 

(vii) Similar effects hold for excited states in other spin sectors (not shown): 
The higher the excitation, the larger the pairing parameter As,b- However, the 



56 



concomittant gain in correlation energy is always less than the kinetic-energy 
cost of having an unpaired electron far from ep. 

(viii) The strong dependence of ^ on s and d for d > A illustrates why in 
this regime a conventional mean-field treatment is no longer sufficient: the sys- 
tem cannot be characterized by a single pairing parameter, since the amount 
of pairing correlations vary from state to state, each of which is character- 
ized by its own pairing parameter. Instead, the present variational approach 
is, roughly speaking, equivalent to a doing a separate mean-field calculation 
for each new choice U of unblocked levels within the Fock space spanned by 
them (i.e. replacing bj {bj — (bj)} + (bj) and neglecting terms quadratic in 
the fluctuations {bj — (bj)}). Indeed, the behavior of As^sid) near d^^^ has 

the standard mean-field form ^1 — d/d^^^, as can be shown analytically [23, 
App. A]. 

To summarize: pairing correlations decrease with increasing d and s and de- 
creasing k. These features survive also in more accurate canonical calculations. 
This is not the case, however, for the abrupt vanishing of Ag^B at dfW^ which 
signals the breakdown of the g.c. approach once d becomes of order A: canon- 
ical methods show that, regardless how large d becomes, some remnants of 
pairing correlations survive and the pairing parameters (A^ B)can do not van- 
ish [Sec. 5], in accordance with the rule of thumb that "in a finite system no 
abrupt phase transition can occur between a zero and nonzero order parame- 
ter." 



4-5 Softening of the H -induced transition to a paramagnetic state 

Since states with nonzero spin are favored by the Zeeman energy but have 
smaller correlation energy due to the blocking effect, a competition arises 
between Zeeman energy and correlation energy. The manifestations of the 
blocking effect can thus be probed by turning on a magnetic fleld; if it becomes 
large enough to enforce a large spin, excessive blocking will destroy all pairing 
correlations. 

The situation is analogous to ultra-thin films in a parallel magnetic field 
[98,99], where orbital diamagnetism is negligible for geometrical reasons and 
superconductivity is destroyed at sufficiently large h by Pauli paramagnetism. 
This occurs via a first order transition to a paramagnetic state, as predicted by 
Clogston and Chandrasekhar (CC) [110,111] by the following argument (for 
bulk systems): A pure Pauli paramagnet chooses its spin s such that the sum of 
the kinetic and Zeeman energies, s"^ /M{e-p) — 2hs, is minimized, and hence has 
spin s = hj\f{ep) and ground state energy —h'^J\f{e-p)- When this energy drops 
below the bulk correlation energy — |A^AA(£f) of the superconducting ground 
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Fig. 15. Thin films in a magnetic field (Figs. 11 and 12 of [98]). (a) Tunneling con- 
ductance from a normal metal through a tunnel barrier into a thin superconducting 
Al film, as function of voltage, for several magnetic fields labeled in increasing or- 
der "a" to "f". The conductance reflects the BCS quasiparticle density of states, 
whose single peak (for a given sign of V) for H = splits into two separate peaks 
for H ^ 0, corresponding to the Zeeman energy difference between quasiparticles 
with spin up and down, (b) Voltage corresponding to the maxima of spin-up and 
spin-down density of states as a function of magnetic field. At the critical field -f^cc 
superconductivity is destroyed and the tunneling threshold drops abruptly to zero. 

state, which happens at the critical field hcc = ^/V^, a transition will occur 
from the superconducting to the paramagnetic ground state. The transition is 
first-order, since the change in spin, from to scc = ^ccA/'(£f) = A/((iv^), 
is macroscopically large {Af^ep) = 1/d oc Vol). 

This transition has been directly observed by Meservey and Tedrow [98,99] in 
ultra-thin (5nm) superconducting Al films (A = 0.38meV), whose density of 
states [Fig. 15(a)] they measured via the tunnel conductance through an oxide 
layer between a normal metal and the film. They found that in a magnetic field 
the BCS quasiparticle peak splits up into two subpeaks, separated in energy 
by [Fig. 15(b)], which simply refiects the Zeeman splitting of quasipar- 

ticles statesP^ with spin up or down (and g = 2). Remarkably, the tunneling 
threshold abruptly dropped to zero at a field of 4.7 T [Fig. 15(b)], which 
they associated with the field Hqq at which the phase transition from the 
superconducting to the paramagnetic ground state occurs. Indeed, Fig. 15(b) 
demonstrates clearly that the transition to the normal state is first order: the 
mean of the spin-up and spin-down spectral gaps, i.e. the pairing parameter 
A, is constant until the critical field Hqc is reached, at which it abruptly drops 
to zero. 

For the case of isolated ultrasmall grains, the above picture of the transition 
needs to be rethought in two respects due to the discreteness of the electronic 
spectrum: Firstly, the spin must be treated as a discrete (instead of continuous) 
variable, whose changes with increasing h can only take on (parity-conserving) 

Recall that the BCS quasiparticles 7]^ = n,cj ^ — avjC- have well-defined spins. 
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integer values. Secondly, one needs to consider more carefully the possibility 
of /i-induced transitions to nonzero spin states that are still pair- correlated 
(instead of being purely paramagnetic), such as the variational states \s,B) 
discussed above. (In the bulk case, it is obvious that such states play no role: 
the lowest pair-correlated state with nonzero spin obtainable from the ground 
state by spin flips is a two-quasiparticlc state, costing energy 2 A — 2h; when h 
is increased from 0, the paramagnetic transition at hcc = A/y^ thus occurs 
before a transition to this state, which would require h = A, can occur.) 

Quite generally, the effect of increasing h from can be analyzed as follows: 
At given d and h, the grain's ground state is the lowest-energy state among 
all possible spin-s ground states \s) having the correct parity p = 2s mod 2. 
Since Ss{h,d) = Ss{^,d) — 2hs, level crossings occur with increasing h, with 
£s' dropping below Ss at the level crossing field 



Therefore, as h is slowly turned on from zero with initial ground state |so = 
p/2), a cascade of successive ground-state changes (GSCs) to new ground 
states |s2), . . .will occur at the fields hg^^s^i hg^^sii ■ ■ - Let us denote this 
cascade by (sq, Si); (si, S2); ■ ■ ■; for each of its GSCs the corresponding level- 
crossing field hs^s'{d) is shown in Fig. 16. Generalizing CC's critical field to 
nonzero d, let us denote the (parity-dependent) field at which the first transi- 
tion {sq, Si) occurs by hcc{d,p) = hg^^siid), which simply is the lower envelope 
of the level-crossing fields hsQ,si in Fig- 16 (shown as bold solid and dashed 
lines for sq = and sq = |, respectively). In the limit 0? — it is numerically 
found to reduce to the Clogston-Chandrasekhar value, i.e. /icc(0,p) = A/y/2, 
as expected. 

In general, the order in which the GSCs occur with increasing h within a 
cascade (i.e. the order of hs g/ lines encountered when moving vertically up- 
ward in Fig. 16) depends sensitively on d, and an infinite number of distinct 
regimes (cascades) I, II, III, . . . can be distinguished: Starting at large d we find 
the typical normal-grain behavior (0, 1); (1, 2); (2, 3); . . . for even grains and 
(|, I); (|, I); . . . for odd grains, with /io,i < (or >) | in regimes I (or II). In 
regimes III and IV of somewhat smaller d, the order of GSCs is (0, 2); (2, 3); . . . 
and (|, |); (|, |); . . ., etc, i.e. the spin Si attained after the first GSC (sq, Si) 
has increased to 2 in the even case. This illustrates a general trend: the spin 
si{d) aft,er the first transition increases with decreasing d and becomes macro- 
scopically large in the d — > limits where Si = hcc/d — A/{dy/2), as explained 
in recounting CC's argument above. 

Furthermore, it turns out that Ag^^i^d) = and therefore i^^*"^ = for all d, 
implying that after the first GSC the new ground state |si) is always (not 
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Fig. 16. d-dependence of the level-crossing fields hs,s'{d)/A [Eq. (71)] at which S^*^^ 
drops below Sf^^ with increasing h. Only those level crossing fields are shown that 
belong to the cascade of (fixed-N) ground state changes (GSCs) (so,si); (si,S2); 
. . . that occur as h increases from at given d. Solid (dashed) lines are used for even 
(odd) grains with integer (half-integer) spins, and some are labeled by the associated 
GSC {s,s'). (In contrast, in Fig. 17 the N-changing tunneling transitions in are 
labeled by |sj) — |s/).) The size \ASsj^j' — ASg^jl of the first jump (occuring at the 
level-crossing field hcciPid) = /iso.si) in the lowest line of the tunneling spectra of 
Fig. 17 is shown by the lowest two (jagged) curves (solid for e ^ o and dashed for 

^ e tunneling spectra), which both approach the CC value 1 — l/\/2 as d — >■ 0. 

only in CC's bulk limit) a purely paramagnetic state, i.e. without any pairing 
correlations in the g.c. framework (canonical calculations would yield some 
weak remnant pairing correlations in the form of fluctuations). In this regard, 
CC's picture of the transition remains valid throughout as d is increased: 
at hccid,p), a transition occurs from the superconducting ground state to a 
paramagnetic, uncorrelated state |si)o, the transition being first-order in the 
sense that Ag^^d) = 0; however, the first-order transition is "softened" with 
increasing d, in the sense that the size of the spin change, si — Sq, decreases 
from being macroscopically large in the hulk to being equal 1 atd':$> A ( regimes 

1 and II). 

To conclude this section, we mention that the above analysis of the paramag- 
netic breakdown of superconductivity has recently been generalized to finite 
temperatures [32] , using the so-called static path approximation [explained in 
Sec. 5.5.3] to treat fiuctuation effects properly. 

4-6 Excitation spectrum in a magnetic field 

In this section we compare the theoretical tunneling spectra for a grain coupled 
to leads, calculated as functions of h and d [22,23,91], and compare these to 
RBT's measurements of Fig. 10. 
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The form of the tunnehng spectrum depends in a distinct way on the specific 
choice of level spacing d and on the electron number parity p of the final 
states I/) of the bottleneck tunneling processes \i) \ f) (or \a') — >• \a.) in the 
notation of Sec. 2.3). However, for the uniformly spaced £j-levels used here, 
particle-hole symmetry ensures that there is no difference between electron 
addition or removal spectra I^Afipi) I/a^)- To calculate the spectrum for 
given d andp, Braun et al. [22,23,91] proceeded as follows: they first analyzed 
at each magnetic field h which tunneling processes \i) — >■ |/) are possible, then 
calculated the corresponding tunneling energy thresholds ASif{h) = £f{h) — 
Si{h) [cf. Eq. (26)] and plotted A£^jj^(/i) — A£^iiiin(0) as functions of h for various 
combinations of \i) and |/), each of which gives a line in the spectrum. Since 
the selection rule sj — Si = ±1/2 holds, only slopes of ±1 can occur. The 
reason for subtracting A£inin(0), the h = threshold energy cost for the 
first (lowest-lying) transition, is that in experiment, this energy depends on 
Vg and hence yields no significant information, as explained in Sec. 2.4.6. 
Neglecting nonequilibrium effects [10-13] (which were minimized in the present 
experiment by tuning Vg, cf. Sec. 2.4.4, and which are discussed in Sec. 6.2), 
the initial state is always taken to be the ground state of a given spin-s sector. 
The appropriate Si{h,d) must be determined from Fig. 16. 

Fig. 17 shows four typical examples of such theoretical tunneling spectra, 
with some lines labeled by the corresponding \i) — > |/) transitions. Whenever 
h passes through one of the level-crossing fields hg-^s , of Eq. (71), the grain 
experiences a ground state change (s^, Sj') , at which the set of allowed tunneling 
transitions changes from — > to — > Therefore, at hs^^s^, 

one set of lines in the tunneling spectrum ends and another begins, producing 
kinks or discontinuities. A kink occurs if one of the new final states coincides 
with one of the old ones, |/') = |/), meaning that it can be reached from both 
|sj) and \si') [i.e. sj — Si = —{sj — Sj/)], in which case ASif{h) and ASif'{h) 
have slopes of opposite sign. However, for most lines this is not the case, so 
that at hg.^s'^ the line \si) — >■ |/) simply ends while new lines — > (/') begin. 
This results in discontinuities (or "jumps") in the spectrum at hg^^s'. of size 
{ASi'f — ASif){hs.^s'^), unless by chance some other final state |/') happens to 
exist for which this difference equals zero. 

Since the order in which the GSCs (sj, Sj') occur as functions of increasing h 
depend on d and p, as indicated by the distinct regimes I, II, III, ... in Fig. 16, 
one finds a distinct kind of tunneling spectrum for each regime, differing from 
the others in the positions of its jumps and kinks. In regime I, where the 
order of occurrence of GSCs with increasing h is (0, 1); (|, |); (1, 2); (|, |); . . ., 
there are no discontinuities in the evolution of the lowest line [see Fig. 17(a)]. 
For example, for the e ^ a spectrum, the lowest |0) 1 1/2) line changes 
continuously to |1) |l/2) at h^ i, since |s/ — = 1/2. However, in all other 
regimes the first change in ground state spin (at /io,si from to Si) is > 1, 
implying a jump (though possibly small) in all e — > o lines, as illustrated by 
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Fig. 17. The theoretical odd-to-even and even-to-odd tunnehng spectra 
{ASif — A£^min(0))/A predicted for an ultrasmall superconducting grain as a func- 
tion of magnetic field h, for two different level spacings: (a) d = 0.67A and (b) 
d = 0.34A (corresponding to regimes I and III of Fig. 16, respectively). Some lines 
are labeled by the corresponding Sj s[ tunneling transitions. Not all possible 
higher lines (corresponding to excited final states \s,j)) are shown. Vertical dashed 
lines indicate those level-crossing fields hg^s' [Eq- (71)] at which kinks or jumps 
occur, with /io,i < /Ji/2,3/2 < ^1,2 < ^3/2,5/2 in (a) and /ii/2,3/2 < ^0,2 < /i2,3 in (b). 

Fig. 17(b). 

The jump's magnitude for the tunneling thresholds, i.e. the lowest e — > o and 
o — > e lines, is shown as function of d in the lower part of Fig. 16. It starts 
at d = from the CC value A(l - l/^/2) measured for thin Al films [98,99], 
and with increasing d decreases to (non-monotonically, due to the discrete 
spectrum). This decrease of the size of the jump in the tunneling threshold 
reflects the fact, discussed in Sec. 4.5, that the change in spin at the first 
ground state change {sq, si) decreases with increasing d (as si — sq ~ hcc/d), 
and signals the softening of the first-order superconducting-to-paramagnetic 
transition. 



The fact that the measured tunneling thresholds in Fig. 10 show no jumps 
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at all, which might at first seem surprising when contrasted to the thresh- 
old jumps seen at hcc in Fig- 15 for thin films in a parallel field [98,99], can 
therefore naturally be explained [22,23] by assuming the grain to lie in the 
"minimal superconductivity" regime I of Fig. 16 (where the jump size pre- 
dicted in Fig. 16 is zero). Indeed, the overall evolution (i.e. order and position 
of kinks, etc.) of the lowest lines of Fig. 10 qualitatively agrees with those of 
a regime I tunneling spectrum, Fig. 17(a). This important result rather con- 
vincingly establishes the phenomenological success of the discrete BCS model. 
It also allows one to deduce the following values for the level-crossing fields 
Hs^^s[ (indicated by vertical dashed lines in Figs. 10 and 17): ifo,i = 4T, 
-f^i/2,3/2 = 4.25T, Hi^2 = 5.25T and -^3/2,5/2 = 6.5T. As corresponding uncer- 
tainties wc take AHg. g'. — 0.13T, which is half the H resolution of 0.25T used 
in experiment. 

By combining the above Hg^^s' values with Fig. 16, some of the grain's less- 
well-known parameters can be determined somewhat more precisely: 

(i) To estimate the grain's "bulk Hcc" , note that since -ffi/2,3/2/-f^o,i — 1-06, 
this grain lies just to the right of the boundary between regions II and I 
in Fig. 16 where d/A ~ 0.63, at which we have /io,i/^cc — 0.95, so that 
Hcc = ^0,1/0.95 ~ 4.2 T. This is quite close to the value Hcc - 4.7 T 
found experimentally [98,99] in thin films in a parallel field, confirming our 
expectation that these correspond to the "bulk limit" of ultrasmall grains 
as far as paramagnetism is concerned. 

(ii) The grain's corresponding bulk gap is A = \/2^bHcc — 0.34 meV, implying 
a coupling constant of A = 0.189 [by Eq. (39)]. A posteriori, these values can 
be regarded as being more appropriate for the present grain than the choices 
A = 0.38 meV and A = 0.194 made in Sec. 4.4.2, though the differences are 
evidently not significant (12% for A and 3% for A). 

(iii) The mean level spacing implied by (i/A ~ 0.63 is d ~ 0.21 mcV. The crude 
volume-based value d ~ 0.45 mcV cited in the caption of Fig. 2 thus seems 
to have been an overestimate. It would be useful if this determination of d 
could be checked via an independent accurate experimental determination 
of d directly from the spacing of lines in the excitation spectrum. Regret- 
tably, this is not possible: the measured levels are shifted together by pairing 
interactions, implying that their spacing does not refiect the mean indepen- 
dent-electron level spacing d. Nevertheless, note that the measured spacing 
of 0.05 meV between the lowest two states of the odd grain agrees quite 

well with the crude BCS estimate \J K'^ + — A [cf. Eq. (56)], which gives 
0.06 meV when evaluated for d — 0.21 meV and A = 0.34 meV. 

The higher lines plotted in Fig. 17 correspond to transitions into spin-sj state 
of the form \sf,k) [cf. Eq. (64) and Fig. 13(d)] (for simplicity these were the 
only ones considered in [22,23,91], though in general others are expected to 
occur too). The jumps in these lines, e.g. in Fig. 17(a) at /ii,2, occur whenever 
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the two final excited states \sf,kf) and |s//,/c//) before and after the GSC at 
hs^^s'. have different correlation energies. (Recall that the correlation energy 
of an excited state \sf,Bf) can be nonzero even if that of the corresponding 
ground state \sf) is zero, since the former's unpaired electrons are further away 
from £f, so that ^sf,Bf > A^^,, see point (vi) of Sec. 4.4.) Experimentally, these 
jumps have not been observed. This may be because up-moving resonances 
lose amplitude and are difficult to follow [10] with increasing h, or because the 
widths of the excited resonances (~ 0.13A) limit energy resolution [11-13]. 

For somewhat larger grains, the present theory predicts jumps even in the 
lowest line, as illustrated in Fig. 17(b). It remains to be investigated, though, 
whether orbital effects, which rapidly increase with the grain size, would not 
smooth out such jumps. 

To conclude this section, we emphasize once again that more than qualitative 
agreement between theory and experiment can not be expected, since both 
the model and our variational treatment thereof are very crude: the model ne- 
glects, for instance, fluctuations in level spacing and in pair-coupling constants, 
and the g.c. wave functions become unreliable for d/A> 0.5. Furthermore, we 
neglected nonequilibrium effects in the tunneling process and assumed equal 
tunneling matrix elements for all processes. In reality, though, random vari- 
ations of tunneling matrix elements could suppress some tunneling processes 
which would otherwise be expected theoretically. 

4-7 Time-reversal symmetry 

In this section we argue that RBT's spectra give direct support for the domi- 
nance of purely time-reversed states in the pairing interaction, implying that 
the sufficiency of using only a reduced BCS Hamiltonian, well established for 
bulk systems and dirty superconductors, holds for ultrasmall grains, too. 

When defining the discrete BCS model in Eq. (37), we adopted a reduced BCS 
Hamiltonian, in analogy to that conventionally used for macroscopic systems. 
In doing so, we neglected interaction terms of the form 

-dJ2 Kh j, i', j')cl+c]_Ci'-Cj'+ (72) 

iji'j' 

between non-time- reversed pairs cj^c]_, following Anderson's argument [45] 
that for a short-ranged interaction, the matrix elements involving time-reversed 
states much larger than all others, since their orbital wavefunctions 

interfere constructively. (This argument can be substantiated using l/gdim as 
small parameter [13], where g^im oi Eq. (4) is the grain's dimensionless conduc- 
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tance, see Sec. 6.1.3.) Interestingly, the experimental results of RBT provide 
strikingly direct support for the correctness of neglecting interactions between 
non-time-reversed pairs of the from (72) at /i = 0: Suppose the opposite, 
namely that the matrix elements A(j -|- k,j,j' + k',j') were all roughly equal 
to A for a finite range of k- and A;'- values [instead of being neghgible for k or 
k' 7^ 0, as assumed in the reduced BCS Hamiltonian (37)]. Then for 2s < k, 
one could construct a spin-s state \s)' with manifestly lower energy [S') than 
that (S) of the state \s) of Eq. (63), namely 

imin + 2.S-l oo 

.n 4 n («5^'^+^jX-+2sHc]_)|Vac), (73) 

^ Jmin 3 Jmin 

where jmin labels the lowest-lying of the interacting levels (i.e. the bottom 
of the interacting band). Whereas in \s) pairing correlations involve only 
time-reversed partners, in we have allowed correlations between non-time- 
reversed partners, while choosing the 2s impaired spin-up electrons that oc- 
cupy their levels with unit amplitude to sit at the band's bottom (see Fig. 18). 
To see that \s)' has lower energy than \s), 

= £;;^°"^+£:f < = , (74) 

we argue as follows: Firstly, S'g^ = S^, since the corresponding uncorrelated 
states |s)o and |s)o are identical (and given by Eq. (70), where B is the 
set of levels i = —s + p/2, . . . , s — 1 + p/2). Secondly, A^, = Aq (implying 
^/cond _ ^cond^^ bccausc the 2s unpaired electrons in |s)' sit at the band's 
bottom, i.e. so far away from ep that their blocking effect is negligible [cf. 
point (vi) of Sec. 4.4.2]. In contrast, A^ < Aq (implying > ^S°°'^), 

because the 2s unpaired electrons in |s) sit around and cause significant 
blocking. Thus the condensation energies satisfy E'^'^'^^'^ < £'^°'^'^(< 0), so that 
Eq. (74) holds, implying that |s)' would be a better variational ground state 
for the interaction (72) than \s). 

Now, the fact that E'^'^'^'^ = E^"""^ is independent of s means that hipping 
spins in |s)' does not cost correlation energy. Thus, the energy cost for turning 
|0)' into by flipping one spin is simply the kinetic energy cost d, implying 
a threshold field h'^i = d/2 [see Eq. (71)]; in contrast, the cost for turning 
|0) into |1), namely {Si — Sq), implies a threshold field /io,i = |(^i — ^o), 
which (in the regime d < A) is rather larger than d/2. [{£i — Sq) is the h — 
"spectral gap" between the first and second lines of Figs. 10(a) and 17(a)]. 
The fact that RBT's experiments [Fig. 10(b)] clearly show a threshold field 
/io,i significantly larger than d/2 shows that the actual spin-1 ground state 
which nature chooses is better approximated by |1) than by |1)', in spite of 
the fact that £[ < Si. Thus the premise of the above argument was wrong. 
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Fig. 18. Schematic representations of the non-time-reversed-pairing state |3/2)' de- 
fined in Eq. (73). The energies ej^h of the single-particle states \ j, =b) are drawn (a) 
for ^ = and (b) for 2h = 3d. We indicated schematically how non-time-reversed 
states are paired according to (uj + Vicy_^^._^cl_) in the BCS-like Ansatz (73), with 
solid or dashed ellipses encircling states that would be completely filled or empty 
in the absence of pairing correlations. 

and we can conclude that those terms in Eq. (72) not contained in the reduced 
BCS Hamiltonian can indeed be neglected, as done throughout in this review. 



4-8 Measurable consequences of the blocking effect: parity effects 

This section is devoted to various measurable manifestations of the blocking 
effect, in the form of parity effects, i.e. differences between a grain with an 
even or odd number of electrons. 



4-8.1 Bulk consequences of blocking 

The most obvious measurable manifestation of the blocking effect is the very 
existence of a spectral gap: "breaking a pair" and placing the two newly un- 
paired electrons in two singly-occupied levels costs a significant amount of 
correlation energy, because the unpaired electrons loose pairing energy them- 
selves and also disrupt the pairing correlations of the other pairs. This, of 
course, is already present in standard bulk mean-field BCS theory via the en- 
ergy cost of at least 2A involved in creating two quasiparticles, and is one of 
the hallmarks of superconductivity. 

In the context of ultrasmall grains, let us denote the pair-breaking energies for 
an even (odd) grain, i.e. the minimum energy cost per electron for breaking a 
pair by flipping a single spin at h = 0, by fie (f^o)- 

Jle = 1(^1 — ^o)h=0, = 1(^3/2 — ^l/2)/i=0 ■ (75) 
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The even pair-breaking gap Jig is of course strikingly visible in RBT's h = 
spectra as a large spectral gap for even grains [cf. Figs. 9 and 10; the lat- 
ter gives Qe = 0.26 meV]. Its presence is direct evidence for the existence 
of pairing correlations in the grain, which in that sense can still be called 
"superconducting" . 

In contrast, the odd pair-breaking gap Qo can not be obtained from h = 
spectra, since in an odd grain the lowest excitation does not involve breaking 
a pair, but simply exciting the unpaired electron, which does not require a 
correlation-induced gap to be overcome. To measure Qq, a finite field is needed: 
by Eq. (71), fie = ^o,i and flo = ^1/2,3/2, hence both spin-flip gaps are equal 
to level-crossing flelds that can be deduced from h ^ data, as explained 
in Sec. 4.6. For Fig. 10 this yields = 0.23 ± 0.01 meV and = 0.24 ± 
0.01 meV [a result further discussed in Sec. 4.8.4]. The reason that the Jig- 
value determined in this way is somewhat smaller than the above-mentioned 
0.26 meV determined at = is presumably that the experimental spectral 
lines are not perfectly linear in h (having a small /^^-contribution due to orbital 
diamagnetism, which should cause the spectroscopic gap to decrease faster 
with h than in our model). 

Another consequence of the blocking effect is that the condensation energies 
~ ~ ^'^^ even and odd grain differ: the unpaired electron 
of an odd grain weakens its pairing correlations relative to an even grain, 
so that El°2^ is less negative than Eq°^^. In the bulk hmit their difference 

approaches E^^jf - E^"""^ ^ A, the energy of a single quasiparticle. For large 

mesoscopic islands (with d/A. <S 1) this energy difference has indeed been 
directly observed: it causes a change from e- to 2e-periodicity in the gate- 
voltage dependence of Coulomb oscillations [62,84-88]. For ultrasmall grains, 
however, ground state energy differences are currently not directly measurable, 
due to experimental difficulties explained in Sec. 2.4.6. 

The parity effects discussed above survive in the bulk limit. Let us now turn 
to parity effects that result from even-odd differences in the d- dependence of 
various quantities. 

4-8.2 Parity- dependent pairing parameters 

As is evident from Fig. 14(a,b), not only the condensation energies Ep/2 are 
parity dependent; as soon as one leaves the bulk regime, the pairing param- 
eters Ap/2 become parity-dependent too, with Aq > A1/2. In the context of 
ultrasmall grains this was first emphasized by von Delft et al. [21], but it had 
been anticipated before by Janko, Smith and Ambegaokar [109] and Golubev 
and Zaikin [112], who had studied the first correction to the bulk limit, finding 
Aq — Ai/2 = (^/2 to leading order in d/A; and this result, in turn, had already 
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Fig. 19. The parity parameters (a) A^^^ of Matveev-Larkin [Eq. (76)] and (b) A^*" 
for tlie pair-breaking energies [Eq. (81)], as functions of d/A, calculated using the 
g.c. variational BCS approach of Sec. 4.4.1 (dashed lines), and Richardson's exact 
solution of Sec. 5.1 (solid lines). In (a), we also show the perturbative result for 
the uncorrelated Fermi sea, (A^'")pert = ^Ad (straight dash-dotted line); and the 
renormalized result (A^L)rcn (i/[2 ln(a(i/A)] of Eq. (80), in its range of validity 
d/A 3> 1 (dotted line). The parameter a = 1.35 is chosen to ensure quantitative 
agreement with the exact result in the limit dj Is.':^ 1. For a summary of the results 
of various other canonical calculations of Ap see Fig. 20. The inset of (a) shows 
the Dyson equation used to calculated the renormalized coupling A in Eq. (78). 

been published by Soloviev in the nuclear physics literature as long ago as 
1961 [104]. 

The g.c. results of Fig. 14(a), in particular the fact that the critical level 
spacing d^J"^ at which Ap/2 vanishes is smaller for odd than even grains 
(^V2^ < ^0^^), suggest that "pairing correlations break down sooner in odd 
than even grains" [21]. However, it should be remembered that the vanishing 
of Ap/2 signals the breakdown of the g.c. approach. A more accurate state- 
ment, that is born out by the canonical calculations reviewed in Sec. 5, is that 
the inequality £^1/2 > ^0 persists for arbitrarily large d (see Fig. 20 in Sec. 5.2), 
i.e. pairing correlations are always weaker for odd than even grains, although 
they never vanish altogether in either. 

Jf..8.3 Matveev-Larkin parity parameter 

To capture the difference between correlations in even and odd grains in terms 
of measurable quantities (which Ap/2 are not), Matveev and Larkin [28] pro- 
posed the parameter (sometimes called "pairing energy" in nuclear physics 
[49]) 

A^^ = £^+^ - ^{£^ + £^+^) (where N is even) , (76) 

i.e. the difference between the ground state energy of an odd grain and the 
mean of the ground state energies of the even grains obtained by removing 
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or adding one electron. Fig. 19(a) shows its behavior as function of d/A. In 
the bulk limit we have Sq ~ ^^o^^^ and Ap^^ ~ A, which is simply the energy 
cost for having an unpaired electron on the odd grain. With increasing rf/A, 
this energy cost decreases since pairing correlations get weaker, hence Ap^ 
initially decreases. It begins to increase again for > A, since then pairing 
correlations are so weak that the behavior of Ap^ is governed by the "self- 
energy" of the one extra pair in Eq'^'^ relative to £q . For example, in the g.c. 
variational BCS result for Ap^, namely 

(A^^)bcs = ^S? - + Ad/2 , (77) 



it is this self-energy which produces the Xd/2 contribution. 

A more careful calculation for the regime d ^ A was performed by Matveev 
and Larkin [28], whose considered the renormalization of A due to "pairing 
fluctuations" about the uncorrelated Fermi ground state \p/2)q. Summing up 
the leading logarithmic vertex corrections [113] [see inset of Fig. 19(a)], they 
obtained a renormalized coupling A given, with logarithmic accuracy, by 



A = A + Aci 
A 



u 



A 



27r [iu — {sj — n)] [—iu — {ej — n)] 
A 



A 



1- AdE 



1 



u 

j 2\s j-f,\ 



1 - Xlog{ujB/d) 



(78) 
(79) 



This result evidently is valid only if o? ^ UDe~^/'^ ~ A/2 (which, inciden- 
tally, is another way of seeing that d ^ A defines the crossover between the 
fluctuation-dominated and bulk regimes) . Matveev and Larkin concluded that 



(A^^)re„^ Ad/2 = d/(2 log d/A) for d>A. 



(80) 



This logarithmic renormalization is beyond the reach of the g.c. variational 
BCS method, but was confirmed using exact methods [33,91,39] (see Sec. 5.2). 
Its occurrence, in a regime that in g.c. variational calculations appears to be 
"uncorrelated" , can be regarded as the "first sign of pairing correlations" , in 
particular since, by Eq. (79), the interaction strength increases upon renor- 
malization only if the interaction is attractive (A < would imply |A| < |A|). 
The pairing fluctuations responsible for this renormalization will be discussed 
in more detail in Sees. 5.2 and 5.3. 

Unfortunately, Ap ^ is at present not measurable in ultrasmall grains, for the 
same experimental reasons as apply to £^1/2 — £0, see Section 2.4.6. 
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4-8.4 Parity effect for pairbreaking energies 

Braun and von Delft [35,23,91] discussed yet another parity effect, based on 

= - , (81) 

the difference between the pair-hreaking energies of an even and an odd grain 
[see Eq. (75)]. Fig. 19(b) shows its behavior as function of d. In the bulk 
limit Vie — Vlo — ^ and Ap*^ ~ 0. The most interesting feature of Ap*^ is 
that it initially becomes negative as d/ A increases; this occurs because in an 
odd grain pairing correlations are weaker and hence breaking a pair costs less 
energy than in an even grain. Ap*^ becomes positive again for d/A ~> 0.5, 
since then pairing correlations are so weak that Ap*^ is governed by the kinetic 
energy cost of flipping a spin, which is 2d for an odd grain but only d for an 
even grain. 

Ap^ is directly measurable in RBT's grains, via the level-crossing fields /io,i = 
Qe and /ii/2,3/2 = [Eq. (71)]. The measured values fig = 0.23 ± 0.01 meV 
and Qo = 0.24 ± 0.01 meV cited in Sec. 4.8.1 give a positive value of Ap^ = 

0. 1 meV, implying that the grain under study was too small to fall in the most 
interesting regime where Ap*^ is negative. Braun and von Delft suggested that 
the latter should be observable in a somewhat larger grain with /ii/2,3/2 < ^0,1, 

1. e. in Regime II of Fig. 16. (This suggestion assumes that despite the increased 
grain size, the complicating effect of orbital diamagnetism is still non-dominant 
in Regime II.) To look for negative Ap^ experimentally would thus require good 
control of the ratio d/A, i.e. grain size. This might be achievable if recently- 
reported new fabrication methods, which allow systematic control of grain 
sizes by using colloidal chemistry techniques [114,186,187], could be applied 
to Al grains. 
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5 Superconductivity: crossover from the bulk to the Hmit of a few 
electrons 



This section is devoted to the question: How do pairing correlations change 

when the size of a superconductor is decreased from the bulk to the limit of only 
a few electrons? In particular, we shall attempt to refine the answer given by 
Anderson [45], namely that superconductivity as we know it breaks down for 
d>A. 

First steps towards a more detailed answer were taken in the early 1970s by 
Strongin et al. [93] and by Miihlschlegel et al. [94] , who calculated the thermo- 
dynamic properties of ensembles small superconducting grains. Experimental 
realizations of such ensembles were, e.g., the granular films studied by Giaver 
and Zeller [80,81]. The interest of theorists was rekindled in 1995 by RBT's 
success in probing individual superconducting grains. Apart from motivating 
the phenomenological theory of Braun et al. reviewed in Sec. 4, these ex- 
periments also inspired a substantial and still growing number of theoretical 
studies [21-44] of how superconducting pairing correlations in such grains are 
affected by reducing the grains' size, or equivalently by increasing its mean 
level spacing d oc Vol~^ until it exceeds the bulk gap A. 

In the earliest of these, von Delft et al. studied the discrete BCS model of 
Sec. 4.2 within a parity-projected g.c. BCS approach [21] closely related to the 
variational BCS method of Sec. 4.4. Their g.c. results suggested that pairing 
correlations, as measured by the pairing parameter or the condensation energy, 
vanish abruptly once d exceeds a critical level spacing d^j^^ that depends on 
the parity (p = or 1) of the number of electrons on the grain, being smaller 
for odd grains {dfj"^ ~ 0.89A) than even grains 

were also found in a number of subsequent papers that used parity-projected 
g.c. methods to study the behavior of the BCS mean- field gap parameter Ag^ 
and related quantities as functions of level spacing [22-32], temperature and 
magnetic field. All these parity effects are consequences of the blocking effect 
(cf. Sec. 4.2.3): for odd grains, the unpaired electron somewhat disrupts the 
pairing correlations of the remaining paired ones, by reducing the phase space 
available for pair scattering. 

A series of more sophisticated canonical approaches [33-41] (summarized in 
Sec. 5.2) confirmed the parity dependence of pairing correlations, but estab- 
lished that the abrupt vanishing of pairing correlations at d^j^ is an artifact 
of g.c. treatments: pairing correlations do persist, in the form of so-called fluc- 
tuations, to arbitrarily large level spacings [28], and the crossover between the 
bulk superconducting (SC) regime {d <^ A) and the fluctuation-dominated 
(FD) regime {d S> A) is completely smooth [36-40]. Nevertheless, these two 
regimes are qualitatively very different [35-40]: the condensation energy, e.g.. 
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is an extensive function of volume in the former and almost intensive in the 
latter, and pairing correlations are quite strongly localized around the Fermi 
energy e-p, or more spread out in energy, respectively. Very recently, Di Lorenzo 
et al. [44] suggested that the remnant pairing correlations in the FD regime 
might be detectable via susceptibility measurements. 

Toward the end of 1998 and after the appearance of most of these works, 
R.W. Richardson pointed out [115] to their various authors that the discrete 
BCS Hamiltonian on which they are based actually has an exact solution, 
discovered by him in 1963 [46] (and independently by Gaudin in 1968 [116]). 
Richardson published his solution in the context of nuclear physics in a series 
of papers between 1963 and 1977 [46-54] which seem to have completely es- 
caped the attention of the condensed matter community. Very recently, the 
model was also shown to be integrable [117,118]. The revival of this remark- 
ably simple exact solution after such a long and undeserved period of neglect 
is perhaps one of the most important consequences of RBT's experimental 
breakthrough: Richardson's solution allows the elucidation and illustration by 
exact means of many important conceptual ingredients of the standard BCS 
theory of superconductivity, such as the nature of pairing correlations, the 
importance of phase coherence, the validity of using a mean-field approxima- 
tion and a grand-canonical formulation for bulk systems, and the limitations 
of the latter approaches for ultrasmall systems. Moreover, it allows the exact 
calculation of essentially all quantities of interest for ultrasmall grains. 

We shall therefore start this section by discussing the exact solution [Sec. 5.1]. 
We then summarize the other canonical approaches somewhat more briefly 
than they perhaps would have deserved had an exact solution not existed, and 
compare their results to those of the exact solution [Sec. 5.2]. Next we analyze 
the qualitative differences between the bulk and FD regimes [Sec. 5.3], then 
discuss the case of randomly (as opposed to uniformly) spaced energy levels 
Ej [Sec. 5.4], and finally discuss finite temperature parity effects [Sec. 5.5]. 
Throughout this section we set = 0, since canonical treatments make no 
reference to a chemical potential. 

5.1 Richardson's exact solution 

In this section we summarize some of the central results of Richardson's exact 
solution of the discrete BCS model. 

5.1.1 General eigenstates 

Consider N = 2n + b electrons, h of which are unpaired, as in Sec. 4.2.3. 
According to the general discussion there, the nontrivial aspect of solving 
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the model is finding the eigenenergies £n and corresponding eigenstates l^ri) 
[Eq. (43)] of the pair Hamiltonian [Eq. (44), in which we set = below] 



Hu = Y.{2ej5,j- Xd)blbj, (82) 



in the Hilbert space of all states containing exactly n pairs fc] = Cj+c^j- of 
electrons, where j runs over the set of all unblocked single-particle levels, 
U — I\B [I is the set of all interacting levels, B the set of all blocked levels] . 
In general, degenerate levels are allowed in /, but are to be distinguished by 
distinct j-labels, i.e. they have Si — Sj for i ^ j. 

Richardson showed that the sought-after eigenstates (with normalization 
{^n\^n) — 1) a-nd eigenenergies have the general form 



^ b\ 



|*n)=A^n^i|0>' ^n = E^- ^th St=E^-f^. (83) 



Here A/" is a normalization constant and the n parameters E^, — 1, . . . , n) 
are a solution of the set of n coupled algebraic equations 



— - E ET + E -E^ ET = , for 1/ = 1, . . . , n , (84) 

Ad ■ 2ej - ^,=l{^l,) - ^'^ 

which are to be solved (numerically, see App. B.2) subject to the restrictions 
E^^ El, ii H ^ V. Richardson originally derived this remarkably simple result 
by solving the Schrodinger equation for the wave- function . . . , jn) of 

Eq. (41). A simpler proof, also due to Richardson [119], may be found in 
Refs. [40,120] and in App. B; its strategy is to verify that {Hu — £„)|^„) = 
by simply commuting Hu past the Bl operators in (83). 

Below we shall always assume the e/s to be all distinct (the more general 
case that degeneracies are present is discussed by Gaudin [116]). Then it can 
be shown explicitly [116] that (i) the number of distinct solutions of Eq. (84) 
is equal to the dimension of the n-pair Hilbert space defined on the set of 
unblocked levels U, namely (^'^^^ , where Nu is the number of unblocked levels; 
and (ii) that the corresponding eigenstates (83) are mutually orthogonal to 
each other, thus forming an eigenbasis for this Hilbert space. This can easily 
be understood intuitively, since there exists a simple relation between the 
bare pair energies 26 j and the solutions of Eqs. (84): as A is reduced to 0, 
it follows by inspection that each solution {Ei, . . . , En} reduces smoothly 
to a certain set of n bare pair energies, say {25^^, . . . ,2ej^}; this particular 
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solution may thus be labeled by the indices ji, . . . ,jn-, and the corresponding 
eigenstate (83) written as |\1'„) = \ji, . . ■jn)u- By inspection, its A ^ limit is 
the state . . .jn)u,o = Il2=i ^'^^^ there is a one-to-one correspondence 

between the sets of all states {|ji, . . . ,jn)u} and {\ji, . . .jn)u,o}- But the latter 
constitute a complete eigenbasis for the n-pair Hilbert space defined on the 
set of unblocked levels U, thus the former do too. 



5.1.2 Ground state 

For a given set of blocked levels B, the lowest-lying of all states of 
the form (40), say |\E'„,i?)G, is obtained by using that particular solution 
lii) • • ■jn)u for which the total "pair energy" Sn takes its lowest possible value. 
The lowest-lying of all eigenstates with n pairs, b blocked levels and total spin 
s — 6/2, say |n, s)g with energy S^{ri), is that B)g for which the blocked 
levels in B all contain spin-up electrons and are all as close as possible to 
Ef, the Fermi energy of the uncorrelated A^-electron Fermi sea IFat). The E,y 
for the ground state \n, s)g coincide at A = with the lowest n energies 
2ej^ {u = 1, . . . ,n), and smoothly evolve toward (initially) lower values when 
A is turned on, a fact that can be exploited during the numerical solution of 
Eq. (84). As A is increased further, some of the E^^s become complex; however, 
they always occur in complex conjugate pairs, so that Sn remains real [51]. For 
details, see Ref. [51] and our App. B.2, where some algebraic transformations 
are introduced that render the equations less singular and hence simplify their 
numerical solution considerably. 



5.1.3 General comments 

Since the exact solution provides us with wave functions, it is in princi- 
ple straightforward to calculate arbitrary correlation functions of the form 
{"^nlbib^j ■ ■ ■ 6j,6^v|^„), by simply commuting all 6's to the right of all b^^s. How- 
ever, due to the hard-core boson commutation relations (46) of the 6's, the 
combinatorics is rather involved. Nevertheless, Richardson succeeded to derive 
[50] exphcit results for the normalization constant M of (83) and the occu- 
pation probabilities v'^ and correlators Cij of Eq. (49), which we summarize 
in App. B.3. The exact result for the Ci/s show that they are all positive, in 
agreement with the requirement (ii) formulated in Sec. 4.3.2. It is also natural 
to ask whether in the bulk limit (o? — > at fixed nd), the standard BCS re- 
sults can be extracted from the exact solution. Indeed they can, as Richardson 
showed in [54] (following unpublished work by Gaudin [116]), by interpreting 
the problem of solving the eigenvalue equations (84) for the E,y as a problem 
in two-dimensional electrostatics (see App. B.2). Exploiting this analogy, he 
showed that in the bulk limit, Eqs. (84) reduce to the well-known BCS equa- 
tions determining the gap and chemical potential at T = [Eqs. (68) and 
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(69)], and the ground state condensation energy EQ°^^{ri) [Eq. (45)] to its BCS 
result, namely — A^/2ci. 

Finally, let us mention that Cambiaggio, Rivas and Saraceno have recently 
shown that the discrete BCS model is integrable and have constructed ex- 
plicit expressions for all its constants of the motion [117]. The latter's relation 
to Richardson's solution was clarified by Sierra [118], who has also explored 
possible connections between the exact solution and conformal field theory. 
It would be an interesting challenge for mathematical physicists to try to ex- 
ploit this integrability to calculate finite-temperature properties exactly — 
although these can in principle be obtained from Richardson's solution by 
"simply" computing the partition function over all states, this is forbiddingly 
tedious in practice for large temperatures, since the eigenenergy of each state 
requires a separate (non-trivial) numerical calculation. 



5.2 Comparison of other canonical methods with the exact solution 

In this section we briefly mention the various canonical methods by which the 
discrete BCS model had been investigated prior to the revival of Richardson's 
exact solution in 1999. All of these studies used a half-filled band with fixed 
width 2ujd of N = 2n + p uniformly-spaced levels [i.e. ej = j d + {1 — p)d/2, 
with p = or 1, as in Eq. (38)], containing N electrons. Then the level 
spacing is d — 2u;b/N and the bulk gap is A = a;D/ sinh(l/A). Following [35], 
we take A = 0.224 throughout this section. To judge the quality of the various 
approaches, we compare in Fig. 20 the results which they yield with those 
from Richardson's solution, for the even and odd (s = 0,1/2) condensation 
energies and the Matveev-Larkin parity parameter Ap^^ [cf. Sec. 4.8.3]. 

In the notation of Sec. 5.1.2, these are given by 

Er''{n) = S^{n) - {F^\H\Fm) , (85) 
Ar(n) = £^/,{n) - [£^{n) + £^{n + l)]/2 . (86) 

Following the initial g.c. studies [21-24,28] of the discrete BCS model, the first 
purely canonical study was that of Mastellone, Falci and Fazio (MFF) [33], 
who used Lanczos exact diagonalization. Despite being limited to n < 12, they 
managed to reach reasonably small ratios of d/ A by using an ingenious scaling 
approach: for a given level spacing d, they increased the coupling constant A 
to about 0.5, thereby decreasing d/A = rf/cjo sinh(l/A) to values as small 
as 0.5. This allowed them to probe, coming from the few-electron side, a 
remarkably large part of the crossover to bulk limit. They found, i.a., that 
the condensation energies are negative for all d, showing that the system can 
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Fig. 20. (a) The even and odd {s = 0, 1/2) condensation energies Eg"^*^ of Eq. (85) 
[in units of A], calculated with BCS, PBCS and exact wave functions [39], as func- 
tions of d/A = 2sinh(l/A)/(2n + 2s), for A = 0.224. For comparison, the dotted 
line gives the "bulk" result E^"^^^ = -K'^/{2d). (b) Comparison [39] of the parity 
parameters Ap ^ [28] of Eq. (86) [in units of A] obtained by various authors: ML's 
analytical result (dotted lines) [A(l — d/2A) for d ^ A, and d/2log{ad/ A) for 
d 3> A, with a = 1.35 adjusted to give asymptotic agreement with the exact result]; 
grand-canonical BCS approach (dash-dotted line) [the naive perturbative result 
is continued to the origin]; PBCS approach (short-dashed line); Richardson's exact 
solution (solid line); exact diagonalization and scaling by MFF (open circles) and 
BH (long-dashed line). 



always gain energy by allowing pairing correlations, even for arbitrarily large 
d. 

Berger and Halperin (BH) [34] showed that almost identical results can be 
achieved with less than 6 pairs, thus significantly reducing the calculational 

effort involved, by first performing a "poor man's scaling" renormalization: 
they reduce the bandwidth from Ud ~ nd to, say, cud ~ nd (with n < 6) and 
incorporate the effect of the removed levels by using a renormalized coupling 
constant. 
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(87) 



The reduced system is then diagonahzed exactly. Note that the renormal- 
ization of Matveev and Larkin [Eq. (79)] corresponds to taking ujd ~ o? in 
Eq. (87), i.e. to integrating out the entire band. Also note that the renormal- 
ization prescription of (87) has the property that it would leave the bulk gap 
invariant in the limit d/A 0, for which Eq. (87) would imply ivBe'^^^ ~ 
ubc'^/^ ~ A. 

To access larger values of n, Braun and von Delft [35] used a fixed-n pro- 
jected BCS approach (PECS), in which BCS-likc variational wavefunctions 
are projected to fixed particle number, as in Eq. (60). The projection integrals 
occurring in Eq. (60) were evaluated numerically for n < 600, using tricks 
developed in the nuclear physics literature by Bayman [121], Dietrich, Mang 
and Pradal [122] and Ma and Rasmussen [123], and summarized in part in the 
book of Ring and Schuck [106] . (A much simpler way of dealing with the pro- 
jection, using recursion relations, was recently found by Dukelsky and Sierra 
[37].) The PBCS method gives condensation energies that (i) are significantly 
lower than the grand-canonical ones [see Fig. 20], thus the projection much 
improves the variational Ansatz, and that (ii) are negative for all d, confirm- 
ing that the abrupt vanishing of the g.c. condensation energies is indeed an 
artifact of the g.c. treatment. The PBCS method is able to fully recover the 
bulk limit, but the crossover is not completely smooth and shows a remnant 
of the g.c. breakdown of pairing correlations: the li-dependence of the conden- 
sation energy (i^cond-jPBCS changes rather abruptly [kinks in the short-dashed 
lines in Fig. 20(a)] from being extensive (~ 1/d) to being practically intensive 
(almost d independent). 

It should be mentioned here that a generalization of the PBCS method to 
finite temperatures has been worked out by Essebag and Egido in the context 
of nuclear physics [124]. The PBCS method has recently also been applied to 
the attractive Hubbard model in one dimension by Tanaka and Marsiglio [43], 
who found even-odd and super-even effects. The latter consist of differences 
between the number of pairs being equal to n = 2m or 2m -|- 1, and arise 
if boundary conditions are used that produce doubly-degenerate levels {e^^ — 
e_^) near the Fermi surface. 

Dukelsky and Sierra [36,37] used the density matrix renormalization group 
(DMRG) (with n < 400) to achieve significant improvements over the PBCS 
results for the discrete BCS model, in particular in the regime of the crossover, 
which they found to be completely smooth. In general, the DMRG approach 
is applicable to systems that can be divided into two pieces, called block and 
environment, which interact via a preferably rather small number of states. 
One starts with a small block and environment, computes their combined 
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density matrix, then enlarges both and recomputes the density matrix, etc, 
until a large part of the system has been treated. Dukelsky and Sierra chose 
the block and environment to consist, respectively, of all particle or hole states 
relative to the Fermi sea (for a detailed description of the method, see [37]). 
Since the pairing correlations involving coherent superpositions of particle 
and hole states are peaked in a rather small regime of width A around the 
Fermi energy [compare Figs. 12 or 22], the "interaction" between block and 
environment is "localized", so that the DMRG can a priori be expected to 
work rather well for this problem. 

Finally, Dukelsky and Schuck [41] showed that a self-consistent RPA approach, 
which in principle can be extended to finite temperatures, describes the FD 
regime rather well (though not as well as the DMRG). 

To check the quality of the above methods, Braun [91,39] computed E^°''^'^ (for 
s = 0, 1/2) and Ap^ using Richardson's solution (Fig. 20). The exact results 

(i) quantitatively agree, for d — > 0, with the leading —i^jld behavior for 
^cond obtained in the g.c. BCS approach [21,35,23], which in this sense 

is exact in the bulk limit, corrections being of order 

(ii) confirm that the even ground state energy always lies below the odd one 
(this had independently been proven rigoroulsy by Tian and Tang [42]); 

(iii) confirm that a completely smooth [36,37] crossover occurs around the 
scale d ~ A at which the g.c. BCS approach breaks down; 

(iv) show that the PBCS crossover [35] is qualitatively correct, but not quan- 
titatively, being somewhat too abrupt; 

(v) are reproduced remarkably well by the approaches of MFF [33] and BH 
[34]; 

(vi) are fully reproduced by the DMRG of [36,37] with a relative error of 

< 10~^ for n < 400; our figures don't show DMRG curves, since they 
are indistinghuishable from the exact ones and are discussed in detail in 
[36,37]. 

The main conclusion we can draw from these comparisons is that the two 
approaches based on renormahzation group ideas work very well: the DMRG 
is essentially exact for this model, but the band-width rescaling method of 
BH also gives remarkably (though not quite as) good results with rather less 
effort. In contrast, the PBCS approach is rather unreliable in the crossover 
region. To study generalizations of the discrete BCS model, e.g. using state- 
dependent couplings of the form dJ2ij \jblbj, the DMRG would thus be the 
method of choice. 
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Fig. 21. Log- log plot [91] of some of the curves of Fig. 20(a) for the even condensation 
energy Eq°^'^ [in units of A], for A = 0.224; its asymptotic — A^/(2d) behavior for 
d/A — *■ is shown by the dotted line. 

5.3 Qualitative differences between the bulk and the few-electron regimes 



Does the fact that the exact condensation energy E'^ is always negative, even 
for arbitrarily large d/A, mean that the system stays "superconducting" even 
if it is arbitrarily small? The answer is certainly no, since in the fluctuation- 
dominated (FD) regime, the pairing correlations are qualitatively different 
than in the bulk, superconducting regime. In this section we shall try to make 
this statement more precise by analyzing the qualitative differences between 
the two regimes, with regard to the A- and d-dependence of E^°^'^, and the 
behavior of the occupation probabilities v"^. 

Fig. 21(a) shows, on a log-log plot, the d-dependence of the even condensation 

energy EQ°^'^{d). Note that even on the log-log plot, the crossover of the exact 
^cond f^Qj^ ^Yie bulk to the FD regime is completely smooth. According to 
Sierra and Dukelsky [37], the exact result for EQ^^'^i^d) can be fitted very well 
to the form 



Er'^id) = {2d) - 7jo{ln 2)ujbX^ + ^o{Ad/2uj^) \og{2uju/d) , 

where rjQ and 70 are constants of order unity [37]. The first term is extensive 
(oc Vol) and dominates in the bulk limit; its standard heuristic interpretation 
[92] is that roughly A/d levels (those within A of e-p) are strongly affected 
by pairing, with an average energy gain per level of —A/2. The second term, 
which is intensive and dominates in the FD limit, is equal (up to the numerical 
factor r^o) to the result from second-order perturbation theory [37], namely 
[\dY Ya=i X]^=n+i(2£i~2£j)~-^. This subleading term's (^-independence (which 
was anticipated in [45,108]) may be interpreted by arguing that in the FD 
regime, the number of levels that contribute significantly to E^'^'^ is no longer 
of order A/d: instead, fluctuations affect all n ~ 2ijjY,/d unblocked levels 
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Fig. 22. The occupation probabilities Cj of Eq. (89) for d/A = 0, 0.27, 1.09, 2.17 
and 4.34 [91]. In all three figures, the thick solid lines give the d = bulk BCS 
result, whereas circles and stars represent Cj-values evaluated for discrete j's using 
the exact solution and PBCS method, respectively. For d = 0.27A, the PBCS and 
exact results are indistinguishable, and are shown in (a) as a single long-dashed line, 
which is also virtually identical to the bulk curve. For small d, pairing correlations 
are evidently localized within a few A of ep- With increasing d more and more 
weight is shifted away from ep into the tails; compared to the exact results, the 
PBCS method somewhat overemphasizes this delocalization, which is one of the 
reasons why it produces a somewhat too abrupt crossover. 

within ud of ep (this is made more precise below), and each of these levels 
contributes an amount of order —{Xd)^/d (corresponding, in a way, to its 
selfenergy). Finally, the third term contains the small parameter A/cu-d and 
thus represents a very small correction. 

The A- and volume-dependencies of Eq"^'^ in Eq. (88) strikingly illustrate 
the qualitative differences between the bulk and FD regimes: in the bulk 
regime, dominated by the first term, Eq"^'^ is nonperturbative in A (since 
A ~ 2a;De~^/^) and extensive, as expected for a strongly-correlated state; in 
constrast, in the FD regime, dominated by the second term, i?Q°°'^ is pertur- 
bative in A and practically intensive (up to the weak log d dependence of the 
third term). 

Perhaps the most vivid way of illustrating the qualitative difference between 
the bulk and FD regimes is to study properties of the ground state wavefunc- 
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tion. We shall consider here the correlators [35] 



(89) 



which measure the probability that a level can be "both occupied and empty" , 
and vanish identically for states without pairing correlations. For the discrete 
BCS model identically equals (b^jbj) — {b]bj)^ = ~ i^J Eqs. (46) and 
(49)], which measures the fluctuations in the pair occupation number of level 
J, and it vanishes for any blocked single-particle level. Note that (7? also equals 
(bjbj) — (c]_,_Cj_,_)(c]_Cj_); this form, which was used in [35] and corresponds 
to the diagonal terms under the sum in Eq. (52) for Acan, can be interpreted 
as the probability enhancement for finding a pair of electrons instead of two 
uncorrelated electrons in a single-particle level \j, ±) . 

When evaluated using the grand-canonical BCS wavefunction, {Cj)BGS is equal 
to UjVj = jA'^/{e'j + A^) [thick solid hnes in Fig. 22, the same function as 
that plotted in Fig. 12]. The (Q)bcs's thus have a characteristic peak of width 
oc A around Ef, implying that pairing correlations arc "localized around ep in 
energy space" , which may be taken to be the defining property of "BCS-like 
correlations". Moreover, in the bulk regime d <^ A, the (C'j)Bcs are virtu- 
ally identical to (Cj)exact [long-dashed line of Fig. 22(a)], vividly illustrating 
why the grand- canonical BCS approximation is so successful: not performing 
a canonical projection hardly affects the parameters Uj and Vj if d <^ A, but 
tremendously simplifies their calculation. 

As one enters the FD regime d > A, the character of the correlator (C'j)exact 
changes [Fig. 22(b), circles]: weight is shifted into the tails far from e-p at the 
expense of the vicinity of the Fermi energy. Thus pairing correlations become 
delocalized in energy space (as also found in [33,36,37]), so that referring to 
them as mere "fiuctuations" is quite appropriate. In the extreme case o? ^ A, 
the (C'j)exact for all interacting levels are roughly equal. 

Richardson's solution can also be used to calculate, for a given set B of 
blocked levels, the li-dependencc of the canonical order parameter A^^((i) 
of Eq. (52). Schechter has found [125] that it can be fit to the form A^^^{d) = 
A(l + 7B(i/A), where 7^ is a positive numerical constant, and the linear term 
essentially refiects the factor of d in the definition of A^^^. The fact that A^^^ 
is a strictly increasing function of d is in very striking contrast to the behavior 
of the grand-canonical pairing parameters As{d) shown in Fig. 14(a). 
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Fig. 23. Exact even and odd condensation energies E'^°'^'^ [in units of A] for equally 
spaced levels (dashed line), and the ensemble-average (E'™'^'^) for randomly-spaced 
levels (solid line) [39]. The height of the fluctuation bars gives the variances SEg°'^'^. 



5.4 Effect of level statistics 



Smith and Ambegaokar investigated the effect of level statistics on the crossover 
between the bulk and FD regimes [24] . In contrast to the uniform level spacing 
used in previous works, they employed a random spacing of levels, distributed 
according to the gaussian orthogonal ensemble. Using a g.c. mean-field BCS 
approach, they found, interestingly, that randomness enhances pairing corre- 
lations: compared to uniform spacings (u.s.) [21], it (i) on average lowers the 
condensation energy Eg°''^'^ to more negative values, {E^"^'^) < i?™'^'^(u.s.), but 
(ii) these still are parity dependent, (i^g""'^) < {Ey2'^). These results can read- 
ily be understood intuitively: pairing correlations become stronger the higher 
the density of levels around ep, where pair- mixing costs the least energy. When 
determining the amount of pairing correlations for a set of randomly-spaced 
levels, fluctuations that increase the level density near ep are thus weighted 
more than those that decrease it, so that randomness enhances pairing corre- 
lations. 

Although the g.c. mean-field treatment of Smith and Ambegaokar breaks down 
for mean level spacings much larger than A, just as was the case in [21-23], 
their main conclusions (i) and (ii) are robust. Indeed, these were recently con- 
firmed by Sierra et al. [39], who used Richardson's exact solution to calculate 
^cond fQj. ensembles of random levels [Fig. 23]. Moreover, they found that 
the blocking effect responsible for (ii) manifests itself in the fluctuations too, 
which likewise are parity dependent: for example. Fig. 23 shows that both 
the variances 6E^°'^'^ = [((E"^""*^)^) — (£'cond^2ji/2 ^^^^^ ^^^e randomness-induced 
changes in condensation energies \{E^°^'^) — £'™°'^(u.c)| were larger for even 
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than for odd grains. 



5.5 Finite temperature parity effects 

Although finite-temperature studies of the discrete BCS model are not of direct 
relevance for spectroscopic measurements of the BRT-type (a finite T would 
simply smear out the discrete spectra, thereby blurring their most interesting 
features) , they are important in their own right for extending our understand- 
ing of superconductivity in ultrasmall grains. We hence review several recent 
finite-T developments below. 

To begin, let us note that parity effects are of course not restricted to the 
T — limit discussed so far. To be observable [84-88], they only require 
the temperature to be smaller than the free energy difference dj-" ^ A — 
k-BTln[Ncs{T)] between an odd and even grain. Here Ncs{T) is the effective 
number of states available for quasiparticle excitations at temperature T, and 

for d <^ A is given by Nes{T) = \J 87rTA / d? [84]. Below the corresponding 
crossover temperature where 5J^ = 0, determined by /cbT^,* = A/ \n[Ncs{T*j.)] 

and roughly equal to A/ In \j8TTA?/d'^, the single unpaired electron begins 
to matter: it causes a crossover from e-periodicty to 2e-periodicity in the 
I-V characteristics of mesoscopic superconducing SET's [84-88], due to the 
ground state energy difference S1/2 — Sq A. Since T*^. becomes of order A in 
nanoscopic grains with ^ A, parity effects should survive to temperatures 
as high as the (bulk) superconducting transition temperature Tc itself. 

Regrettably, the canonical methods discussed in the preceding sections be- 
come impractical at finite temperatures, since the number of states that need 
to be considered increases rapidly for T > d,A. On the other hand, g.c. 
finite-T methods, some of which wc review below, are, in principle, inher- 
ently unreliable for d > A. This applies in particular to the simplest of these, 
parity-projected mean-field theory [109,112] (Sec. 5.5.1) and certain varia- 
tional generalizations thereof [25,26] (Sec. 5.5.2): they yield the same sharp 
phase transition as function of temperature for finite systems as for bulk sys- 
tems, whereas on general grounds no sharp transitions are possible in finite 
systems. The reason for this problem is that they neglect fluctuations in the 
order parameter, which become very important in the transition region. The 
sharp transition is smoothed out once fluctuations are included. A rather effi- 
cient way of doing this is the so-called static path approximation (Sec. 5.5.3). 
Its use is illustrated in Sec. 5.5.4 for a calculation of the spin susceptibility, 
which shows an interesting parity effect that should be measurable in ensem- 
bles of ultrasmall grains. 
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5.5.1 Parity-projected mean- field theory 

The simplest finite-T approach that is able to keep track of parity effects 
is parity-projected mean-field theory, first used in nuclear physics by Tan- 
abe, Tanabe and Mang [126] , and, indepedently, introduced to the condensed- 
matter community by Janko, Smith and Ambegaokar [109] and Golubev and 
Zaikin [112]. One projects the g.c. partition function exactly onto a subspace 
of Fock space containing only even or odd (p = 0, 1) numbers of particles, 
using the parity- projector Pp. 

Zf = Tr^'Ppe-^^''-^''''^ , Po,i = |[1 ± (-l)"*] • (90) 

One then makes the mean-field replacement bj — > {bj — {bj)p} + {bj)p , neglects 
terms quadratic in the fluctuations represented by { }, and diagonalizes H 
in terms of the Bogoljubov quasiparticlc operators 'jja of Eq. (58). The self- 
consistency condition Ap/2 = ^dJ2j{bj)p, evaluated in a parity-projected g.c. 
ensemble according to Eq. (90), leads to a gap equation of the standard form, 

J = d E ^(l-E^-)' E,^^{e,-i,y + \Ap/,\\ (91) 

\ej\<0Ju ^ ^ ^ ^ 

which is parity-dependent, via the occupation function fpj„ = {'^]fj")ja)p for 
quasiparticles. Since their number parity is restricted to be p, fpj„ differs from 
the usual Fermi function /?^. The condition 2n-\-p— {N)p fixes the chemical 
potential /i to lie exactly half-way between the last filled and first empty levels 
if p = 0, and exactly on the singly-occupied level if p = 1, implying /j, — in 
both cases [by Eq. (38)]. 

von Delft et al. [21] apphed this approach to the discrete BCS model with 
uniformly- spaced levels, and solved Eq. (91) for the parity parameter Ap/2{d, T) 
as function of both level spacing and temperature. Fig. 24 summarizes their 
results. At zero temperature, Ap/2{d,0) corresponds to the spin-dependent 
parity parameters As=p/2 discussed in Sec. 4.4.2 [cf. Fig. 14(a)], and drops to 
zero at a critical level spacing df^^. The Ap/2 — > limit of Eq. (91) defines 
the parity-dependent "critical temperature" Tc^p{d), which can be viewed as 
another measure of how rapidly pair-mixing correlations break down as func- 
tion of level spacing (although ultrasmall grains of course cannot undergo a 
sharp thermodynamic phase transition, which can only occur if n — > oo). In 
both the even and odd cases, the behavior of Tc,p{d) shows direct traces of 
the parity projection: 

In the even case, Tcfl{d) [Fig. 24, curve D] is non-monotonic as function of 
increasing d, initially increasing slightly before dropping to zero very rapidly 
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Fig. 24. d-and T-dependence of the pairing parameter Ap/2{d, T), calculated using 

parity-projected mean-field theory [21]. Curve A gives the bulk gap A(0, T), with 
A(0,0) = A; curves B-E give Ap/2id,T)/ A as a function of d/A and T/A for p = 
(B, D) and p = 1 (C, E). The critical spacings d^^^ = 3.56A and d^jf = 0.890A 
given here differ somewhat from those in Fig. 14(a), because the present mean-field 
approach differs in minor details (via terms that vanish when d ^ 0) from the 
variational approach of Sec. 4.4.2. 

as d ^ dQ*^^. The intuitive reason for the initial increase is that the difference 
between the actual and usual quasiparticle occupation functions is fpj^r — 
fj^ < for an even grain (becoming significant when ~ A), reflecting the 
fact that exciting quasiparticles two at a time is more difficult than one at 
a time. Therefore the quasiparticle-induced weakening of pairing correlations 
with increasing T will set in at slightly higher T if d ~ A. 

In the odd case, the critical level spacing dfj^(T) [Fig. 24, curve E] is non- 
monotonic as a function of increasing T, first increasing to a maximum before 
beginning to decrease toward dfjl^{Tc) = 0. The intuitive reason for this is 
that for < Ai/2 <S 7", d, the odd j — function fpOa{T) becomes somewhat 
smaller than its T = value of 1/2, because with increasing T some of the 
probability for finding a quasiparticle in state j "leaks" from j = to higher 
states with j ^ 0, for which < Eq^ in Eq. (91). Thus, the blocking-of- 
pair-scattering effect of the odd quasiparticle becomes slightly less dramatic 
as T is increased, so that dfj^ increases slightly. 

It should be noted, however, that although the non-monotonicities of Tcfi{d) 
and (if^^^(T) are intuitively plausible within the g.c. framework in which they 
were derived, their physical significance is doubtful, since they fall in the 
regime where d/Ag > 1 and the g.c. approach is unreliable, due to its ne- 
glect of fluctuations. 

5.5.2 Variational eoctensions of BCS theory 

The above-mentioned results of von Delft et al. [21] were reproduced and ex- 
tended to flnite magnetic flelds by Balian, Flocard and Veneroni, using a more 
general g.c. variational BCS approach [25,26]. It is designed to optimize the 
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characteristic function ^{^) = InTrPpC A{^), where Pp is the parity 

projector of Eq. (90) and A{^) = expi^—J^-y^-yQ-r), with being observables 
of interest (e.g. the total spin) and the associated sources (e.g. the magnetic 
field). This approach goes beyond the usual minimization of the free energy 
[101], since it optimizes not only thermodynamic quantities but also equi- 
librium correlation functions, which can be obtained by differentiating e'^^^^ 
with respect to However, its g.c. version also suffers from the drawback 
of yielding abrupt, spurious phase transitions even though the systems are 
finite. Presumably this problem would be cured if an exact projection to fixed 
particle number were incorporated into this approach, but this is technically 
difficult and has not yet been worked out. 

5.5.3 Static path approximation 

For finite systems, in contrast to infinite ones, fluctuations of the order pa- 
rameter about its mean-field value are very important in the critical regime, 
causing the phase transition to be smeared out; conversely, the spurious sharp 
transition found in the g.c. approaches above is a direct consequence of the 
neglect of such fluctuations. A rather successful way of including fluctuations 
is the so-called static path approximation (SPA), pioneered by Miihlschlegel, 
Scalapino and Denton [94] and developed by various nuclear theorists [127- 
143], while recently an exact parity projection has also been incorporated 
[29-32]. A detailed and general discussion, including a complete list of rele- 
vant references, was given very recently by Rossignoli, Canosa and Ring [30]. 
We therefore confine ourselves below to stating the main strategies of the SPA 
and illustrating its capabilities by showing its results [Fig. 25] for the quantity 

^ln = m'J:[C^,-iC.,)x=o]■ (92) 

ij 

Acan is reminiscentP^ of Acan of Eq. (52), and measures the increase in pairing 
correlation energy due to a nonzero coupling strength A. 

One starts by decoupling the quartic interaction Hj-cd of Eq. (37) into a qua- 
dratic form using a Hubbard- Stratonovich transformation with a complex 
auxiliary field A(r) = A+(r) -|- iA~(r), with Matsubara-expansion A± r = 
X)n 6^^™"^^^ the interval r G (0,/3). The parity-projected partition func- 
tion of Eq. (90) then has the following path integral representation (our nota- 
tion is deliberately schematic; see [30] for a precise version): 



The definitions of Acan and Acan differ by terms of order (d/w/))^/^; for example, 
when evaluating both using |BCS) of Eq. (47) and comparing to Age of Eq. (48), 

one finds (Acan)BCS = Age = (Acan)BCS + 0[{d/uJD)^^^]- 
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Fig. 25. Temperature dependence of the pairing correlation energy Acan of Eq. (92) 
[in units of level spacing d], as calculated in [30] using parity-projected mean-field 
BCS theory (dotted lines), SPA (dash-dotted), CSPA (dashed) and exact diagonal- 
ization (solid lines). A system of 10 equally-spaced, doubly-degenerate, levels was 
studied, containing 10 (left panel) or 11 (right panel) electrons. Acan(7') is shown at 
four fixed values of d/A (thus this figure elaborates Fig. 24), namely 0.60, 0.91, 1.7, 
and 15, labeled by a,b,c and d, respectively (calculated using A = 0.55, 0.45, 0.35 
and 0.2). CSPA data are shown only above the CPSA breakdown, which occurs at 
TcSPA < A/4 for the cases considered. The absence of dotted lines for the cases d 
(even) and c,d (odd) means that for these no nontrivial mean-field BCS solution 
exists. Of course, (Acan)exact is nonzero nevertheless. The abrupt BCS transition 
is completely smeared out for the SPA, CSPA and exact results, for which the 
asymptotic decay at T S> A can be shown [30] to be Acan ~ \{X^d/T)^/^. 

h[A] = J2{ej -fi- Ac^/2)c],c,, - ^(6] A + A*b^) + ' ^ " 



JO- 



Xd 



(94) 



The path integral can be treated at several levels of sophistication: 

(i) In the simplest, one uses a "fixed-phase saddle-point approximation" for the 
"static" n = modes and neglects all n ^ modes, i.e. one fixes the phase of 
Aq + IAq = I AqI e"^" by, say, setting 0o = 0, so that / (IAq dAg is replaced by 
/ d|Ao|, and approximates this integral by its saddle-point value. The saddle- 
point condition for maximizing 2p[| Ao|] then yields the gap equation (91), thus 
this approach simply reproduces the parity-projected mean-field approach of 
Sec. 5.5.1, including its sharp phase transition (Fig. 25, dotted lines). 

(ii) The next-best approximation is obtained if one writes / dAg dAQ — 
/o°° |Ao| d|Ao| Jq'^ d0o and performs the phase integral fully. Remarkably, "lib- 
erating" the phase degree of freedom in this way already suffices to smooth 
out the phase transition [141,30], even if the / d|Ao| integral is again replaced 
by its saddle-point value, provided that the latter is found by now maximizing 
I Ao|.Zp[| Ao|] (i.e. including the factor | Ao| from the integration measure). This 
yields a modified gap equation with a nontrivial solution for arbitrarily large 
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T, i.e. no abrupt transition. 



(iii) For finite systems, fluctuations about the saddle become large in criti- 
cal regions. To obtain an improved description of the latter (Fig. 25, dash- 
dotted lines), the static path approximation (SPA) [94,127-133] incorporates 
all static fluctuations exactly, via a (numerical) evaluation of the full integral 
/o°° |Ao| d|Ao| Jo'' d(/)o over all "static paths". 

(iv) In the so-called correlated static path approximation (CSPA) (also called 
SPA+RPA), small-amplitude quantum fluctuations around each static path 
are included too, by performing the remaining / dA^_^Q integrals in the gaus- 
sian approximation [134-140,29-32]. The CSPA yields qualitatively similar 
but quantitatively more reliable results (Fig. 25, dashed lines) than the SPA, 
but breaks down below a temperature Tcspa, below which the fluctuations of 
the A^_^Q modes become large at unstable values of |Ao|, causing the gaussian 
approximation to fail. 

(v) Finally, in the so-called canonical CSPA one projects the partition function 
not only to flxed number parity (as done throughout above) but also to flxed 
particle number, by performing an integration over the chemical potential 
{before performing any of the A^ integrals) [124,131-133]. However, this too is 
usually done only in the gaussian approximation (and would produce negligible 
corrections to the CPSA results for the quantities shown in Fig. 25). 

Comparisons with exact diagonalization results [142] (Fig. 25, solid lines) show 
that in its regime of validity (T > Tcspa), the CSPA produces results that 
are qualitatively completely similar and also quantitatively very close to the 
exact ones, whereas the quantitative agreement is significantly worse if only 
the SPA is used. Since the CSPA is conccptTially simple, well-documented [30] 
and straightforward to implement, it seems to be the method of choice for 
not too low temperatures. A possible alternative is a quantum Monte Carlo 
evaluation of the path integral (93) [144,145], but the numerics is much more 
demanding than for the CPSA, while the convergence at low T is in general 
rather poor, due to the famihar sign problem of Monte Carlo methods. 

The development of canonical finite-T methods that remain quantitatively 
reliable for d > A and arbitrarily small T is one of the open challenges in this 
field. It would be very interesting if progress in this direction could be made 
by exploiting the integrability [117,118] of the model, using Bethe Ansatz 
techniques. For the FD regime, another possibility would be to develop a 
finite-T generalization of the self-consistent RPA approach of Dukelsky and 
Schuck [41]. 
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Fig. 26. Spin susceptibility xo (Xi) of ^ven (odd) superconducting grain as func- 
tion of T, plotted in units of its bulk high-T value XPauli = (a) Bulk limit 
(d/A <SC 1): the even (solid) and odd (short-dashed) curves were calculated using 
parity-projected mean-field theory, the long-dashed curve using standard (unpro- 
jected) BCS theory, (b) Fluctuation-dominated regime (d/A > 1) for x\ (main fig- 
ure) and xo (inset): All curves were calculated using the static path approximation, 
except the short-dashed curves in the main figure, calculated using Richardson's 
exact solution, and the long-dashed curves for the non-interacting case (A = 0). 

B.B.Jf. Re-entrant spin susceptibility 

For grains so small that d ^ A, the spectroscopic transport measurements of 
BRT are not able, in principle, to reliably detect the effect of pairing correla- 
tions, since in this regime these cause only small changes to the eigenspectrum 
of a normal metallic grain, whose spectrum is, however, irregular to begin with. 
In contrast, thermodynamic quantities do have the potential to measurably re- 
veal the existence of pairing correlations for d ^ A. Since very recently parity 
effects for the spin susceptibility have been observed experimentally for an 
ensemble of small, normal metallic grains [146], it is an interesting and exper- 
imentally relevant question to investigate how pairing correlations affect its 
behavior in superconducting grains. 

This question was worked out in detail by Di Lorenzo et al. [44] , whose results 
are summarized in Fig. 26. The spin susceptibility for an isolated grain is 
defined as 



Xv{T) = 



a".Fp(r, H) 
dm 



(95) 



H=0 



where — —k^T In is the free energy of a grain with parity p and 
is the canonical partition function. 

In the bulk limit [Fig. 26(a)], it is well-known [147] that the spin susceptibil- 
ity decreases below its Pauh value XPauU — '^I^B^/d once T drops below the 
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superconducting transition temperature Tc, since the electrons tend to bind 
into Cooper pairs, which are spin singlets and do not contribute to the spin 
susceptibility. Interestingly, however, the spin susceptibility becomes parity- 
dependent as T is lowered below the crossover temperature T*^. mentioned in 
the opening paragraphs of Sec. 5.5: In the even case, xo exponentially drops to 
zero for sufficiently small temperatures, T <C max(A,ci), for reasons that are 
intuitively obvious in the two limits A ^ d (all electrons bound into Cooper 
pairs) and A ^ d (no Cooper pairs, but all levels doubly occupied). In con- 
trast, in the odd case Xi shows a re-entrant behavior, in that it increases as 
/ib^/T for low temperatures, due to a Curie-like contribution from the un- 
paired odd electron. As a result, Xi(^) has a minimum somewhat below T^,*, 
which can be viewed as a "smoking gun" for pairing correlations, since it is 
absent for odd normal grains. For the latter, Xi(T) [long-dashed A = curve 
in Fig. 26(b)] also has the Curie-hke increase at very low T, but lacks the 
initial pairing-induced decrease as T is reduced below T^.. 

Remarkably, Di Lorenzo et al. found that this re-entrance of Xi survives also 
for d > A [Fig. 26(b)]: although pairing correlations survive here only as 
fluctuations, these are evidently sufficiently strong to still significantly reduce 
Xi(T) relative to XPauU [by several percent even for d/A ~ 50 (!)], before 
the Curie-like increase sets in at low T. Di Lorenzo et al. established this 
result by considering the limits T ^ d and T <^ 2d analytically, using a 
static path approximation to capture the crossover numerically, and checking 
the results for T < d using Richardson's exact solution (they considered all 
eigenstates with excitation energy up to a cutoff A ~ AOd, for grains with 
iV < 100 electrons). This check shows that the static path approximation 
somewhat underestimates the amount of pairing correlations (its minima for 
Xi(T) are too shallow), but in general is in good qualitative agreement with 
the exact results, confirming that it is a useful and qualitatively reliable tool 
for describing the crossover regime. 
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6 Nonequilibrium effects 



So far, we have always assumed that the gate voltage has been tuned such 
that the system is close to a Coulomb-blockade degeneracy point, where the 
Coulomb-blockade barrier min[5£'j^Qt] of Eq. (7) is small (of order d) and 
nonequilibrium effects can be neglected. However, since this requires fine- 
tuning V^, it is a rather non-generic situation. In this section, we discuss 
the opposite case in which this Coulomb barrier is so large that transport 
occurs only in the nonequilibrium regime, for which several interesting new 
phenomena have been observed [11,10,12,13,15]. Among these, we focus 

(i) in subsection (6.1) on the observation of clusters of closely-spaced levels 
in normal grains; 

(ii) in subsection (6.2) on the occurrence of prethreshold structures in the 
excitation spectra of even superconducting grains; and 

(iii) in subsection (6.3) on a direct observation of the crossover from a discrete 
to a continuous excitation spectrum at energies of order -Exhouiess- 

All of these can be understood only by going beyond the simple theoretical 
framework presented in Sec. 2.3: (i) and (iii) require consideration of correc- 
tions to the "orthodox" model for treating the Coulomb interaction, which are 
of order d/gdim and hence become important for sufficiently small grains; and 
(ii) requires consideration of cotunneling processes, whereas we had hitherto 
restricted our attention to sequential tunneling. 

The general significance of the results presented below is that they provide 
a very direct confirmation and illustration of the general theoretical picture 
[66.64] (summarized in Sec. 6.3) for the nature of excitations in disordered 
interacting systems that has emerged during the last decade. The development 
of this picture was initiated by studies of semiconductor quantum dots [66], 
but is here beautifully confirmed for metallic grains too. 

6. 1 Clusters of resonances 
6.1.1 Experimental results 

Figs. 27(a,b) display the excitation spectrum, d// dV versus energy, of two 
different norroaZ metallic grains, of radii roughly 2.7 nm and 3.6 nm (if assumed 
hemispherical). The spectra display four clear features: 

(1) The low resonances of the differential conductance are grouped in clusters. 

(2) The distance between nearby clusters is of order the mean level spacing 
d of the noninter acting electrons in the dot. 
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Fig. 27. (a,b) The excitation spectra (at T = 30 mK, = 0) of two ultrasmall Al 
grains with volumes f» 40 nm^ (a) and ^ 100 nm^ (b) [11]. The first resonance is 
isolated while subsequent resonances are clustered in groups. The distance between 
nearby groups of resonances is approximately the single-particle mean level spacing 
d. (c) Differential conductance for a model system of 7 equally spaced levels, occu- 
pied alternately by 4 or 5 electrons in a current-carrying steady state, calculated [11] 
from Eq. (18) by solving the master equation (17) numerically. The tunneling rate FJ^ 
into each level was chosen to be uniform, F^ = Ft^n, and a temperature T = d/100 
and SU = d/5 were used, where SU is the variance of the fluctuations [Eq. (110)] in 
the interaction energy, see text. Dashed line: inelastic relaxation rate F™*^^ of each 
level was chosen to be larger than the tunneling rate into it (F'"''' = 5Ftun); solid 
line: no inelastic processes (F™*^^ = 0). 
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(3) The first cluster contains only a single resonance. 

(4) Higher clusters consist of several resonances spaced much more closely 
than d. 

Agam, Wingreen, Altshuler, Ralph and Tinkham [11] have shown that these 
features are manifestations of the interplay between nonequilibrium effects and 
electron-electron interactions, where the latter must be treated more carefully 
than in the usual orthodox model for charging effects. The general idea is 
simple: if the threshold voltage for the onset of current through the grain is 
large, the transfer of one electron through the grain by sequential tunneling 
can leave behind several possible particle-hole excitations on the grain. If suf- 
ficiently long-lived, each of these excitations will modify the threshold energy 
for the next electron that tunnels into the grain in a slightly different way, 
and thus cause what would have been a single conductance resonance to split 
up into a cluster of subresonances. 

6.1.2 Reaching nonequilibrium via sequential tunneling 

To explain the interpretation of Agam et al. in more detail, wc shall begin our 
analysis within the framework of the orthodox model defined by Eqs. (8) to 
(10) of Sec. 2.3, where the notation used below is introduced and explained. 
The present grains are so small that the influence of superconducting pairing 
may be neglected, since the single-particle mean level spacing, d ~ 1 meV, 
is much larger than the BCS superconducting gap for Al, which is 0.18 meV. 
We thus describe the grain by the normal-state single-particle Hamiltonian of 
Eq. (13) and (14), 

HZ^^ = £^pot(^ex) + E^fcjc^ = E(e^DA^ex,a + S^)\a){a\ , (96) 

p a 

|a) = |K}), £„ = EciV|^,, + E%<^ (97) 

p 

where, for notational simplicity, we have here used the roman label p for the 
combined labels laT) of Eq. (13). H^^^ incorporates the effect of electron- 
electron interactions only through the electrostatic contribution £'pot(-/Vex) [de- 
fined in Eq. (6)], which, as wc shall sec below, is too simplistic to explain the 
clustering properties (1-4). Corrections to the orthodox model will be intro- 
duced [Sec. 6.1.3 below] once its insufficiency has become apparent. 

For definiteness, suppose that the overall ground state of the grain has N 

electrons and that the measured excitation spectra arise via tunneling tran- 
sitions, with the left barrier acting as bottleneck, between N and {N + 1) 
electron states [see Fig. 28]. Denoting their respective ground states by |G)Ar 
and |G)jv+i = Co|G)jv, where cj is the creation operator for the lowest unfilled 



93 




Fig. 28. Schematic depiction (using the conventions of Fig. 4) of the tunnehng 
sequences, starting from the ground state |G)jv, that lead to clusters of resonances 
in nonequilibrium transport. Filled circles depict the electron configuration of the 
intermediate 1-particlc states (with odd), solid arrows the bottleneck tunnel- 
ing transitions into them; levels drawn close together represent time-reversed pairs 
that are degenerate in the absence of a magnetic field, (a) Spp and Spf [Eq. (100)] 
are two possible sequences via level p that leave the grain in an excited final state, 
\pP)n or \pp')n- (b) Sqq^pp or (c) Sqq^ppi [Eq. (104)] are two possibilities, differing 
only in the hole position {p vs. p'), for the next tunneling sequence via level q; their 
thresholds differ by an amount eSV^'^^^^ [Eq. (109)], since the interaction energies 
of two electrons in levels q and p' or q and p differ [Eq. (110)]. The figure is drawn 
such that p = and q = I are the two lowest-lying empty levels of |G)Ar- 

level of |G)jv, the corresponding ground state energy difference is 



where 5E^^^ = Epot{Nex + 1) ~ -^pot(^cx) is the electrostatic energy cost for 
entering the grain. All iV-electron excited states are particle-hole excitations 
w.r.t. \G)n- We shall denote them by \pP)n = c],Cp|G)jv and \qq,pp)N = 
c^Cgcj,Cp|G) AT, etc., and within the orthodox model of Eqs. (96) and (97), their 
excitation energies are 



i^pp )orth — £p £p : i^qq,pp )orth — Eq + £p £p , (99) 

etc. Analogous definitions hold for + 1. Now, sequential tunneling through 
the grain occurs via sequences of transitions of the type [Fig. 28(a)] 



Spp: \G)n ^ |pO)iv+i ^ \pP)n. (100) 

In step (i), an electron tunnels via the left (bottleneck) barrier into an empty 
level p (with Sp > Eq), leaving the grain in an intermediate state |pO)Ar+i, 
which we shall call a "quasiparticle state" (following [64]); and in step (ii). 
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one tunnels out via the right barrier from a filled level p, leaving the grain in 
the excited (if p ^ p) final particle-hole state \pp)n- According to Eq. (26), 
steps (i) and (ii) become possible only once the following inequalities hold, 
respectively (with Vr = Vr — Vd): 



(i) eV^>S^-^'-S^ , (ii) -eV^>S^-S^-^'. (101) 
Adding Eqs. (101.i) + (98) and Eqs. (101.i) + (101.ii) gives, respectively. 



(i)' eV/2^eVj^>SE;,, + e,, (ii)' eV > S^^ - , (102) 

which have obvious physical interpretations: Condition (i)' makes explicit that 
transport is possible only if the bias voltage is large enough that an electron 
leaving the left lead can overcome the Coulomb blockade barrier SEp^^; the 
lowest-threshold conductance resonance evidently occurs for p = 0. Condition 
(ii)' states that an excited final state \pp)n can be produced whenever the 
energy supplied by the bias voltage exceeds the corresponding excitation en- 
ergy — 6q. Now, (i)' imphes that if the Coulomb blockade barrier is large 
{SEp^^ ^ d), i.e. for transport far away from a degeneracy point, the thresh- 
old bias voltage for the current to begin to fiow will be large {eV ^ d); but 
this would automatically guarantee that (ii) is fulfilled, at least for low-lying 
final-state excitations with an energy cost of only a few d. We can thus draw 
the following important conclusion: Since sequential tunneling transport away 
from a degeneracy point requires a large threshold bias voltage, it will always 
generate excited final states once this threshold is reached. 

Whether an excited final state \pp) n affects transport or not depends on the 
ratio between F^^^^ its inelastic decay rate, and F^^j^, the average tunnehng 
rate across the bottleneck barrier r into an individual discrete level near ^p- 
In all previous sections we had assumed that r™f' ^ ^Imi^ i-^- ^^^^ \pp)n will 
decay to |G)Ar long before the next electron tunnels in across the left barrier. 
In this case, the distance between two so-called "quasiparticle conductance 
resonances", associated with the onset of sequential tunneling through two 
different intermediate quasiparticle states |gO)iv+i and \pO)n+i, is determined 
by the bottleneck condition (101. i) [or equivalently from Eq. (30)] to be 



eSV^^ ^ S^,-^' - S^,-^' = e, - e, . (103) 

However, estimates of F^f' (discussed in Sec. 6.3) show that it is quite possible 
for an ultrasmall grain to have have F^^?' < F^^^. In this case, the grain will still 
be in nonequilibrium when the next electron tunnels, so that new tunneling 
sequences of the type [Fig. 28(b)] 
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Sqq,pp- \m>)N kO,pp)jv+i -H \qq,pp)N, (104) 

(iii) eV^>^qoU-^PP^ i^) ' eV^ > ^Lp ' ^ (105) 

become possible once the corresponding conditions (iii) and (iv) are met. These 
differ from Sqq by the presence of additional particle-hole excitations pp. The 
change in eV^ needed to reach threshold (105, iii) from threshold (101, i) is 



[see Fig. 28(b)], and each time this is positive, the conductance should show 
another resonance (if e6V^^^^ < 0, the sequence Sqq^pp simply increases the 
peak height associated with the sequences Spp through level p, without causing 
another peak, i.e. these do contribute to the current, although they are spec- 
troscopically "hidden"). But within the orthodox model, eSV^^^^'' = Sq — Sp [by 
Eqs. (99)], which is the same as c^Vl^ of Eq. (103). Thus, the orthodox model 
predicts that even in the presence of nonequilibrium final states, the spacing of 
conductance peaks will be of order d. It is therefore unable to explain the ob- 
served occurrence in Figs. 27(a,b) of clusters with mean inter-cluster spacing 
d and an intra-cluster spacing much less than d. 



6.1.3 Beyond the orthodox model: interaction fluctuations 

The reason for this failure, identified by Agam et al. [11,13], is that the or- 
thodox model treats the electron-electron interaction too simplistically. Its 
general form is 



-^el-el = 2 51 5Z Uijklc\^c]^,Cka'Cla, (107) 
ijkl aa' 

Um = / dri dr2 C/(ri, r2)^*(ri)V';(r2)^Kri)V'fe(r2) , (108) 

where f/(ri,r2) is the (screened) interaction potential. The orthodox model 
keeps only "diagonal" terms in which two pairs of indices are identical {Uijji, 
Uiijj, Uijij), namely 



4T-el = ^C^ex - XdJ2cl+cl_Cj.Cj+ + X,d ^ 4c]^,Q^.C,> + 0{d/ Qdir.) , 

ij ijjCra' 

and neglects both "off-diagonal terms" Uijki (no equal indices) and ij-dependent 
fluctuations in diagonal terms. (A and are dimensionless, volume-indepen- 
dent interaction constants in the pair and spin channels; we took A = As = 
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in this section, since we are not interested in superconducting or magnetic cor- 
relations here). Now, by studying the correlations of eigenfunctions in chaotic 
systems (and an irregularly shaped ultrasmall grain may be viewed as chaotic), 
the neglected terms can be shown [13,67-71,148] to be small as d/g^i^, where 
Qdim = -E'Thouioss/c^ is the dimeusiouless conduct anc^^ and -Exhouicss the Thou- 
less energy [cf. Eqs. (4) and 5]. Intuitively, one may say that the smaller 
1/fi'dim, the more uniform are the wavefunctions and the smaller the fluctu- 
ations in the interaction energy. Indeed, the neglect of these fluctuations is 
parametrically justified for bulk samples (l/^'dim — 0) and in mesoscopic sam- 
ples (1/5'dim ^ !)• For ultrasmall grains, however, l/(7dim can become of order 
unity, so that interaction fluctuations are expected to become measurable. 

Agam et al. [11,13] attributed the above-mentioned clustering properties (1-4) 
to a combination of precisely such interaction fluctuations and the presence 
of nonequilibrium excitations on the grain: they argued that the energy of 
the added electron in level q of the intermediate state |gO,pp)Ar+i depends, in 
general, on precisely which particle- hole excitation pp is present. In particular, 
two intermediate states that differ only by the interchange of a hole between 
levels p and p' , say |gO, pp) at+i and |gO, pp') will have thresholds that differ, 
according to Eq. (105, iii), by [see Fig. 28(b,c)] 

e6V^^^^ = i^^^^l — £pp) — i^^^p' — £ppi) ■ (109) 

Agam et al. argued that the magnitude of this expression may be estimated 
using the Hartree term of the interaction (107), which gives Ugp'p'g — Ugppg, i.e. 

eSV^piPP')^ J dri| dr2|^,(ri)pt/(ri,r2)(|V^,Kr2)p-|^p-(r2)|2) . (110) 

Let 6U denote the variance of (110), evaluated either by ensemble or energy 
averaging. If diffusive dynamics is assumed, 6U can be shown [68] to be of order 
cd / Qdim^^'V^ where c is a geometry-dependent constant of order unity. Thus, 
each quasiparticle resonance associated with the set of sequential tunneling se- 
quences SqQ^pp through a level q should split up into a cluster of subresonances, 
spaced by d/g^im, associated with different nonequilibrium particle-hole config- 
urations pp. Since in Fig. 27(a,b) there typically are about 5 subresonances 



Strictly speaking, in this context one should use ^dim = hji/d, where 71 is the first 
non- vanishing Perron-Frobenius eigenvalue [67,13], which depends on the geometry 
and the amount of disorder present. Eg. for a pancake-shaped grain with thickness 
z, radius r and diffusive dynamics, Agam [13] gives gdim {hv-p /2rd){z/r)'^ . 
^"^ This estimate does not take into account a change in the electrostatic potential 
due to the extra electron in level q. Blanter, Mirlin and Muzykanskii [69] showed 
that the latter effect will lead to even stronger fluctuations, namely 5U ~ d/y/gdim- 
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per cluster, we can take gdim — 5 to be an "experimental estimate" for the di- 
mensionless conductance, which is entirely reasonable for grains of the present 
size .[3 

Whether the lowest-lying quasiparticle resonance splits into subresonances or 
not depends on the threshold voltage for the onset of transport and on the 
grain's number parity. Let Eq and Ei denote the two lowest-lying empty levels 
in |G)Ar. The lowest-lying resonance involves tunneling sequences Sop via level 
p = 0. A necessary condition for it to split into a cluster is that the level p 
from which the second electron leaves the grain [step (ii) of Sop] should come 
from a Kramers doublet not involving the level Eq into which the first electron 
tunneled [step (i) of 5*0^], implying Eo — Ep > d, which requires a threshold 
[by Eq. (102, ii)'] of eV > — Sq > d. This condition is evidently not met 
in either of Figs. 27(a,b), for which the threshold voltage at the first peak 
is significantly smaller than the inter-cluster spacing, which explains in both 
figures why the first peak not split into a cluster. 

If this condition is met, then the lowest-lying peak will or will not split into 
subresonances, depending on whether Eo and Ei are degenerate or not, respec- 
tively, as can be seen by choosing p = and g = 1: If is even so that Eo = Ei 
(by Kramers degeneracy, see Sec. 3.1), the lowest quasiparticle resonance will 
have a cluster of sub-resonances on its upper shoulder, generated by those tun- 
neling sequences 5*15,0^ through level 1 for which edV^^^^^ > 0. On the other 
hand, if N is odd so that Ei — Eo — d, the threshold difference eSV^^^^^ is of 
order d too, hence the Siq^op cluster will be separated from the lowest-lying 
resonance Sop, which remains unsplit, by an energy of order d. 

In general, when M available states below the highest accessible energy level 
(including spin) are occupied by M' < M electrons, there are (^^^ different 
occupancy configurations. The typical width of a cluster of resonances in this 
case is of order W^^'^d/ g^^i^, where W = min(M — M',M'). The width of a 
cluster of resonances therefore increases with the bias voltage. The distance 
between nearby resonances of the cluster, on the other hand, decreases as 
{W^^'^d/gdim)/y^l^- This behavior can be seen in Fig. 27(c), which shows 
the differential conductance obtained in a model calculation by Agam et al. 
[11] and explicitly demonstrates the splitting of resonances induced by fixed 
fluctuations 6U in the interaction energy. 



° For comparison, note that a similar order of magnitude can obtained from 
Eqs. (2), (4) and (5,ii), for a pancake-shaped grain with Volume ~ Anr'^z and 
a ~ (r/z)'^ [cf. footnotes 3 and 4], for which ^dim ~ 50z^/r(nm)~^: for radius 
r ~ 3 nm and thickness z ~ 0.7 nm this would give gdim ~ 6. It is clear, though, 
that such estimates depend very sensitively on the relative magnitudes of r and z. 
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Fig. 29. The tunneling spectra [10], for five different values of the gate voltage 
V^, of a largish Al grain with (a) an odd and (b) an even number of electrons 
in its 1^ = ground state. The top curve in each case corresponds to at a 
degeneracy point, i.e. tuned such the peaks occur at the lowest possible values for 
the bias voltage As is changed to shift them to larger \V\ values (by amounts 
S\V\ = 0.2, 0.4, 0.6 and 1.2 mV), the peaks broaden and develop substructure. The 
corresponding curves have been artificially shifted in energy to align peaks due to 
the same eigcnstatcs. The sample is the same as in Figs. 2, 7, 9 and 10. All curves 
were taken at T = 50 mK and, to improve spectroscopic resolution, at H = 0, so 
that the Al leads were superconducting [this is the reason for the dip to the right 
of the first peak in the topmost curve, cf. Sec. 2.4.3 and Fig. 6(b)]. 

6.2 Prethreshold structures in superconducting grains 

Nonequilibrium effects also cause interesting anomalies in Al grains large 
enough for superconducting pairing correlations to be important [10]. Figs. 
29(a) and 29(b) show the evolution with gate voltage Vg of the odd-to-even 
and even-to-odd tunneling spectra, respectively, of the same grain as that dis- 
cussed in Sec. 4.1. The top curves in both (a) and (b) correspond to 1^ at a 
degeneracy point, and show features familiar from Sec. 4.1 and indicative of 
pairing correlations: a distinct gap (of 0.55 meV) between the first and second 
main peaks in the odd-to-even case, and its absence in the even-to-odd case 
[cf. Figs. 9 and 10]. As Vg is changed away from the degeneracy point, we note 
the following salient features: 

(a) Odd-to-even spectrum [Fig. 29(a), iV — >• iV + 1 transitions, with N odd]: 
(al) The height of the first peak rapidly decreases and a structure of "prethresh- 
old" subresonances develops on its low-voltage shoulder, i.e. below the 
threshold of the first main peak. 

(a2) The characteristic spacing between these subresonances is comparable to 
that of the first few peaks of the even-to-odd spectrum in Fig. 29(b). 

(b) Even-to-odd spectrum [Fig. 29{h), N -\- 1 ^ N transitions]: 
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(bl) For Vg sufficiently close to the degeneracy point (top three curves), the first 
few peaks of the differential conductance are merely shifted by increasing 
Vg, without changing their shape. 

(b2) For sufficiently large changes of Vg (lowest two curves), the resonances are 
rapidly smeared out by increasing V^. 

Properties (bl) and (b2) are easily understood: due to the large spectral gap of 
0.55 meV for creating an excitation on an even grain, the even-to-odd spectrum 
will be affected by nonequilibrium excitations only once the bias voltage e\V\ 
exceeds this value. Therefore, the low- lying peaks of the even-to-odd spectrum 
are smeared out by Vg only for those curves for which Vg has shifted the bias 
threshold to 5\eV\ > 0.55 meV, namely the lowest two of Fig. 29(b). 

Property (al), however, stands in direct contrast to the predictions of the se- 
quential tunneling mechanism of Sec. 6.1.3, namely that for odd A^, the lowest 
peak of the spectrum does not split into subresonances [cf. third paragraph 
after Eq. (110)]. Moreover, the existence of the prethreshold peaks mentioned 
in (al) suggests that nonequilibrium is reached even while \eV\ still is below 
the threshold for the onset of equilibrium sequential tunneling through the 
bottleneck barrier (say L), namely 



eVL > «^G^' - ■ (111) 

Agam and Aleiner [12] have therefore proposed that in the present case, 
nonequilibrium is established not by sequential tunneling, but instead by co- 
tunneling. A cotunneling process with \G)n as initial state involves a coherent, 
second-order transition of the form 



|G)^ A |/3);v+i ^ \a)N , (112) 

in which an electron tunnels into level p via the bottleneck barrier (say L), 
producing a virtual intermediate state \P)n+i, and one tunnels out of level 
p via the other barrier (R), leaving the grain in the final state \a)N- (Since 
the grain is pair-correlated, |Q;)Ar and \P)n+i are, in general, of the form of 
Eq. (40).) The energy conditions for cotunneling are 



(i) eV^<S^^'-S^, (ii) eV>S^-S^. (113) 

Here (n) is the analogue of Eq. (102, n)' for sequential tunneling, and (i), 
which should be contrasted to Eq. (101, i), is the condition that the |G)Ar — > 
\P)n+i transition is virtual, not real (once the inequality (113,i) is violated, 
cotunneling turns into sequential tunneling). Since the conditions (113) are 
much less stringent than (111), cotunneling sets in long before equilibrium 
sequential tunneling. 
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The current due to pure cotunneling is of course typically very small, since the 
cotunneling rate Tpp oc T^T^ is smaller than that for cither step of a sequential 
tunneling sequence by one factor of F [the latter are defined in Eq. (15)]. 
However, if the excited state \a) ^ produced by cotunneling is sufficiently long- 
lived, it can enable the next electron to enter the grain via sequential tunneling 
whenever the condition 

eVL > S^^' - £^ (114) 

is satisfied, which will cause a step-like increase in current as eV is tuned past 
this "nonequilibrium sequential tunneling threshold". Now, since the latter 
lies below the threshold of Eq. (Ill) for equilibrium sequential tunneling by 
an amount — £q > 0, a set of subresonances will occur on the low-bias 
shoulder of the first main resonance, precisely as observed in (al). Moreover, 
the spacing between these subresonances will rcficct that between consecutive 
excited states £^ of the odd grain, i.e. between the peaks measured in the 
even-to-odd excitation spectrum, in accordance with property (a2). 

Agam and Aleiner [12,13] worked out these ideas quantitatively and concluded 
that the rate of creation of nonequilibrium excitations by cotunneling is indeed 
sufficient to explain the data of Fig. 29(a): they showed that the smallness of 
the cotunneling rates due to their being quadratic in F^ is offset to a large 
extent by the increase in the density of states near e-p of both the grain and 
the leads, due to superconducting pairing correlations in both. 

6.3 Discussion of relaxation rates 

As emphasized above, nonequilibrium can only be established if relaxation 
processes are slower than tunneling processes, F^"*^' < ^tmi for £ — o?, hence 
it is important to understand the mechanism contributing to F™'^'. Agam et 
al. [11] discussed several possibilities: (i) The electron-phonon interaction; (ii) 
the electron-electron interaction beyond the Hartree-Fock approximation; (iii) 
Auger processes in which an electron in the dot relaxes while another one in the 
lead is excited; (iv) the relaxation of an electron in the dot as another electron 
tunnels out to the lead; and (v) thermalization with the leads via tunneling. 
Of these, it is easy to see that (iii) to (v) are at best small corrections: (iii) is 
exponentially small in the ratio w/x the barrier width w to the screening 
length X, and moreover the interaction between electrons on both sides of the 
tunnel junction can take place only within a very limited volume; and (iv) 
and (v) are small, since they clearly happen on time scales larger than the 
tunneling time, (i) and (ii), however, require a careful discussion. 
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6.3.1 Electron-phonon interaction 



Since the temperature of T ~ 30 mK at which the experiments are being 
done is much smaller than the mean level spacing, the probability of phonon 
absorption is negligible, and only emission may take place. The sound velocity 
in aluminum is f sound = 6420 m/sec, therefore the wavelength of a phonon 
associated with relaxation of energy e ^ d = 1 meV is approximately 5 nm, 
the same as the system size. In this regime, Agam et al. estimate the phonon 
emission rate to be 



where e-p is the Fermi energy (11.7 eV in Al), and p is the ion mass density 
(2.7 g/cvo? in Al). This rate is that of a clean metal but reduced by a factor of 
Te\d/h, where Tei is the elastic mean free time [149] which, in ballistic systems, 
equals the traversal time across the system of an electron at the Fermi level. 
Assuming ballistic motion this factor is of order 10~^. The resulting relaxation 
rate for e ~ d is therefore of order F^lcj'^ ^ 10® sec~^ which is similar to the 
tunneling rate Ftun ~ 6 • 10® sec"^ (corresponding to an increase in current 
of 10""^'^ A at the first current step). Thus, by increasing the resistance of the 
tunnel junctions, it should be possible to cross over to the near-equilibrium 
regime shown by the dashed line in Fig. 27(c). 

6.3.2 Electron- electron interaction beyond Hartree-Fock 

The effect of the electron-electron interaction on the resonance widths mea- 
sured in a tunneling experiment, or more generally on the life-time of quasi- 
particle excitations in a finite, disordered system of interacting electrons, has 
been elucidated by Altshuler, Gefen, Kamenev and Levitov (AGKL) [64] . Since 
their conclusions are directly relevant for the interpretation of nonequilibrium 
tunneling spectra, we briefiy review them here. 

Consider again tunnehng transitions from the to (A" -|- 1) particle Hilbert 
spaces, which probe the excitations of the latter. In Sec. 6.1.2, we took these 
excitations to be simply particle-hole excitations with respect to the Fermi 
ground state |G)Ar+i. such as \pp)n+i- \qQ-PP)n+i, etc.; in the nomenclature 
of AGKL, these arc called Hartree-Fock states, |pO)jv+i is called a "quasiparti- 
cle", \qO,pp)N+i two quasiparticles and a quasihole, etc. Now, while these are 
exact eigenstates of the orthodox model Hamiltonian i?5?orth' [Eq- (96)], they 
will decay in the presence of the more general interaction i^ei-ei of Eq. (107), 
e.g., a quasiparticle will decay into two quasiparticles and a quasihole, etc. 

A quasiparticle state with energy e will thus acquire a lifetime l/r™^\ and in 




(115) 
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general should be viewed as a wave packet, with packet width F™^^ (we set 
h— 1 in this section), constructed from the exact eigenstates of the full Hamil- 
tonian. In a finite-sized system with a discrete spectrum, a quasiparticle will 
be well-defined if F™^^ < d (note that this condition is more stringent than the 
one for infinite systems, namely F^"^' <S e). Moreover, the quasiparticle wave 
packet width F^"*^' can be identified with the measured width of a resonance 
in the single-particle tunneling density of states (DOS) probed by an incident 
tunneling electron, since the mathematical description of a single-electron tun- 
neling process involves projecting a single-particle state onto exact eigenstates 
of the system [this is made explicit in Eq. (12)]. In principle, in a finite sys- 
tem every exact eigenstate should produce an infinitely sharp conductance 
resonance. However, in reality only a small fraction of the exact eigenstates 
overlap sufficiently strongly with the one-particle excitations produced by the 
incoming electron to be detected as strong peaks by a measurement of finite 
sensitivity. The nature and spacing of these strong peaks depends on the en- 
ergy £ of the incoming electron; under certain conditions, clarified by AGKL, 
they group into clusters of combined width F™®' that can be interpreted as 
quasiparticle peaks. 

The standard golden rule (GR) approach for estimating F^"*^^ gives T^^^ ~ 
£^/£f for bulk systems [150] and Ff q*^ ^(e/dg'dim)^ for finite systems [66], 
but is valid only when the quasiparticle effectively decays into a continuum 
of final states (otherwise the GR will not give the decay rate, but rather just 
a first-order perturbation correction to the energy of a given eigenstate). For 
sufficiently low e, the accessible final states will not form a continuum; for 
example, the density of final states for the decay of a quasiparticle into two 
quasiparticles and a quasihole, namely /2cP, is much smaller than that of all 
many-body states. To analyze how matters change with decreasing e, AGKL 
considered the decay of a quasiparticle into two quasiparticles and a quasihole, 
which in turn can decay into five quasiparticles and three quasiholes, etc., and 
thus mapped the problem onto one of delocalization in Fock space, in which the 
"distance" between two states is measured, roughly speaking, by the number 
of particle-hole excitations by which they differ. A state that is "localized" 
or "delocalized" in Fock space is associated with a quasiparticle peak of zero 
and finite width, respectively. Assuming the dimensionless conductance to be 
large, g'dim ^ 1, AGKL exploited analogies to Anderson localization to study 
the localization-delocalization crossover quantitatively, and constructed the 
following picture. 

The nature of an excitation with energy e depends on the relation of e to the 
hierarchy E** < E* < -Bxhouiess of energy scales, where 



E** = fl'dim/ In(g'dim) , E* = d^/g^ , -EThouless = 5'dim ■ (116) 
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E** separates the localized from the delocalized regime, and altogether four 
cases be distinguished: 

(i) e < E** (localized phase): Here the true many-body excited states are 
just slightly perturbed Hartree-Fock states, i.e. their overlap with simple 
particle-hole excitations of the Fermi ground state is very close to unity. 
In other words, quasiparticles do not decay significantly, and the single- 
particle DOS will consist of a few 5-function-like resolved peaks, which may 
have weak satellites due to coupling to many-particle states involving more 
particle-hole pairs. If the latter are purposefully created by driving the grain 
into nonequilibrium, the satellites will become strong and each quasiparticle 
peak will appear as a cluster. — Since g^im — 5 for Fig. 27, at least its 
first few quasiparticle excitations fall into this regime; note that AGKL's 
conclusion that these are well-described by simple particle-hole states fully 
justifies using only such states in the analysis in Sees. 6.1.2 and (6.1.3). 

(ii) E** < e < E* (delocalized phase with non-Lorentzian peaks): A given quasi- 
particle state is connected to several others far away from it in Fock space, 
so that the number of satellites in a quasiparticle peak rapidly increases with 
e. This causes the quasiparticle peaks in the DOS to have non-Lorentzian 
shapes, but their widths are still -C F^o'g. 

(iii) E* < e < -Exhouiess (delocalized phase with Lorentzian peaks): A given 
quasiparticle state is connected to so many others that it effectively decays 
into a continuum, so its decay rate is given by the GR result, Ff ^ 
d{e / dgdim)^ ■ Since this is still < d for e < i^^Thouiess, the peaks in the DOS 
can still be well-resolved, and have a Lorentzian shape. 

(iv) -EThouiess < £ (quasiparticles not defined): In this regime F^'og > d, hence the 
quasiparticles are not well-defined, and the DOS is a featureless continuum. 

Experimental results strikingly reminiscent of this picture have recently been 
obtained by Davidovic and Tinkham [15], see Fig. 30. Fig. 30(b), in partic- 
ular, seems to show the progression, with increasing energy, of the measured 
tunneling spectrum through all the crossovers. Estimates for ii^Thouiess, cited 
in the caption of Fig. 30, agree within a factor of 2 with the voltagef^ eV at 
which the measured spectra become continuous. Note that in order to observe 
the entire regime up to and beyond the Thouless energy, all within the first 
step of Coulomb blockade, one needs to have -Exhouiess -C Ec- In this respect 
ultrasmall metallic grains, with their huge charging energies, have an edge 
over disordered semiconducting quantum dots, studies of which had initiated 
the analysis of the role of -^Thouless as the upper scale for the observability of 
level discreteness [66]. 



Davidovic and Tinkham argued that -EThouioss should be compared to the voltage 
eV itself and not to the corresponding (capacitatively corrected) excitation energy, 
because the number of possible nonequilibrium excitations is controlled by eV, cf. 
Eq. (102,ii)'. 
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Fig. 30. Excitation spectra of three different nm-scale Au grains measured by Davi- 
dovic and Tinkham [15], at T = 30 mK and if = 1 T. In order to convert from 
bias voltage to excitation energy, V must be multiplied by e times the capacitative 
conversion factor [cf. Eqs. (30) and (31)] indicated for each graph. Assuming the 
grains to be hemispherical, the grains parameters were estimated to be as follows 
for samples 1, 2 and 3, respectively, : grain radius (estimated from Ci + C2 and 
a capacitance per unit area of 0.05 aF/nm^): r = 4.5, 2.4, and 3 nm; estimated 
single-electron level spacing (from Eq. (2), using r as input): d = 0.65, 4.6, and 
2.1 meV; measured single-electron level spacing d (taken as distance between the 
two lowest energy peaks): d = 1, 7, and 1.2 meV; Thouless energy (estimated as 
E'Thouless = hvp/Gr' [Eq. (5) with a = 3, cf. footnote 3], where r' is the grain radius 
estimated from Eq. (2), using the measured d- values as input): -Exhouiess = 37, 75, 
and 40 meV. 

7 Spin-orbit interaction 



For several of the normal-state grains investigated so far, the magnetic-field 
dependence of the excitation spectra showed marked deviations from the sim- 
ple behavior described in Sec. 3.1: g factors differing from 2 were measured 
[7,14,15] (in Au grains, some were as small as 0.28 [15]), and deviations from 
linear iJ-dependence, in the form of avoided level crossings, were observed 
[7,14]. These features were attributed to the spin-orbit interaction, which is 
the subject of this section. A summary of early work on the spin-orbit in- 
teraction in small metallic grains, mainly with regard to thermodynamical 
properties, may be found in the reviews by Buttet [151], Perenboom, Wyder 
and Meier [19], and Halperin [20]. 
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The information on spin-orbit scattering obtained from the analysis summa- 
rized below is complementary to that obtained from studies of bulk systems. 
There, the quantity of primary interest is the average scattering rate (Fgo) 
for an electron, assumed to be initially in a pure spin-up or spin-down state, 
to be scattered into a continuum of states with opposite spin. This rate can 
be measured in studies of weak localization in disordered metals [152] or of 
tunneling between thin superconducting films in a parallel magnetic field [99], 
or it can be related to the effective g factors [153] measured by electron spin 
resonance. In ultrasmall grains, spin-orbit scattering can be studied at a more 
detailed level via its effects on individual, discrete eigenstates of the grain, 
such as the above-mentioned occurrences of anomalous g factors and avoided 
crossings. Moreover, the effects of mesoscopic fluctuations [55] on phenomena 
governed by spin-orbit scattering can be observed directly. 



7.1 Experimental results for weak spin-orbit interaction 

Fig. 31 shows the results of measurements published by Salinas, Gueron, 
Ralph, Black and Tinkham [7,14] on a grain of estimated radius r ~ 3 nm, 
made from nominally pure Al. It presumably contained an unintended defect 
or impurity, because its behavior was similar to that observed in grains inten- 
tionally doped with 4% of Au impurities. The following salient features can 
be identified: 

(i) This is an even-to-odd excitation spectrum, since all three the lowest peaks 
in Fig. 31(a) Zeeman-split in an applied magnetic field H. 

(ii) Orbital effects are negligible: For each pair of peaks the splitting is symmetric 
about their H = resonance energy, and the average of their resonance 
energies shows hardly any //-dependence (fitting the latter to a linear form, 
g'^^^fiBH, yields g^^^ < 0.1, which is negligible). This was to be expected, 
since the largest value of H used here (7T) is much smaller than the scale 
Horh [Eq. (36)] beyond which orbital effects become significant. The latter 
will thus be neglected below. 

(iii) The effective g^^ factor (for low fields) of each Zeeman-split resonance is 
defined by writing its Zeeman splitting, to linear order in H, as 6Ej — 
gfl^BH [cf. Eq. (119)]. Then Fig. 31(b) gives gf = 1.84 ± 0.03, gf = 
1.68 ± 0.08, and gf = 1.76 ± 0.05. 

(iv) An avoided crossing occurs between resonances If and 2 J, in Figs. 31(a,b); 

(v) the areas under the corresponding peaks 1| and 2 J, in Fig. 31(a), plotted 
in Fig. 31(c), are strongly H -dependent: the smaller one (2|) grows at the 
expense of the larger one (1|) in the avoided-crossing region, while their 
sum remains approximately constant. 
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Fig. 31. The effect of the spin-orbit interaction on the excitation spectrum of an 
Al grain of radius r ~ 3 nm [14] : (a) Differential conductance for a range of apphed 
magnetic fields, at T = 50 mK. The curves are offset in dl/ dV for visibility. 
Orbital state no. 2 gives small but visible resonances at low H. Small changes 
in offset charge (Qq) occurred between the 0.1 and 1 T scans and between the 6 
and 7 T scans, shifting peaJi positions. The 0.1 and 7 T scans have therefore been 
shifted along the voltage axis, to give the best fit to a linear dependence for peak 
1|. Dotted lines are guides to the eye. (b) Magnetic- field dependence of the peak 
positions of (a), which have been converted from voltage to energy by a conversion 
factor of eCi^/{C-^ + Cr) = e(0.53 ± 0.01), R being the bottleneck barrier here [cf. 
Eqs. (30)]. Dots represent data points. Dashed lines are extensions of the low-field 
linear dependence of the energies on H. Solid lines show the result of the spin-orbit 
model of Sec. 7.2, describing the avoided crossing between levels |1|) and |2J,). (c) 
Areas under the first three conductance peaks, which are equal to the current steps 
51 at the corresponding resonances. Note the exchange in weight between 5Ir^ and 
(5/2|, which occurs near the avoided crossing between the eigenstates |1|) and |2 j) in 
(b). (d) Magnetic-field dependence of the coefficients x(i^) and r]{H) in Eq. (125), 
which relate the eigenstates |lt) and |2J,) to the basis states |lT)p and |2J,)p. (e) 
and (f): Markers: Tunneling rates and V]^ for transitions out of the grain's 
eigenstates |lt) and |2|), extracted from the measured current as discussed in text, 
assuming that the relaxation rate of nonequilibrium excitations of the grain is either 
smaller (e) or larger (f) than the tunneling rates. Lines: fits using the predictions of 
the spin-orbit model. 
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Properties (iii) to (v) contrast strongly with the behavior observed for pure 
Al grains [cf. Fig. 10]; the latter typically have (yf-factors close to g''^'^^'^ = 2 [see 
also footnote 12 in Sec. 4.2] that vary only weakly from level to level (if at 
all), their spectra show no clear avoided crossings, and they have essentially 
if-independent peak areas. 

7.2 A simple model for the spin-orbit interaction 

Salinas et al. [14] have attributed anomalies (iii) to (v) to spin-orbit scattering 
due to some defect in the grain, and have constructed a simple model within 
which their measurements could be analyzed quantitatively. Using a somewhat 
more general notation than they did and guided by the reviews of Refs. [19,20], 
we shall consider a Hamiltonian ^ftot ~ Hq + Hg^, to be called the spin- orbit 
model, of the following simple form: 

Ho = Yl ^j<^fi 4<^,0 S>,0 ' = £j,0 + |(7g'P"'">B-H" , (117) 

jo- 

H,o = J2<j4a,0Cj,,0, With = ^J^* = -V], , v]^ = 0. (118) 

ija 

Since time-reversal symmetry guarantees each orbital level to be two-fold de- 
generate in the absence of an applied magnetic field {H — 0), we have adopted 
a single-particle basis of pairs of time-reversed states |jcr)o, labeled by a dis- 
crete orbital index j and a spin index a = (t, i) = (+!^)- The orbital ei- 
genenergies Ejfi are meant to include the effects of all spin-independent forces 
and interactions, and are taken to have a mean spacing of d. A magnetic field 
Zeeman-splits these levels by ^ag^^^'^n^H. In principle, it can also produce 
orbital diamagnetic effects, which could cause £j,o to be if-dependent and, 
via spin-orbit coupling, contribute to reducing the effective g factors [154,55]. 
We shall neglect these here, however, in accordance with property (ii). H^o 
describes the spin-orbit interaction and is taken to couple only states of oppo- 
site spin, cr and a = —a. The conditions specified in Eq. (118) for the matrix 
elements ensure that Hgo is hermitian and invariant under time-reversal 
[which maps c]^ — >• crc]^ and vf^ — > (vfj)*], i.e. that represents a sym- 
plectic matrix, as required on general grounds for the spin-orbit interaction 
[19,20]. 

The average spin-orbit scattering rate in this model, expressed in terms of 

the average (over i,j) of the matrix elements in Hgo, is Fgo = 27r|f ^'^•[^/(M). 
The effects of spin-orbit scattering on an individual eigenstate of iJtot will 
be "weak" or "strong" depending on whether the dimensionless parameter 
Ago = hTgo/d is <^ 1 or 1, respectively. Since the rate Fso is an intensive 
(volume-independent) quantity (implying that v^, must scale as Vol"^''^), it 
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should be possible to reach the regime of weak spin-orbit scattering in grains 
with sufficiently large d. Indeed, the grain of Fig. 31 is an example of this 
case, since the minimum splitting at the avoided crossing, which is a measure 
of l^^l, is much less than the mean level separation, d. 

Since Htot is quadratic, it can in principle be brought into the diagonal form 
-^tot = Hjcj ^ja-c]a(^ja & linear transformation. The time-reversal symmetry 
of ^^tot for H = implies that the exact eigenstates again come in time- 
reversed pairs, say \j<j), degenerate at H = 0. Since this degeneracy is lifted 
by the magnetic field. 



(119) 



the labels (+, — ) = (t, i) identify states with "predominantly spin-up or 
down" character (though no longer eigenstates of o"^). The linear term in 
Eq. (119) defines the effective g factor g'^^ for the pair of states |jcr), which, 
being superpositions of spin-up and spin-down states, will haveP^ gf^ < g^^'^'^. 
For weak spin-orbit scattering (Ago <^ 1), this reduction can be calculated per- 
turbatively: using the lowest-order perturbative approximations (designated 
by the subscript P) to and Eja to evaluate gf^ , one finds [20] 
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(120) 
(121) 
(122) 



(Eq. (122) follows by expanding (121) to linear order in H.) Since the bare 
energies Ej^o and matrix elements vfj are subject to mesoscopic fluctuations, 
gj^ will be too. This is in accord with the observation of j- dependent, reduced 
gj^ values cited in property (iii) above. 



7. 3 Detailed analysis of an avoided crossing 



In order to analyze the avoided crossing in Figs. 31(a,b) quantitatively, Salinas 
et al. truncated the problem to an effective four-state system, associated with 
the four lowest levels in Fig. 31(b). The idea is to treat the effect of all other 

In situations where orbital effects are important, such as semiconductor quantum 
dots, they may reverse this trend and give g^^ > g^^^*^, see [55]. 
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levels on the lowest four using perturbative expressions evaluated in the low- 
field limit, while treating the spin-orbit interaction among the lowest four fully, 
for arbitrary H . To this end, one introduces four "bare" basis states for the 
truncated problem, |li)p, |lT)p; |2i)p, |2t)p, defined by the H = limit of 
Eqs. (120), with "bare" eigenenergies eja,p given by (121), expanded to linear 
order in H. In order not to overcount the interaction among these states, the 
Y^i sums in Eqs. (120) and (121) should now exclude i = 1 and 2 (the tilde 
on the label p indicates this restriction). The effective Hamiltonian for the 
truncated system therefore is p^ 

Htot= X/ X/ ^i'^'P'^io-,? So">p ~'~ X/ ^12'^1o-,p'^2ct,p + h.c. (123) 
i=i,2 o- o- 

pureol i |2 

^j'T,p-^j,f> + y9j +2^gj)l^BJ^, Qj = gj +7 • (124j 

Note that Qj^ is the contribution of the z 7^ 1, 2 levels to the total effective g^^ 
factor of state |jcr)p. Moreover, simply for convenience in the fitting, a linear 
term gY'^^^sH has been allowed too, to model any shift in the average energy 
of the Zeeman-split pairs [the fits turn out to give (7j™~0, cf. property (ii)]. 

It is now straightforward to solve this truncated problem explicitly and to 
find the eigenstates |li), |1T)) |2i), |2t) of -fftot- Fitting the corresponding 
eigenenergies to the lowest four lines of Fig. 31(b), the level repulsion between 
I If), |2|) is reproduced and an excellent fit obtained for |t>j;2| = 73 ± 4 /xeV, 
gf^ = -0.03 ± 0.04, ^f™ = -0.10 ± 0.06, gf = 1.90 ± 0.04, and gf = 
1.74 ± 0.04. This yields a difference of cjj^ — gj^ = 0.06 for j = 1, 2, which is 

consistent with Eq. (124), since 4|t>j;2p/(£2,o — ^1,0)^ = 0.06. This amount is the 
contribution to g^^ of spin-orbit coupling between states |lcr) and \2a), and can 
be used as estimate for 4\vfj\'^/(P = (2/7r)Aso; the result. Ago — 0.09, confirms 
that we are in the regime of weak spin-orbit scattering. Since the amount 
of 0.06 accounts for only 40% (or 20%) of the total reduction g^^^"^ — gf^ of 
0.16 (or 0.32) for orbital state 1 (or 2), the remainder must come from spin- 
orbit coupling to other {i 7^ 1, 2) states, whose contribution is therefore very 
significant. 

The eigenstates of H^^t have the general form 

|1T) = 7|lT)p + r7|2i)p , |2i) = -ry*|lt)p + 7*|2i)p , (125) 
with I7P + Ir^p = 1 (and similarly for jlj), |2t))- The coefficients ^{H) 

When deriving by reexpressing Htot in terms of Cjo-^p-operators, off-diagonal 
terms beyond linear order in v"- were dropped, since they contain further factors of 
the small quantity vlJiejafi — £1^,0)1 with i = 1,2 and j 7^ 1,2. 
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and r]{H) are fully determined by the parameters obtained from the above- 
mentioned fit, and are plotted in Fig. 31(d). Evidently, their relative magni- 
tudes are "interchanged" as H is tuned through the avoided crossing, and the 
character of, e.g., |lt) changes from being predominantly |lT)p for small H 
to predominantly |2|)p for large H. 

This feature is also the reason for property (v). Salinas et al. were able to 
confirm this quantitatively by analyzing the if-dependence of, say, the rates 
= T^f^ and = r^"^ for a tunneling-off transition via barrier L from one 
of the odd eigenstates |l|)Ar+i or |2|)7v+i into the even ground state \G)n. 
According to Eq. (24), these rates have the general formp^ 

r^i = ErLI(G|Q.|iT)P ^ + I^I't^ ^ (126) 

Icr 

= ErLI(G|Q.|2i)p ^ \r^\'t\, + bpf^ , (m) 

Icr 

^ ErlrKGlQtliDpP , ^ Erlil(G|Q,|2i)pp . (128) 
I I 

The if-dependence of F^^ and F^ is therefore contained completely in the fac- 
tors I7P and |?7p, since F^^ and F^ are if-independent (and fully determined 
by if = data). Salinas et al. compared the predictions of the spin-orbit 
model for F^. and [solid lines in Figs. 31(e,f)] with estimates [markers in 
Figs. 31(e,f)] that were extracted from the measured conductance peak areas 
mentioned in property (v). These estimates involved inverting Eq. (18) for the 
current in order to express the rates F^^^ and F^ occurring therein in terms 
of the (measured) current and the probabilities that the grain is in eigen- 
state I a), and determining the Pa by solving a corresponding master equation 
[Eq. (17)]. The analysis is complicated by the fact all accessible states |a) 
must be considered, not only |1|) and |2|), including nonequilibrium excita- 
tionsp^ that can be created when the bias voltage \eV\ is sufficiently large. 
To keep the model tractable nevertheless, some simplifying assumptions were 
made about unknown model parameters, e.g. that the ratios F^/F^ are the 
same for all a (see Ref. [14] for details). The results depend, in particular, on 
whether the relaxation rates F™''^ for nonequilibrium excitations are assumed 
smaller or larger than the average tunneling rate Ttnn- The "slow-relaxation" 
assumption [F™°' ^ Ftun, Fig- 31(e)] yields good agreement between the pre- 
diction (solid lines) of the spin-orbit model and the experimental estimates 

The rightmost expressions in Eqs. (126) and (127) were obtained by neglecting 
terms containing factors of the form p(j<7|c|- |G). These are not completely zero 
(because, by Eq. (120), \jcr)p is not a pure spin-cr state and contains some spin a 
components), but are smaller than those kept by a factor of order vfj / {ejafl — £ia,o)- 

The shifting of single-particle excitation energies due to nonequilibrium excita- 
tions discussed in Sec. 6.1.3 could not be resolved for the present grain, however. 
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Fig. 32. Properties of ensemble- averaged g factors of ultrasmall metallic grains, as 
function of the strength Ago = hVso/d of the spin-orbit interaction, (a) Sketch [55] of 
{{g)) as function of Ago, for a magnetic field applied in a random direction; the limits 
Ago ^ 1 and ^ 1 illustrate Eqs. (129) and (131), respectively, (b) Distributions 
of the principal g factors g\, g2, 53 for weak, intermediate and strong spin-orbit 
coupling [56]. The data points were obtained from numerical simulations of a ran- 
dom-matrix version of the spin-orbit model of Sec. 7.2 (see [56] for details; the 
parameter Ago corresponds to A^ there, up to numerical coefficients). 

(markers). For the "fast-relaxation" assumption [F™^^ » Ftun, Fig- 31(f)] the 
agreement is not so good, though qualitatively still reasonable. The fact that 
the slow-relaxation assumption works better is in accord with estimates [11] 
(see Sec. 6.3) that the energy relaxation rate due to phonons is F™®^ 10*s~^, 
an order of magnitude less than the estimates of Ftun — 10^s~^ for the tunnel- 
ing rates extracted from the detailed analysis. 

7.4 Distributions for the effective g factors 

The statistics of the fluctuations in g'^^ have recently been investigated theoret- 
ically by Matveev, Glazman and Larkin [55], and independently by Brouwer, 
Waintal and Halperin [56], who studied the dependence of g^^ on the direc- 
tion of the applied field. In this section, we summarize those aspects of these 
results that are relevant for ultrasmall metallic grains. 

7.4.- 1 Distribution of gf^ for a random field direction 

Matveev, Glazman and Larkin [55] studied the statistics of g'f , as defined by 
Eq. (119), for a magnetic field applied in a random direction, in the two limits 
that the spin-orbit interaction is weak or strong. Their results are as follows: 

(i) For weak spin-orbit scattering (Ago -C 1), the average (^g)) of the effective 
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factors, averaged either over an ensemble of grains or over different levels 
of a single grain, can be calculated from the perturbative expression (122). 
MGL's calculation corresponds to replacing P by X^odP 1271 and performing 
the sum by assuming the bare energy levels Sj^Q to obey Wigner-Dyson 
statistics in the orthogonal ensemble. This leads to the following result, in 
which the logarithm arises due to level repulsion 



^pure 



1 Ago In - — 

12 Ago. 



(129) 



(ii) For strong spin-orbit scattering (Ago S> 1), Eq. (122) is not an adequate 
starting point. Instead, MGL calculated the linear term in Eq. (119) pertur- 
batively in the magnetic field if, exploiting general properties of the exact 
H = eigenstates |jcr), which were treated using random matrix theory for 
an ensemble of symplectic matrices. They found the distribution P{g) of g'^^ 
factors to be given by 




exp 



Sg' 



V?>T^Y"{{g'f'\ (130) 
(131) 



Here 6 is a geometrical constant of order unity and r is the system size (e.g. ra- 
dius for a spherical or hemispherical grain); for a diffusive grain, /tr = S-Ddifr/'^F 
is the transport mean free path, and for a ballistic grain, 4r should, roughly 
speaking, be replaced by r. Fig. 32(a) schematically illustrates the behavior of 
{{g)) as function of Ago- The parametric dependencies of the two contributions 
to (((7^)) in Eq. (131) can be interpreted intuitively as follows [55]: the linear-in- 
H contribution to the energies ej± of Eq. (119) of the time-reversed states |j±) 
is given by (j ± \M^\j±)^iBH = ±{M^)jfiBH, where = U + g^'^^'s^ is the 
total magnetic moment, and Iz and Sz are the orbital angular momentum and 
spin, respectively (all in units of h). Therefore, by Eq. (119), gf^ = 2{Mz)j. 
Now, the angular momentum of an electron traversing the "closed trajectory" 
corresponding to a discrete quantum level can be estimated as the typical 
(directed) area Atyp covered by its trajectory divided by the period h/d of its 
motion, (Iz) ~ mAtypd/h'^, and its spin as {sz) ~ 1/VN, where N ^ {h/d)Tso 
is the number of spin-flips due to spin-orbit scattering which it undergoes dur- 
ing one period. It follows that (s^)^ ~ 1/Ago, which corresponds to the first 



'^^ Level repulsion suppresses the probability po of two levels being very close by a 
factor poi^ifl — Sj,o) = vr^|ej — £j\/6(P; the logarithmic divergence in the sum 
is cut off above by d, and below at the energy scale |ej^o — ^j,o\ ~ \vfj\ oc dy/Xso, 
because of the additional level repulsion caused by the weak spin-orbit coupling. 
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term of Eq. (131). Moreover, using Atyp r'^^/gdim [cf- Sec. 3.2] and Eqs. (4) 
and (5) for the dimensionless conductance gdim, one finds that (f^)^ ^ hv/"^ 
[or 0(1)] for the diffusive [or baUistic] case, which corresponds to the second 
term of Eq. (131). 

Note that, in contrast to the case of weak spin-orbit scattering considered 
in most of this section, the orbital diamagnetic contribution, though small, 
cannot be neglected for a grain with Ago 3> 1, since the spin contribution to 
{{g))"^ is small as l/Ago- MGL suggested that for such a grain the parametri- 
cally small factor li^/r distinguishing the diffusive from the baUistic case [as 
illustrated in Fig. 32(a)] can be used to judge which of the two cases applies, 
depending on whether its tunneling resonances are measured to have effective 
g factors of gf^ <^ 1 or 1 < gf^ < 2, respectively. 

Note also that if one defines the correlation field (i^orb) at which orbital dia- 
magnetism becomes dominant by equating the orbital level splitting to the 
mean level spacing, 2{lz)jlJ.BHorh ~ d, and substitutes (Iz) ~ niAtypd/tf, one 
arrives at iforb ~ $o/(7'"^typ), which is (up to a factor of n) the criterion used 
in Sec. 3.2 to derive Eq. (36). If, more formally, i^orb is associated with the 
field at which the crossover between the symplectic and unitary ensembles, 
driven by the orbital effects of the magnetic field, is complete [154], the result 
is again Eq. (36). 



7.4-2 Distribution for anisotropic tensor of g factors 

Very recently, Brouwer, Waintal and Halperin (BWH) [56] pointed out that 
since spin-orbit coupling provides a mechanism for the spin to "notice" aniso- 
tropies in the orbital wave functions, the splitting of a Kramers doublet Sja, 
and hence the effective g^^ factor of a given conductance resonance, should 

— * 

measurably depend on the direction of the applied magnetic field H. In general, 
the splitting of the Kramers doublet has the form [155] 



(%+ - s,.y = ^ib'H -GrH^ I^B^glHl + glHl + glHi) . (132) 

Here (Qj) is a 3 x 3 tensor, the g"^ (a = 1,2,3) are its eigenvalues, called 
"principal g factors" of level j, and the Ha are the components of H along the 
level's "principal axes" (defined as the coordinate axes which diagonalize Qj). 
Qj is isotropic (all three ga — '2) only in the absence of spin-orbit scattering. In 
its presence, it becomes anisotropic, is subject to mesoscopic fluctuations, and 
for large Ago decreases as l/Ago- BWH calculated the distribution P{gi, g2, g^) 
with respect to an ensemble of small metallic grains of roughly equal size, 
neglecting orbital effects (i.e. assuming (/^) = 0, as applicable for a diffusive 
grain with Ztr/^ <S !)■ For sufficiently strong spin-orbit scattering (Ago ^ 1), 
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they found it to be of the form 

P{9i,92,9s) ^U\9l- 9me^'''^'^^''^^ ■ (133) 

a<b a 

Here g"^ — ^{gj + g2+ gi) is the average of {ej+ — £j_)^/(/iB|if|)^ over all di- 
rections of H, and its average over the ensemble of grains, {{g"^)), is oc l/Ago [in 
agreement with Eq. (131)]. Importantly, the factor \gl — g'^\ in Eq. (133) sup- 
presses the probability for the Qj tensor to be isotropic (i.e. gaS all equal), and 
hence favors an anisotropic response of a given level Sja to an applied magnetic 
field. Fig. 32(b), which shows the distributions for the gaS for weak, interme- 
diate and strong spin-orbit coupling [56], illustrates both how increasing Ago 
leads to smaller ^f^'s, and the tendency for gi, g2 and g^ to be unequal. By 
changing variables from gi, g2 and g^ (with the convention that gi < g2 < gs) 
to g, ri2 = \gi/g2\ and r23 = \g2/g3\, in terms of which Eq. (133) becomes 

P(9. r,„ r.) oc 'jfj;' ^f^l^' '^'^ fe-''IHft . (134) 

the anisotropy is seen to be rather strong: the typical value for ri2 (or r23) is 
of order 1/3 (or 1/2), so if the direction of H is varied arbitrarily, changes in 
g"^^ by a factor of 8 to 10 can be expected (independently of the actual value 
of Ago, as long as it is » 1). 

Eq. (133) can be used to recover the results of MGL [55] by associating their 
g with [{Qj)zzY^'^ I the gij^-factor for a magnetic field in the 2;-direction (which 
is a random direction with respect to the grain's principal axes). For g thus 
defined, the resulting distribution P{g) is found to be given by Eq. (130). 

7.5 Experimental results for strong spin-orbit interaction 

Experimentally, three reports of grains in the regime Ago ^ 1 have been pub- 
lished to date, all involving Au, whose large atomic number leads to a much 
stronger spin-orbit interaction than for Al: 

(1) Salinas et al. [14] observed values of gj^ in the range 0.5-0.8, as well as the 
occurrence of avoided crossings, in large Al grains doped with 4% of Au 
impurities. 

(2) Davidovic and Tinkham [15] observed g^^ — 0.28 and 0.44 in two Au grains 

whose estimated radii (assuming hemispherical shapes) were 4.5 nm and 
3 nm, respectively [these grains correspond to samples 1 and 3 in Fig. 30]. 

(3) Davidovic and Tinkham [16] observed g^^ values between 0.2 and 0.3 in 
an Au grain whose estimated radius (assuming hemispherical shape) was 
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1.5 nm. 

The grain of case (3) had the further interesting feature that the measured 
level spacing was very much larger than the free-electron estimate, and that 
the spacing between three subsequent levels increased with applied field in the 
large- regime, in a way reminiscent of a spin multiplet. Davidovic and Tink- 
ham have suggested the possibihty that this might perhaps reflect a ground 
state having a total spin s larger than 1/2, favored by the Coulomb interac- 
tion: if, e.g., the spacing between two orbital energy levels at e-p is smaller 
than the Coulomb interaction, say it, between two electrons in the same or- 
bital state, then, in analogy to Hund's first rule in atomic physics, the ground 
state would be a s = 1 spin triplet in which both orbital states are singly occu- 
pied by electrons with parallel spins, because this allows their coordinate wave 
functions to be antisymmetrized, reducing their Coulomb interaction energy. 
This is well-known to occur in cylindrically-shaped semiconductor quantum 
dots [156], and was predicted to be possible in metallic grains too by Brouwer, 
Oreg and Halperin [157]: they calculated the statistical distribution P{s) of 
the ground-state spin s for an ensemble of small normal metallic grains with- 
out spin-orbit coupling, and found an appreciable probability for values other 
than or 1/2 already for interaction strengths well below the Stoner crite- 
rion, i.e. u well below d (see also [158] for similar work for quantum dots). 
For example, already at the quite modest interaction strength of -u/rf ^ 0.4, a 
ground-state spin of s = 1 was found to be more likely than s = 0. The prob- 
ability to find non- minimal s is reduced if spin-orbit coupling is present, but 
may still be appreciable if the Coulomb interaction is strong [159]. More work 
is required to demonstrate definitively that the clustering observed in [16] is 
due to excitations within spin multiplet and not a non-equilibrium effect (of 
the sort described in Sec. 6). 



7.6 Spin-orbit interaction in superconducting grains 

The Au-doped Al grain mentioned in case (1) of the previous section was suf- 
ficiently large that its odd-to-even tunneling spectrum, shown in Fig. 33(a), 
exhibited the sizeable spectroscopic gap typical of superconducting pairing 
correlations [cf. Figs. 9, 10]. Interestingly, this allowed Salinas et al. [14] to 
study how pairing correlations are affected by spin-orbit scattering, by com- 
paring the magnetic-field dependence of this spectrum to that of a pure Al 
grain of similar size, shown in Fig. 33(b) [the same grain as that discussed in 
Sees. 4.1 and 6.2]. We note the following features: 

(i) N onequilihrium broadening: For the Au-doped grain (a), which had no 
gate, the resonance peaks are somewhat broader than in the undoped grain 
(b). This is believed not to be related to the Au impurities, but is instead 
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Fig. 33. (a) Odd-to-even tunneling spectrum [14] for an Al grain containing 4% 
Au impurities, for a sequence of magnetic fields from 0.03 to 9 Tesla in 1 Tesla 
increments, at T=15 mK. (b) Comparison data for a pure Al grain [10], the same 
as in Figs. 2, 7, 9, 10 and 29. The curves are artificially offset for visibility. 

a nonequilibrium effect, since the spectrum of grain (b), which had a gate, 
showed a similar broadening when purposefully tuned into nonequilibrium [cf. 
Fig. 29(a)]. 

(ii) Reduced g^^ factors: In contrast to grain (b), for which g'^^ = 2 ± 0.05, 
for grain (a) the first two peaks move at low H with slopes gf^ /2 + g^^^ — 
0.41 ± 0.03 and -gf/2 + gf^ = -0.27 ± 0.03, suggesting values for g^^ in the 
range 0.5-0.8. This behavior is typical of strong spin-orbit scattering. 

(iii) Spectroscopic gap unchanged: At H = 0, the spectroscopic gaps between 
the first two states are very similar for grains (a) and (b), namely 2ile — 0.26 
and 0.25 meV, respectively. This nicely illustrates that spin-orbit couphng 
docs not disrupt the superconducting pairing of time-reversed states, since it 
does not break time-reversal symmetry. 

(iv) Increased critical field: Whereas for grain (b) the spectroscopic gap be- 
tween the lowest two states has disappeared already at 4 T, for grain (a) it 
decreases much slower with H, due to the reduced g^^ factors, and is non-zero 
even at 9 T. This implies that the critical field for the paramagnetic breakdown 
of superconductivity [cf. Sec. 4.5] is much increased by spin-orbit scattering, 
an effect familiar from thin films in a parallel magnetic field [99] . 

(v) Avoided crossing: At fields above 6 Tesla, the slope of the energy vs. H 
curve of the first peak of grain (a) changes sign (with the energy decreas- 
ing with increasing H at high fields), suggesting an avoided crossing (with a 
minimum gap of 130 /leV) with the higher- lying levels [similar to Fig. 31(a)]. 
Interestingly, this implies that spin-orbit scattering modifies the details of the 
paramagnetic breakdown of pairing correlations: in contrast to a pure grain, 
where this breakdown is expected to occur rather abruptly when the Zeeman 
energy of a spin-1 state |1) (with one broken pair) crosses below that of the 
spin-0 state (with no broken pairs), here the first two levels never really cross. 
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and are not pure spin states to begin with. Therefore the disruption of pairing 
correlations must occur gradually, as the spin content of the particle's ground 
state changes continuously in the avoided crossing region. 
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8 Ferromagnetic grains 



If the leads and/or island of a single-electron transistor are made from an 
itinerant ferromagnetic material such as Co, Fe or Ni, which have different 
densities of states for spin-up (majority-band) or spin-down (minority-band) 
electrons, transport becomes spin-dependent. This leads to a number of inter- 
esting new phenomena [160], such as a tunneling magnetoconductance (the 
conductance of the SET depends on the relative orientation of the magnetic 
moments of leads and island) or Coulomb oscillations as a function of the 
applied magnetic field (which can shift the chemical potential of the island, 
analogous to the effect of Vg, if spin- up and spin-down electrons have dif- 
ferent density of states near ep). Such effects had previously been studied in 
micron-size ferromagnetic islands [160-169] , and also in nm-scale cobalt grains 
[170,171] at or above hehum temperatures. 

This work has recently been extended to smaller size and/or lower tempera- 
tures by Gueron, Deshmukh, Myers and Ralph [18], who studied individual nm 
scale Co grains using tunneling spectroscopy. They found that these differed 
in several interesting ways from non-magnetic grains, showing, in particular, 
hysteretic behavior and a largcr-than-expected number of low-energy exci- 
tations. Since a complete understanding of the observed phenomena and a 
reliable theoretical framework for their analysis is still lacking, we shall here 
just briefly mention the most important experimental features, referring the 
reader to Ref. [18] for details. 



8.1 Experimental results 

Gueron et al. studied Co grains connected to Al leads via AI2O3 tunnel junc- 
tions, without a gate. Grain radii were estimated to vary in the range 0.5-2 nm, 
which, assuming hemispherical grains and a magnetic moment of 1.7/UB/atom 
(the value for bulk Co), would imply magnetic moments of 90-6000 /ib- The 
tunneling spectra showed well-resolved peaks due to tunneling via discrete 
electronic levels within each particle, quahtatively similar to those of Al and 
Au grains, but with a very different magnetic-field dependence. The latter is 
summarized in Fig. 34, which exhibits the following representative features: 

(i) Hysteresis: Fig. 34(a) shows the small- if behavior of the transition ener- 
gies for the lowest three conductance resonances of a Co grain (sample 1). 
Each line shows hysteresis as H is swept from —0.45 T to 0.45 T (thick lines) 
and back (thin lines); it changes smoothly with H, except at two "switching- 
fields" ±ifsw, at which all three transition energies simultaneously show sud- 
den jumps. The qualitative explanation for these features is that the direction 
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H(T) 



Fig. 34. Magnetic-field dependence of the tunneling spectra of several Co grains 
[18]. (a) Hysteresis curves showing the dependence of tunneling energies on H for 
sample 1, at T = 20 mK. (b,c) Lowest-energy transitions calculated using the Ha- 
miltonian -^spin of Eq. (135), for S = 50 and H oriented 45° from the easy axis, for 
the case where S increases during tunneling (b) and decreases (c). The qualitative 
features are independent of the value of S. (d) Voltage threshold for tunneling in 
one sample which exhibits antihysteretic behavior, (e) Tunneling energies over a 
larger range of H, for three different samples. 

of the magnetic moment vector m of the grain changes as H is increased from 
—0.45 T: initially rh is parallel to H, then rotates toward and past the easy 
direction (say z) as the field approaches zero and changes sign, and at Hg^ 

— * 

abruptly changes direction to align itself with H again. We shall call this 
abrupt directional change "magnetization reversal" . 

(ii) Antihysteresis: Fig. 34(d) shows data for a different grain which exhibits 
antihysteresis, i.e. magnetization reversal occurs before H changes sign. This 
can be explained by the influence of a dipolar magnetic field oriented opposite 
to the applied H, produced by a second magnetic grain adjacent to the one 
through which electron tunneling occurs. The reversed field from the second 
particle can shift the hysteresis curve of the first so that its value of ifgw can 
be negative, while the non-monotonic shifts at large positive H (~ 0.2 T) are 
understood as the magnetization reversal of the second particle. 

(iii) Strong spin-up-down asymmetry: Fig. 34(e) shows the large-iJ behavior 
of the Co tunneling spectra. Once H is significantly larger than ifgw (say 
if > 2 T), all measurable transition energies within a given sample have the 
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same sign of slope for 5E[H) vs. H. This contrasts with Al and Au grains, 
which have hnes of both slopes, with comparable conductance amplitudes, and 
with a degeneracy at if = 0. Thus, for Co there is a strong asymmetry between 
the tunneling probability for spin-up or spin-down electrons. 

(iv) g*^^ factors: The effective g^^ factors for individual transitions vary be- 
tween 0.8 and 1.9 and fluctuate quite strongly from level to level. 

(v) Many low-energy excitations: The energy spacing between tunneling peaks 
is somewhat less than 0.2 mcV, much smaller than expected from the mean 
level spacing d for non-interacting electrons: For grains of radii 0.5-2 nm, the 
calculated density of states (including both sp and d bands) in Co, namely 
0.88 eV~^atom~^ [172], imphes d-values between 0.75 and 40 meV. The level 
spacing for individual sp or d bands would be even larger. Note also that the 
peak spacings of the three samples in Fig. 34 are surprisingly similar to each 
other. 



8.2 In search of a model 

At present, a theoretical framework within which all these features can be 
consistently understood is still lacking. The hysteretic behavior of point (i) 
can be understood qualitatively in terms of a simple model Hamiltonian [18]: 

^spin = -^^Vb^ • S - K^'^fXBSl/^SiS + l) , (135) 

It describes a quantum- mechanical spin S (representing m), Zeeman-coupled 
to a magnetic field, with an easy-axis anisotropy in the z direction, where K^^^ 
is an anisotropy energy per unit |m|. The magnitude of S of the spin is assumed 
not to change with H, since the exchange energy U in Co is huge (several eV), 
and U determines the energy cost for changing S (analogous to the role of 
the charging energy Eq for changing N). By diagonalizing ifspin for a grain 
with N or N ±1 electrons (and spin 5" or 5" ± 1/2), the tunnehng spectrum 
S£{H) = £^^^ — 8^1 can be calculated for spin-increasing and -decreasing 
transitions (see [18] for details). The results, shown in Fig. 34(b,c), do exhibit 
hysteresis and jumps reminiscent of the measured ones. In particular, they 
show that the scale of the anisotropy constant K^^^ can be estimated as 
Ijl-qK^^^ £:j"™P, where £'^^^ is the size of the jump in the tunneling energy 
at ifsw This yields HbK^^^ ~ 0.05 meV for the jumps of Fig. 34(a). Thus, for 
ultrasmall Co grains the "characteristic magnetic energy scale" , which governs 
the low-lying excitation spectrum in the meV range (and in that sense plays a 
role analogous to the pairing parameter A for superconducting grains), is the 
anisotropy energy (and not the exchange energy U). 
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This model fails to capture an important detail, however: it predicts that the 
jumps occurring for 68 (H) as \H\ increases through H^^ will always have 
the same sign (±) as the ensuing large- 1 if | slope of SS{H), but the lowest 
line of Fig. 34(a) is a counterexample (jump goes downward, then line moves 
upward) . Gueron et al. proposed that the latter type of behavior could be ex- 
plained by assuming either that K^^^ changes significantly, due to mesoscopic 
fluctuations, if the grain's electron number changes, or that the magnitude 
of S (classically m) does in fact change with H (despite the largeness of f/), 
but no theoretical estimates for the likelihood of these possibilities exist yet. 
Making such estimates would have to involve a somewhat more microscopic 
approach [173], perhaps in the spirit of a Stoner model, in which the degrees of 
freedom are individual electrons in discrete energy levels, not just their total 
spin S. 

Next, note that the spin asymmetry of point (iii) is not surprising for Co, 
since the exchange field breaks the symmetry between spin-up and down. 
However, since Co is not fully spin-polarized (P ~ 30%, [99]), a priori lines of 
both slopes would be expected to occur, even if in unequal numbers and with 
different amplitudes. Recent data (unpublished) do occasionally show a hne 
with a slope quite different in magnitude from that of the majority of lines, 
although the sign of the slope is still the same. The reason why all observed 
transitions shift as a function of magnetic field with the same sign of the 
slope is still a mystery. Its resolution may involve a better understanding of 
whether the states involved in tunneling are primarily sp or primarily d states 
(whose density of states near e-p differ vastly), whether there is a difference 
in the nature of the matrix elements for spin-up or spin-down electrons to 
enter the many-body eigenstates of the grain, or whether spin-orbit coupling 
within the magnetic nanoparticle is sufficiently strong that the states involved 
in tunneling are coherent mixtures of spin states (rather than being purely 
spin-up or spin-down). 

A possible candidate for explaining the small level spacing of point (iv), apart 
from nonequilibrium effects (which cannot be ruled out in ungated devices), 
are "spatially uniformf^ {k = 0) spin-wave modes" : classically speaking, these 
involve fluctuations in the direction of m about its ground state direction (say 
mo), quantum-mechanically speaking, they involve different Sm eigenvalues of 
the operator mo ■ S. The spin-wave excitation energy can be estimated as 
~ ~ 0.1 meV, which is comparable to the observed inter-peak spac- 

ing. Note, though, that spin-waves of different quantum numbers Sm could 
be expected to have s^-dependent slopes for 6£{H) vs. H, whereas no sys- 



The contribution of exchange energy to the lowest-energy non-uniform spin-wave 
modes can be estimated by quantizing the spin-wave dispersion curve of Co within 
the size of a nanoparticle. This gives an energy (300 meV)(a/2r)2 where a is a lattice 
spacing and r is the grain radius [102], or ~ 1 meV for a 2-nm grain. 
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tematic tendencies for the slopes are discernable in Fig. 34(e). This might be 
held against the spin-wave interpretation, but might also be a consequence of 
mesoscopic fluctuations in g^^ factors [cf. point (iv)]. 

8.3 Dynamics of magnetization reversal 

As the above discussion indicates, a number of features of the tunneling spec- 
tra of Co grains still await detailed clarification. This reflects the fact that 
the electronic structure of Co is much more involved than that of Al or Au. 
However, independently of whether all open questions can be resolved or not, 
the fact that magnetic anisotropy causes each resonance energy in the tunnel 
spectrum to shift reproducibly by on the order of 0.1 meV as if is swept about 
the hysteresis loop has a very interesting potential application: it could be used 
as a tool to perform detailed studies of the dynamics of magnetization reversal 
in individual nm-scale grains, since the jumps occuring in the tunneling spec- 
tra upon reversal of the magnetization allow one to monitor precisely when, 
as function of ramped applied field or of waiting time, this magnetization re- 
versal occurs. Thus this method would be complementary to magnetic force 
microscopy [174], Hall magnetometry [175], and SQUID techniques [176,177] 
for studying magnetization reversal. 

Overcoming the energy barrier between two different directions for the mag- 
netization occurs by thermal activation for large temperatures, or quantum 
tunneling for sufficiently small T. For Co grains of the present size (m ^ 90- 
GOOO/xb), the latter case would be an example of macroscopic quantum tun- 
neling (MQT) of the magnetization, which has been studied in great detail 
theoretically [178]. It should be very interesting to use the relevant informa- 
tion in the literature in order to predict the feasibility for observing MQT 
in an individual Co grain in a SET geometry, and to estimate what grain 
sizes would be optimal for this purpose. Wernsdorfcr et al. [176] found devi- 
ations from purely thermally activated behavior when studying an individual 
Co grain with a diameter of 20 nm, using highly sensitive SQUID techniques. 
Since the grains of Gueron et al. are much smaller, it is quite likely that they 
would show such deviations too. It will be important, however, to clarify to 
what extent damping, due to the coupling to the leads, reduces the chances for 
seeing MQT. Conversely, if it turns out that MQT can indeed be observed, it 
might be possible to investigate the effects of damping on MQT in a controlled 
way by purposefully tuning the grain out of equilibrium. 
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9 The Kondo box: a magnetic impurity in an ultrasmall grain 

In this section we briefly discuss what is expected to happen if a normal 
ultrasmaU grain contains a single magnetic impurity This system, to be called 
a "Kondo box", was studied theoretically by Thimm, Kroha and von Delft 
(TKvD) [57]. The Kondo box is another example of the general rule of thumb 
that we repeatedly encountered in previous sections: when the mean level 
spacing d becomes larger than the energy scale characterizing the system's 
correlations, in this case the Kondo temperature Tk, interesting new effects 
occur. TKvD showed that for > Tk, level discreteness strongly affects the 
Kondo resonance, in a parity- dependent way, and predicted that this should 
lead to measurable anomalies in the conductance through the grain. 

For the impurity concentrations of 0.01% to 0.001% that yield a detectable 
Kondo effect in bulk alloys, an ultrasmall grain of typically W — 10^ atoms 
will contain only a single impurity, so that inter-impurity interactions need 
not be considered. TKvD thus considered a single impurity in an ultrasmall 
metallic grain, described by the impurity Anderson model [179,180] with a 
discrete conduction band: 

H = Ho + edJ2 ^da^d^ + ^ E(c]crCda + ^da^ja) + U c\^C^^c\_C^_ (136) 

where Hq = X]j> ^jcj^c^-^. Here a denotes spin and the c^^ create conduc- 
tion electrons in the discrete, delocalized eigenstates \ja) of the "free" system 
(i.e. without impurity). Their energies, measured relative to the chemical po- 
tential II, are taken uniformly spaced for simplicity: sj = jd + eq — fi. 
creates a spin-cr electron in the localized level of the magnetic impurity, which 
has bare energy far below e-p, and a Coulomb energy cost U{-^ 06) for 
being doubly occupied. The hybridization matrix element v between the lo- 
cal impurity level and the conduction band is an overlap integral between a 
localized and a delocalized wave-function, and, due to the normalization of 
the latter, scales as v (x VoP^^^. Thus the effective width of the impurity 
level, r = TTV^/d, is volume-independent, as is the bulk Kondo temperature, 
Tk = \j2VD/n exp(— 7rerf/2r), where D is a high energy band cutoff. To distin- 
guish, within the grand canonical formalism, grains for which the total number 
of {cja and Cda) electrons is even or odd, is chosen either on {ji = Eq) or 
half-way between two {/i — £o + d/2) single-particle levels, respectively [21,22]. 

The system's correlations can be characterized in terms of the spectral density 
Afiai^^) of the impurity Green's function Gdait) = —i9{t){{cda{t), cli^{0)}) , 
which TKvD calculated using the standard noncrossing approximation (NCA) 
[181,182]. Their results are summarized in Fig. 35. 
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Fig. 35. Results for a magnetic impurity in an ultrasmall metallic grain. (a,b) Impu- 
rity spectral function Ad„{uj) for various values of d at T = O.STk; (a) even, (b) odd 
total number of electrons. The individual curves are vertically offset by one unit 
each, (c) Even/odd dependence of the average peak height of the Kondo resonance, 
as function of T/Tk- For an even box (o), Ada was averaged over a range d centered 
on its central sub-peak, for an odd box (•) over a range 2d centered on its central 
two sub-peaks [as indicated by arrows in (a) and (b)]. The inset shows the same 
quantity as function of d/T^- Numerical Tinccrtaintics are smaller than the symbol 
sizes, (d) Anomalous temperature dependence of the weights Wi (circles) and W2 
(diamonds) of the first two conductance peaks of an even grain (o,0) and an odd 
grain (•,♦), for d = 3Tk and d = Tk (inset). 

(i) Bulk limit: For d ^ T, the impurity's spectral density is indistinguishable 
from the bulk case {d — > 0) [lowest lines in Figs. 35(a,b)]: it shows the 
familiar Kondo resonance near cu — 0, with width of order Tk when T <^Tk, 
which arises due to coherent virtual transitions between the impurity level 
and the conduction band. 

(ii) Subresonances: When d is increased well beyond T, however, the Kondo res- 
onance splits up into a set of individual sub-peaks, reflecting the discrete- 
ness of the conduction band [higher lines of Figs. 35(a,b)]. Nevertheless, 
the Kondo resonance retains its main feature, namely significant spectral 
weight within a few Tk around the Fermi energy, up to the largest ratios of 
(i/max(T, Tk) (— 5) that were considered. This implies that the Kondo cor- 
relations induced by the spin-flip transitions between the impurity level and 
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the lowest-lying unoccupied j-levels persist up to remarkably large valuesP^ 
of (i/max(T, Tk). 

(iii) The Kondo correlations do weaken systematically with increasing d, how- 
ever, as can be seen in the inset of Fig. 35(c), which shows the average peak 
height of the Kondo resonance (which quantifies the "strength" of the Kondo 
correlations) as function of d at fixed T: the peak height drops logarithmi- 
cally with increasing d once d becomes larger than about T. Conversely, at 
fixed it drops logarithmically with increasing T once T becomes larger 
than about O.Sd [main part of Fig. 35(c)], thus reproducing the familiar bulk 
behavior. Quahtatively, these features are readily understood in perturba- 
tion theory, where the logarithmic divergence of the spin flip amplitude, 
t{uo) oc Z^jVw is cut off by either T or d, whichever is largest. 

(iv) Parity Effects: For d ^ the even and odd spectral functions Ada in 
Fig. 35(a) and (b) differ strikingly: the former has a single central main peak, 
whereas the latter has two main peaks of roughly equal weight. This can 
be understood as follows: For an even grain, spin-flip transitions lower the 
energy by roughly Tk by binding the d electron and the conduction electrons 
into a Kondo singlet, in which the topmost, singly-occxvgied j level of the 
free Fermi sea carries the dominant weight, hence the single dominant peak 
in Ada- For an odd grain, in contrast, the free Fermi sea's topmost j level 
is doubly occupied, blocking energy-lowering spin-flip transitions. To allow 
the latter to occur, the two electrons in this topmost level are redistributed 
with roughly equal weights between this and the next-higher- lying level, 
causing two main peaks in A^a and reducing the net energy gain from Tk 
by an amount of order d. This energy penalty intrinsically weakens Kondo 
correlations in odd relative to even grains; indeed, the average A^a peak 
heights in Fig. 35(c) are systematically lower in odd than in even grains, 
and more so the larger d and the smaller T. 

(v) Anomalous conductance: The Kondo correlations for Adai^^) affect the prop- 
erties of the conduction electron density of states on the grain, in a way 
which should be detectable by using single-electron-tunneling spectroscopy 
and measuring the differential conductance G{V) through the grain. TKvD 
showed that contributions to G{V) due to the bare conduction electron den- 
sity of states interfere with contributions coming from the Kondo resonance 
in Adaiu^), in a way reminiscent of a Fano resonance [185]. This interference 
causes the conductance to acquire a distinct, anomalous T dependence: the 
weights Wj under the individual peaks (labeled by j) of G(y) become T 
dependent. (In contrast, the weight Wq under an individual peak of the 
bare conductance, say GoiV), is T independent, since the T dependence of 
the peak shapes of Go is determined solely by the derivative of the Fermi 
function.) This anomalous behavior is illustrated in Fig. 35(d), which shows 
the T dependence of the weights Wi and W2 of the 1st and 2nd conduc- 



Actually, remnants of Kondo correlations have been demonstrated to exist even 
in molecular systems [183,184]. 



126 



tancc peaks that occur when V is increased from 0: When T decreases at 
fixed d = Tk- both Wi and W2 decrease, while at fixed d = 3Tk, Wi de- 
creases whereas W2 increases. The fact that the weights can either increase 
or decrease with decreasing T results from the constructive or destructive 
Fano-like interference effects mentioned above. Moreover, at the larger value 
for d, both Wi and W2 develop a parity effect in the strength of their T de- 
pendence. 

(vi) Strength of anomalies: TKvD estimated that the Kondo-induced T depen- 
dence in peak weights in Fig. 35 should be strong enough to be experimen- 
tally detectable. For e.g. an Fe impurity in an even Cu grain of size (3nm)^ 
{d = 30K, Tk ~ lOK), cooling from T = 2Tk to 0.5Tk should change Wi 
by ~ 7%. (TKvD expect Wi to change some more as T is lowered further, 
but their numerics became unreliable in this regime.) 

To conclude this section, we note that since the effects discussed above result, 
above all, from the discrete density of states near ep, they should be generic 
for ultrasmall grains. In other words, they should be robust against including 
randomness in the model, like j-dependent level spacings dj and hybridization 
matrix elements Vj. 
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10 Summciry and outlook 



The technique of single-electron-tunnehng spectroscopy on ultrasmall metaUic 
grains has proved to be a very fruitful way of probing electron correlations in 
metals, and the way in which these are modified by mesoscopic fluctuations 
and level discreteness. Let us briefly summarize the main conclusions reached 
in each of the preceding sections: 

Sec. 3: For normal grains, the effect of an applied magnetic held is simply to hft 
the Kramers degeneracy of time-reversed states. This can be used to determine 
the grain's number parity. In ultrasmall grains, the effect of a magnetic fleld 
on orbital motion is negligible. 

Sec. 4: For largish Al grains, the observation of a distinct spectral gap in even 
grains and its absence in odd grains is clear evidence for the presence of su- 
perconducting pairing correlations. These can be satisfactorily described using 
the simple discrete BCS model introduced in Sec. 4.2. The blocking of some 
levels by unpaired electrons leads to various measurable parity effects; among 
these, a pairbreaking-energy parity effect should be observable in experiments 
of the present kind, provided the grain size can be better controlled. The 
dominant mechanism by which a magnetic fleld destroys pairing correlations 
in ultrasmall grains is Pauli paramagnetism. Decreasing the grain size softens 
the flrst-order transition observed for thin fllms in a parallel fleld, by reducing 
the number of spins flipped from being macroscopically large for d ^ A to 
being of order one for d ~ A. The grand-canonical variational BCS approach 
fails for > A; nevertheless, it yields a useful framework for a qualitative 
analysis of the experiments, which had d < A. 

Sec. 5: The crossover of the behavior of superconducting pairing correlations 
from the bulk limit (d <^ A) to the fluctuation-dominated regime {d ^ A) 
is parity dependent and completely smooth. This remains true for systems 
with non-uniform rather than uniform level spacings. - Very remarkably, the 
discrete BCS model has an exact solution, due to Richardson, with which 
T = properties can be calculated rather easily. Finite-temperature properties 
for flnite-sized systems can be calculated quite reliably with the correlated 
static-path approximation (provided T > Tcspa)- Thus, the only remaining 
open problem is the development of canonical flnite-T methods that remain 
quantitatively reliable for d > A. - The spin susceptibility x(T) of an odd grain 
shows an interesting re-entrant behavior even for d ^ A, which might be a 
way to detect remnants of pairing correlations in the fluctuation-dominated 
regime. 

Sec. 6: The excitation spectra of nonequilibrium grains show clusters of reso- 
nances. The spacing between clusters is comparable to the mean level spacing. 
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but the spacings between subresonances of the same cluster are much smaller 
than d. Each chister of resonances can be identified with one excited single- 
electron state of the metal grain. Subresonances refiect energy shifts (of or- 
der d/gdim) due to different nonequilibrium occupancy configurations of other 
single-electron states. These energy shifts are manifestations of mesoscopic 
fluctuations of the matrix elements of the electron-electron interaction. - Su- 
perconducting grains that are tuned out of equilibrium by a gate voltage show 
subgap structures in odd-to-even but not in even-to-odd tunneling spectra; 
these features can be attributed to nonequilibrium excitations generated by 
cotunneling processes. - As a function of increasing excitation energy, the na- 
ture of the excitation spectrum has been observed to change from consisting of 
discrete peaks to being continuous, the crossover scale being e ~ i?Thouicss; the 
observed details of this crossover are in accord with theoretical expectations 
for the behavior of quasiparticle lifetimes in disordered interacting systems. 

Sec. 7: The spin- orbit interaction leads to reduced g^^ factors for the split- 
ting of Kramers doublets in a magnetic field, and to avoided crossings. - For 
grains with superconducting pairing correlations, the critical field for which 
these are expected to break down is increased significantly relative to pure 
samples; the H — spectroscopic gap remains unchanged, however, as ex- 
pected due to time-reversal symmetry. - The gj^ factors exhibit mesoscopic 
fiuctations, which can be studied theoretically using random matrix theory. 
Interestingly, the Zeeman-response to a magnetic field has been predicted to 
be anisotropic in the presence of spin-orbit scattering, implying that upon 
changing the direction of the apphed field, the measured gf^ factors should 
vary strongly (by a factor of up to 8-10). 

Sec. 8: Ferromagnetic Co grains have hysteretic tunneling spectra. These can 
be used as a tool to measure the switching field at which the grain's magne- 
tization undergoes a sudden change in direction as an applied magnetic field 
is ramped. There is a strong asymmetry for the tunneling probability of spin- 
up versus spin-down electrons. The spacing between low-lying excitations is 
much smaller than expected from the estimated free-electron mean level spac- 
ing d (and remarkably uniform from one grain to the next). However, it is 
comparable to the grains' anisotropy energy (per unit magnetization), which 
for ultrasmall Co grains is the "characteristic magnetic energy scale" govern- 
ing the low-lying excitation spectrum in the meV range. Many features of 
these measurements are not yet understood in detail, and further work, both 
experimental and theoretical, will be needed to clarify them. 

Sec. 9: For an ultrasmall grain containing a magnetic Kondo impurity, the 
Kondo resonance is strongly affected, in a parity-dependent way, when the 
mean level spacing in the grain becomes larger than the Kondo temperature. 
The tunneling spectrum is predicted to show Kondo-induced Fano-type res- 
onances of measurable size, with an anomalous dependence on temperature 
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and level spacing. 

Prospects for future work: 

Experiment: In the current generation of experiments, the grain's actual size 
and shape cannot be determined very accurately. It would be a great advance 
if fabrication techniques could be developed to the point that grains can be 
used which have been custom-made, by chemical techniques, to have well- 
defined sizes and shapes (e.g. spherical). This would significantly reduce the 
uncertainties which one presently encounters when estimating characteristic 
parameters of the grain, such as the single-particle mean level spacing (i, the 
dimensionless conductance (^dim or the Thouless energy. Moreover, it would 
allow systematic studies of the dependence of various quantities on grain size 
or mean level spacing [for example, it would be interesting to try to do this for 
the pairbreaking energies ilg, of Fig. 19(b)]. Encouragingly, the feasibility 
of using chemically-prepared grains in SETs has already been demonstrated 
several times [114,186,187], though the resulting devices have not yet been 
used for single-electron-tunneling spectroscopy. 

Moreover, most of the devices that have been studied so far had no gate, 
because it is technically very challenging to add one. It would, however, be 
very useful if ways could be found for routinely fabricating gates, because that 
allows both the grain's number parity and the amount of nonequilibrium to 
be adjusted controUably. 

Theory: The behavior of superconducting pairing correlations in an individual 
ultrasmall grain can now be regarded as a subject that is well understood. 
It would be interesting to try to use the insights that have been gained for 
a single grain in order to now study systems of several coupled grains: what, 
for example, is the fate of the Josephson effect between two coupled grains as 
their sizes are reduced to the point that d ~ A? 

The experiments on nonequilibrium effects and spin-orbit interactions revealed 
the importance of mesoscopic fluctuations in remarkably direct ways, and are 
well understood in terms of theories that treat these fluctuations via random 
matrix theory. The success of the latter approach stems from the fact that 
for chaotic systems such as irregularly-shaped grains, "details do not matter" 
much, so that a description in terms of random matrices is appropriate. 

At present, the most intruiging open theoretical problem seems to be that 
of flnding a suitable framework within which to analyze and interpret the 

experiments on ferromagnetic Co grains. The observed phenomena are rather 
rich and complex, and it is currently not at all clear which details "matter" 
and which do not. More experimental studies on this system will be needed 
as guidelines for theoretical attempts to construct a useful model for it. 
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A further stimulus for theoretical work on ferromagnetic grains might come 
from the prospect of exploiting the hysteretic properties of their tunneling 
spectra to study the dynamics of magnetization reversal: it would be inter- 
esting to estimate the feasibility for seeing macroscopic quantum tunneling in 
this system. In particular, one should try to determine whether (or to what 
extent) this phenomenon can survive when a current flows through the grain. 
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A Superconducting leads 



For the case of superconducting leads in the absence of an external magnetic 
field, the analysis of Sections 2.3.2 and 2.3.3 has to be modified as follows (for 

generality, we continue to allow for the presence of a gate; if none is present, 
set Cg = below): Firstly, the lead Hamiltonian (8) must be replaced by the 
BCS Hamiltonian for Bogoljubov quasiparticles. 



k<j 



Iklr 
Vkr 



^kr^klr ~l~ '^krC_j^'^j. , 



Ekr/ Ekr) 



,1/2 



(A.l) 

(A.2) 
(A.3) 



and i?tun of (11) should be rewritten in terms of the 7's. Secondly, the golden 
rule expressions for the tunneling rates now are [instead of (20) to (23)]: 



2n 



^aa' — 2^ 



r* 1 2 



f{E}^r)\'>J'kr\'^5{,£-a — ^a' — Ekr — eVr) 



(A.4) 



ka 



+ (1 - f{E_kr))\VkrfS{Sa - Ea' + E_kr - eVr)] (A.5) 
f{£^ - - eVr) K,^^{Sa - <^a' - 6^) F^^ , (A.6) 
1 - f{Ekr))\Ukr\^5{£o, - E^, + Ekr - 6^) (A.7) 



27r 

yr- _ _ IT--?- |2 
k(T 



+ f{E^kr)\Vkr\^^{Ea — Ea' — E^kr — cK) 
= f{Ea - Ea' + eVr)] Mr,^^{E^ - E^' + 6^) F^^, , 



(A.8) 
(A.9) 



where J^r,qp{^) is the BCS quasiparticle density of states in lead r [see Eq. (33)]. 
To arrive at (A.6) and (A.9), we took the limit T < [so f{Ekr) ^ 0] and 
used [instead of Eq. (15)] 



y E^^Mv } - Ekr) = v\^KM 

k I \'^kr\ 

(the £- antisymmetric part of Iw^rP and \ukr\^ does not contribute). 
Eqs. (A.6) and (A.9) are the results used in Eq. (32) of Sec. 2.4.3. 



(A.IO) 
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B Richardson's exact solution of discrete BCS model 

This appendix summarizes some of the main features of Richardson's exact 
solution of the discrete BCS model. We begin in App. B.l by deriving in 
pedagogical detail, following [40], some of the key results of Richardson's exact 
solution of the discrete BCS model of Sec. 4.2. In App. B.2 we transform 
the eigenvalue equation derived below [Eq. (B.17)] into a less singular form 
more convenient for numerical solution, and in App. B.3 we present analytic 
formulas needed to calculate the correlators {h]hj) and exactly. 

B. 1 Derivation of eigenstates and eigenvalues 
B.1.1 The hard-core boson problem 

For presented purposes it suffices to consider the pair HamiltonianP^ [cf 
Eq. (44)] 

u 

Hu = J2 C^^j^ij ~ 9) bjbj , (where g = Xd), (B.l) 

where j runs over a set U of Nu of non-singly-occupied single-particle levels 
U (degenerate levels with = Sj for i ^ j are allowed), on which the pair- 
creation operators fe] = c]_|.c]_ satisfy "hard-core boson" relations, 

&f = 0, [h^,h],] = 5,,,{l-2h%l [h%h],]=5,yh]. (B.2) 

We seek eigenstates that contain n pairs, satisfy {Hu — £n)\'^n) and are of the 
general form 

U n 

\^n)= E ^(ji,...,Jn)n^llVac). (B.3) 

In his original publications [46-48] , Richardson derived a Schrodinger equation 
for iplJi, . . . ,jn) and showed that its exact solution was simply a generalization 
of the form that ip{ji, . . . , jn) would have had if the 6's had been true (not hard- 
core) bosons. With the benefit of hindsight, we shall here follow an alternative, 
somewhat shorter root, also due to Richardson [119]: we ffist consider the 

In this appendix the shorthand g = Xd is useful, since Richardson's solution 
depends only on this combination. 
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related but much simpler case of true bosons and write down the generic form 
of its eigenstates; we then clarify why this form fails to produce cigenstates 
of the hard-core boson Hamiltonian; and having identified the reason for the 
failure, we show that (remarkably) only a slight generalization is needed to 
repair it and to obtain the sought-after hard-core-boson eigenstates. 



B.1.2 True bosons 

Let bj denote a set of true bosons (i.e. = Sjj'), governed by a Ha- 

miltonian Hij of precisely the form (44), with bj bj. This problem, being 
quadratic, can be solved straightforwardly by any number of methods. The 
solution is as follows: Hu can be written as 

Hu = Y, EjB\Bj + const. (B.4) 
J 

where the new bosons B\ (with normalization constants Cj) are given by 



and the boson eigenenergies Ej arc the roots of the eigenvalue equation 



This is an equation of order Nu in Ej. It thus has Nu roots, so that the label 
J runs from 1 to Nu. As the coupling g is turned to 0, each Ej smoothly 
evolves to one of the bare eigenenergies 2ej. A general n-boson eigenstate of 
Hu and its eigenenergy thus have the form 



lb 11/ 

I^n) = n^lj0)' £n^Y.Ej., (B.7) 



where the n indices Ji, . . . , J„ that characterize this state need not all be 
distinct, since the Bj are true bosons. 



B.1.3 Complications arising for hard-core bosons 

Let us now return to the hard-core boson Hamiltonian Hu- Its eigenstates 
will obviously not be identical to the true-boson eigenstates just discussed. 
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since matters are changed considerably by the hard-core properties of bj. To 
find out exactly what changes they produce, it is very instructive to take an 
Ansatz for similar to (B.7) (but suppressing the normalization constants 
and taking all boson indices to be distinct), namely 



u 



|*n) = n^i|0)' with Bl^Yl 



(B.8) 



and to check explicitly whether or not it could be an eigenstate of Hjj, i.e. 
to check under what conditions {Hu — would equal zero, where Sn = 

El,. To this end, we commute Hu to the right past all the Bl, operators in 
l^n), using 



Hu. n Bl 



n I / u—1 



v=l I \ri=l 



\.IJ,=v+l 



To evaluate the commutators appearing here, we write Hu as 



(B.9) 



u 



Hu^Yl '^^S^j - gBlBo , where B^^Yl ^1 ' 



(B.IO) 



and use the following relations: 



[%,Bl 



2£j — El, 



u 1 _ O/itfc. 



fB.Ul 



[Hu,Bl] = E,,Bl + Bl 



j 2£j - El, 



(B.12) 



Inserting these into (B.9) and using Hu\0) = and Sn = Y^H, Ei,, we find 



/^C/|*n)=<^n|*n) + E 



1 



U 

E 



i?j=i 



u=l 



^ ^gBl b% 

. 2s j - El, 



E„ 



Bl\ n Bl]\0) 

^r,= l(/^) 



n Bl 



|0) 



(B.13) 



Now, suppose we do the same calculation for true instead of hard-core bosons 
(i.e. run through the same steps, but place a ~ on Hu, bj, Ey and E^). Then 
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the second hne of (B.13) would be absent (because the b'jbj terms in the 
second of Eqs. (B.2) and (B.ll) and in (B.12) would be absent); and the first 
line of (B.13) would imply that {Hu — £^„)|\1'„) = provided that the term 
in square brackets vanishes, which is nothing but the condition that the 
satisfy the the true-boson eigenvalue equation of (B.6)! In other words, we 
have just verified explicitly that all true-boson states of the form (B.7) are 
indeed eigenstates of Hu, provided that the Ei, satisfy (B.6). Moreover, we 
have identified the term in second line of (B.13) as the extra complication 
that arises for hard-core bosons. 



B.1.4 The cure: a generalized eigenvalue equation 

Fortunately, this extra complication is tractable: first, we note that 



u 

E 



u 

E 

j 



bl 



2£j — Ei^ Isj 



E,^ 



(B.14) 
(B.15) 



The second line follows via a partial fraction expansion, and remarkably, con- 
tains only operators and no more bjbjS. This enables us to eliminate the 
bjbjS from the second line of (B.13), by rewriting it as follows (we commute 
its term in square brackets to the right, using a relation similar to (B.9), but 
with the commutator (B.14) instead of [Hu^Bj^]): 



n ( / u-1 \ n ( / 
i/=l l\??=l / lJ,=u+lWr]'=u+l 



t Bl - Br 



Eu — Efj_ 



n Bi, 



■|o) 



E 



E 



2g 



"[ El, — Efj 



Bi\ n ^j|io) 



E 

i/=i 



E 



^g 



/i=i/+i Bu E^ 



Bi\ n ^;)io) 



E 



E 



2g 



Bi[ n ^i)io) 



(B.16) 



(The last line follows by renaming the dummy indices v ^ ji in the second 
line.) Substituting (B.16) for the second line of (B.13), we conclude that {Hu — 
^n)\^n) will equal zero provided that Eq. (84) holds, namely: 
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u „ "2^ 



E 



/^o P + E -FT-^-FT^O, for u = l,...,n. (B.17) 



This consitutes a set of n coupled equations for the n parameters Ei, . . . , 
En, which may be thought of as self-consistently-determined pair energies. 
Eq. (B.17) can be regarded as a generahzation of the true-boson eigenvalue 
equation (B.6). It is truly remarkable that the exact eigenstates of a compli- 
cated many-body problem can be constructed by such a simple generalization 
of the solution of a quadratic (i.e. non-interacting) true-boson Hamiltonian! 



B.2 Solving the eigenvalue equation numerically 



In this section we discuss some algebraic transformations invented by Richard- 
son to render Eqs. (84) or (B.17) less singular and thereby simphfy their nu- 
merical solution. 

First note that the solutions of Eqs. (84) can be given an electrostatic inter- 
pretation [116,54], since they are the extrema w.r.t. E^, of the function 

W[{E^}] = E ^ + ^ E Elog |2£, E ^og\E, - E^\, (B.18) 

1/^9 ^ u j v<ix 



where E^ is regarded as a complex number, say E^.^ + iEi,y^ and \Ey\ — 
\jEl^ + Ely, etc. iy[{£^i/}] can thus be interpreted as the potential energy 
of n free unit charges at the locations (Eyx^^vy) in the xy-plane (actually 
lines of charge perpendicular to the plane), under the combined influence of 
a uniform external field —{\l2g)x parallel to the x-axis, and the field of a 
set of fixed charges of strength —1/2, located at the points 2£j on the x-axis. 
The task at hand is to find the (unstable) equilibrium positions of n free unit 
charges. 

We shall consider only the case that all the Sj are non-degenerate (degeneracies 
are considered in [116]). Then every solution of Eqs. (84), i.e. every set {E^), 

coincides at (7 = with a set of n bare energies, say {2£j^} (by inspection; or 
in the electrostatic analogy, the field —l/2g is so strong that the only stable 
configurations are those where each positive unit charges sits infinitely close to 
some fixed —1/2 "host" charge). Richardson showed [51] that as is turned on, 
all E'l.'s initially smoothly evolve toward lower values (as the field weakens, the 
free charges leave their hosts, pulled in the negative x-direction by the field); 
but as g increases further, singular points are reached at which two roots 
approach both each other (and a bare energy), say E^ = E^{= 2ej), and then 
turn complex, becoming a complex conjugate pair, so that £n remains real 
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(as the field weakens some more, the repulsion between the positive charges 
begins to become important; it tends to push some of them off the ,x-axis into 
the xy plane, but does so without breaking reflection symmetry about the 
X-axis) . 

Due to this complication, it is convenient to parametrize the roots in terms of 
purely real variables, proceeding as follows: Denote roots destined to become 
conjugates by {E2a-i,E2a) [with g = values (£i2a-i> ^i2a)' say], with a = 
1, . . . n/2 for even n, with one further purely real root, say £'0, for odd n. Let 
us write ^sa-i = Ca-iVa, = Ca + iVa, where and r)l are purely real, with 
1]'^ > for g > Qa, the critical value where this pair of roots becomes complex. 
Rewriting Richardson's eigenvalue equation (84) in terms of these purely real 
variables by adding and subtracting the equations for i/ = 2a — 1 and u = 2a, 
one readily finds, for each a, 



1_ 2ej - 2Xa{l + Vg) 



{Eo - + ril (41 + Vl + Vl? - Hril 

(1 - Vl) , V- Va 

-vt + v a 



where we introduced the further real variables 



Xa = U- ^2a-l - e2a , Va = vl/ i^l - pi) , (B.21) 

Pa — ^2a — ^2a-l , Cba = Cb — Ca , (B.22) 

and factored out an overall factor of (x^ — Pa) from Eq. (B.20). Moreover, the 
terms involving Eq occur only if n is odd and should be omitted if n is even. 

Compared to Eqs. (84), Eqs. (B.19) and (B.20) have the advantage that their 
roots are always real and that no singularities occur any more. They can 
thus be solved straightforwardly for the XaS and yaS by standard numerical 
techniques (e.g. Broyden's algorithm [188] in the form given in Ref. [189]), 
using the set TZ = {{sj2a-ij^j2a)y^jo} as "initial solution". Note, however, that 
the choice of initial pairings in TZ is crucial: an incorrect choice results in the 
failure of solutions with real Xa, ya to exist beyond a certain value. Since 
making this choice requires knowledge of which roots end up as conjugate 
pairs, some trial and error may be involved in finding the correct TZ (see, e.g., 
Ref. 15 of [39]). 
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B.3 Correlation functions 



In Ref . [50] , Richardson derived the following explicit results for the normal- 
ization constant J\f of (83) and the occupation probabilities Vj and correlators 
Cij of Eq. (49): 



A^=[det(M)]-'/% 
where we used the following abbreviations: 



-g\dEjdg), 



(B.23) 
(B.24) 

(B.25) 



M, 



C^2 — 23^2 ior II — 1/ 

Ejfj, = {2ej - E^)~^ , E^i, = {En - E^y^ , 

n\iJ, n\/u 



Cij,2 = X] -^i ' ^i^ = Yl ' 



Iii'{ij) = S'„_i 



1 2E12E112' 

1 D22'{ij)-2S2S2' 



1 



2En2En'2' 



2-Eln-£'l'n' 
2E2nE2'n' 

Dnn'{}j) ~ 2SnS„ 



n\ij,,u 

Df,u = Ejf^Ejy = -4:El^ + 2E,inEni, , 



Here M is a matrix with matrix indices v. The notation Y.^'' means a sum 
over all = 1, . . . , n, excluding v. The coefficients a^j, in Eq. (B.24) can most 
conveniently be found by solving the set of n algebraic equations 



(B.26) 

(B.27) 
(B.28) 



(B.29) 

(B.30) 
(B.31) 



n\^l 

(C^2 - 2Sn2)an + Y El^ar, = 1 , (a^ = 1, . . . , n) , (B.32) 

which follow from differentiating Eq. (84) with respect to g. In Eq. (B.29), S'^_i 
is an operator which symmetrizes the (n — 1) primed indices r}' = 2', . . . ,n' 
occuring in the determinant for 7ii/(ij), and then sets rj' = ri; in other words, 
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S' stands for the sum over the (n — 1)! permutations of these indices, and for 
f]' = 2', . . . n', the prime on f]' is only a mark to distinguish it from r] for the 
purposes of symmetrization. (To obtain simply make the replacements 

!<->// and 1' <-> in the determinant written in Eq. (B.29) for then 
symmetrize over the indices r]' — 1\ . . .n\ excluding //'.) The need to evaluate 
(n— 1)! permutations makes Cij practically impossible to calculate numerically 
for n much larger than, say, 10. However, Richardson showed that Eq. (B.25) 
simplifies to the (tractable) 

C,,-~^a^£;,^£;,-^, (B.33) 

IJ,V 

if D^v{ij) is approximated by D^^ in Eq. (B.29). The accuracy of this approx- 
imation can be jugded by using the fact [51] that when used to evaluate £n — 
{^n\Hu\^n)^ it produces, instead of zero, the result {l/g) YZ ^Ad"^^^ 
1], which can easily be evaluated numerically. 
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